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OF  THE  AQUIA  AND  PINEY  POINT-NANJEMOY 
AQUIFER  SYSTEM  IN  SOUTHERN  MARYLAND 

by 

Francis  H.  Chapelle  and  David  D.  Drummond 


ABSTRACT 


The  hydrogeology  and  ground-water  geochemistry  of  the  Aquia  and  Piney  Point-Nanjemoy  aquifers  in 
Southern  Maryland  have  been  investigated.  This  study  was  made  in  order  to  evaluate  the  availability  and 
chemical  quality  of  water  from  this  aquifer  system. 

The  Aquia  aquifer  is  Paleocene  in  age  and  is  a  medium-  to  fine-  grained  quartz  sand.  It  is  highly  glauconitic 
and  contains  abundant  carbonate  shell  material.  The  transmissivity  of  the  Aquia  ranges  from  200  to  2,000 
feet  squared  per  day  and  tends  to  increase  to  the  northeast  along  strike.  The  Aquia  is  overlain  by  the  Marlboro 
Clay  and  lower  Nanjemoy  Formation  which  together  act  as  a  confining  bed.  This  confining  bed  ranges  from 
100  to  250  feet  thick  and  exhibits  vertical  hydraulic  conductivities  that  range  from  10-7  to  10~10  feet  per  sec¬ 
ond.  In  1980,  approximately  6  million  gallons  of  water  per  day  was  produced  from  the  Aquia  aquifer  in 
Southern  Maryland.  Several  cones  of  depression  ranging  from  20  to  80  feet  below  sea  level  have  developed  in 
response  to  this  pumping  stress. 

The  Piney  Point-Nanjemoy  aquifer  is  Eocene  in  age  and  is  a  coarse-  to  fine-grained  glauconitic  quartz  sand 
which  contains  abundant  carbonate  shell  material.  The  transmissivity  of  the  Piney  Point-Nanjemoy  aquifer 
ranges  from  100  to  500  feet  squared  per  day.  The  Piney  Point-Nanjemoy  is  overlain  by  Chesapeake  Group 
sediments  which  act  as  a  confining  bed.  This  confining  bed  ranges  from  150  to  250  feet  thick  and  exhibits  ver¬ 
tical  hydraulic  conductivities  which  range  from  10-7  to  10~10  feet  per  second.  In  1980,  approximately  2  million 
gallons  of  water  per  day  was  produced  from  the  Piney  Point-Nanjemoy  aquifer  in  Southern  Maryland. 

Hydrogeologic  information  including  water  levels,  altitude  of  top,  thickness,  and  transmissivity  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifer  system  is  presented  as  a  series  of  maps.  Hydrogeologic  information 
such  as  thickness,  vertical  hydraulic  conductivity,  and  specific  storage  of  confining  beds  is  also  presented. 
This  information  combined  with  geologic  cross  sections  obtained  during  bridge  construction  across  the 
Chesapeake  Bay  and  the  Patuxent  River  demonstrate  that  Pleistocene  river  channels  have  in  places  truncated 
aquifers  and  confining  beds.  This  erosional  truncation  is  a  major  control  on  the  natural  head  distribution  of 
the  Aquia  and  Piney  Point-Nanjemoy  aquifer  system. 

A  quasi  three-dimensional  digital  model  of  this  aquifer  system  was  constructed.  The  Aquia  and  Piney  Point- 
Nanjemoy  aquifers  were  modeled  as  confined  aquifers  separated  by  semipermeable  confining  material. 
Recharge  to  these  aquifers  was  considered  to  occur  by  leakage  from  the  overlying  Pleistocene  water-table 
aquifer.  The  pumping  history  of  the  aquifer  system  was  simulated  from  1890  to  1980  during  model  calibra¬ 
tion.  The  model  was  calibrated  by  adjusting  the  vertical  hydraulic  conductivity  of  confining  beds  and  match¬ 
ing  calculated  water  levels  to  measured  water  levels.  Recharge  rates  to  the  aquifer  system  during  the 
simulated  pumping  history  varied  from  0.22  to  0.57  inches  per  year.  Simulations  with  the  calibrated  model  in¬ 
dicate  that  the  aquifer  system  can  sustain  the  1980  rates  of  pumpage  to  the  year  2000  with  water-level 
declines  of  about  0.5  feet  per  year  in  the  Aquia  aquifer.  Simulations  that  increase  pumpage  according  to  pro¬ 
jected  population  growth  predict  water-level  declines  on  the  order  of  1.5  feet  per  year  in  the  Aquia  aquifer  to 
the  year  2000.  A  simulation  which  assumes  that  all  holders  of  Gound-water  Appropriation  Permits  pump 
their  maximum  allocation  for  10  years  predicts  water-level  declines  of  8  feet  per  year  near  Lexington  Park  and 
Piney  Point  in  southern  St.  Marys  County,  Maryland. 

The  principles  of  equilibrium  chemistry  were  utilized  to  investigate  the  ground-water  geochemistry  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers.  Trends  of  calcium,  magnesium,  sodium,  and  bicarbonate  concen¬ 
trations  plotted  versus  distance  along  the  flowpath  delineate  three  regions  in  the  aquifers  where  different  pat¬ 
terns  of  concentration  changes  occur.  Chemical  models  were  constructed  and  tested  for  each  region  to  account 
for  the  observed  changes  in  water  chemistry.  Mass  balance  calculations  based  on  the  verified  models  suggest 
that  observed  calcite  cementation  of  the  Aquia  aquifer  is  post-depositional  and  has  resulted  from  the  reaction 
of  ground  water  with  aquifer  material. 


1 


INTRODUCTION 


Purpose  and  Scope 

An  adequate  supply  of  fresh  water  for  use  in 
households,  industry,  and  agriculture  is  an  impor¬ 
tant  factor  in  the  economic  well  being  of  all  com¬ 
munities.  Many  regions  in  the  United  States  rely  on 
surface-water  reservoirs  such  as  natural  lakes,  ar¬ 
tificial  lakes,  or  rivers  for  fresh-water  supplies.  In 
Southern  Maryland,  the  extensive  use  of  surface 
water  is  not  feasible  due  to  the  brackish  nature  of 
most  surface-water  bodies.  Southern  Maryland, 
therefore,  has  come  to  rely  on  ground  water  for  most 
of  its  water  needs.  The  Aquia  and  Piney  Point- 
Nanjemoy  aquifers  of  Paleocene  and  Eocene  age  re¬ 
spectively  are  two  major  sources  of  fresh  ground 
water.  The  Aquia  aquifer  is  used  in  conjunction  with 
deeper  Cretaceous  aquifers  in  Anne  Arundel,  Prince 
Georges,  and  Charles  Counties.  In  St.  Marys  and 
Calvert  Counties,  the  Aquia  and  Piney  Point- 
Nanjemoy  aquifers  are  the  principal  sources  of 
potable  water. 

Since  the  first  artesian  wells  were  drilled  in 
Southern  Maryland  at  the  end  of  the  19th  century, 


the  use  of  ground  water  has  increased  rapidly.  This 
increased  pumpage  has  caused  widespread  declines 
of  water  levels  in  many  aquifers.  Proper  manage¬ 
ment  of  ground-water  resources  in  future  years  will 
require  an  adequate  understanding  of  available 
aquifers.  The  purpose  of  this  report  is  to  provide 
county  planners,  elected  officials,  builders,  contrac¬ 
tors,  engineers,  well  drillers,  and  individual  citizens 
with  information  on  the  availability  and  chemical 
quality  of  ground  water  from  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers. 

The  first  section  of  this  report  deals  with  the 
physical  geology  and  hydrogeology  of  the  Aquia 
and  Piney  Point-Nanjemoy  aquifer  system.  Infor¬ 
mation  on  aquifer  and  confining-bed  stratigraphic 
relationships,  thickness,  depth,  areal  distribution, 
lithology,  and  water  levels  has  been  assembled  and 
presented.  This  information  will  be  useful  to  engi¬ 
neers,  geologists,  and  well  drillers  interested  in  the 
practical  aspects  of  siting  and  constructing  wells. 
This  section  includes  information  on  the  Aquia  and 
Piney  Point-Nanjemoy  aquifers  which  has  previous¬ 
ly  been  scattered  among  many  reports. 


Figure  1.— Location  of  study  area  and  boundary  of  modeled  area. 
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The  second  section  describes  a  quantitative  eval¬ 
uation  of  the  future  availability  of  water  from  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers  in  South¬ 
ern  Maryland,  using  a  three-dimensional  digital 
model.  This  technique  was  chosen  because  it  has 
been  shown  that  a  properly  calibrated  digital  model 
can  predict  the  water-level  response  of  an  aquifer 
system  to  expected  rates  of  future  pumpage  with  a 
fair  degree  of  accuracy.  These  predictions  can  be  of 
great  value  to  water  planners  in  deciding  on  future 
sources  of  water. 

The  third  section  describes  the  ground-water  geo¬ 
chemistry  of  the  Aquia  and  Piney  Point-Nanjemoy 
aquifers.  The  chemical  quality  of  ground  water  is  an 
important  consideration  in  the  management  of  an 
aquifer  system.  This  chemical  quality  results  from 
the  chemical  reaction  of  recharge  water  with  gases 
and  minerals  in  the  aquifer.  An  understanding  of 
these  chemical  processes  will  aid  in  the  management 
of  the  aquifers  and  contribute  to  an  understanding 
of  the  system  as  a  whole. 

Location  and  Extent  of  the  Study  Area 

Southern  Maryland  lies  within  the  Atlantic 
Coastal  Plain  physiographic  province.  This  pro¬ 
vince,  which  extends  from  South  Carolina  to  New 
York,  is  characterized  by  generally  low  topography. 
The  Atlantic  Coastal  Plain  is  underlain  by 
predominantly  unconsolidated  clastic  sediments  of 
Lower  Cretaceous  to  recent  age,  which  thicken  to 
the  southeast. 

The  project  area  is  shown  in  figure  1.  The  area  in¬ 
cludes  the  southern  half  of  Anne  Arundel  County 
and  all  of  Calvert  and  St.  Marys  Counties.  The  boun¬ 
dary  of  the  digital  model  is  also  shown  in  figure  1. 
The  area  modeled  was  considerably  larger  than  the 
project  area  in  order  to  simulate  the  lateral  extent 
and  distant  boundaries  of  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers.  Included  in  the  model  are 
Anne  Arundel,  Prince  Georges,  Charles,  St.  Marys, 
Calvert,  Dorchester,  Talbot,  Caroline  and  Queen  An- 
nes  Counties. 

Methods  of  Investigation 

The  first  phase  of  this  study  was  to  define  the 
hydrogeologic  framework  of  the  Aquia  and  Piney 
Point-Nanjemoy  aquifer  system  in  Southern  Mary¬ 
land.  Data  were  collected  for  the  topography  and 
stratigraphy  of  Southern  Maryland  and  the  Eastern 
Shore;  depth  to  top,  thickness,  physical  boundaries, 
and  historical  water  levels  of  the  aquifers;  and  areal 
distribution  and  thickness  of  the  confining  beds. 

Some  of  this  information  was  available  from 
previously  published  reports.  However,  much  had  to 


be  gathered  in  the  field.  This  fieldwork  consisted  of: 

(1)  Establishing  an  observation-well  network  in 
Southern  Maryland  and  on  the  Eastern  Shore 
for  both  the  Aquia  and  Piney  Point-Nanjemoy 
aquifers.  This  network  included  a  total  of  160 
wells.  Fifteen  observation  wells  were  drilled 
for  this  project  where  water  levels  were  not 
otherwise  available. 

(2)  Obtaining  geologists’  logs  and  geophysical 
logs  of  15  new  observation  wells. 

(3)  Obtaining  undisturbed  cores  of  confining-bed 
material  for  laboratory  determination  of 
confining-bed  vertical  hydraulic  conductivity. 

(4)  Obtaining  aquifer-test  data. 

(5)  Obtaining  pumpage  information  from  users  of 
greater  than  5,000  gal/d. 

The  second  phase  was  to  construct  and  calibrate  a 
three-dimensional  digital  model  of  the  Aquia  and 
Piney  Point-Nanjemoy  aquifer  system  in  Southern 
Maryland.  This  was  done  by: 

(1)  Developing  a  finite-difference  grid. 

(2)  Collecting  information  on  the  areal  extent, 
transmissivity,  and  storage  coefficient  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers. 

(3)  Translating  potentiometric,  transmissivity, 
and  confining-bed  thickness  maps  into 
computer-manipulatable  data  sets. 

(4)  Calibrating  the  digital  model  using  historical 
pumpage  and  water  levels. 

The  calibrated  model  was  then  utilized  to  predict 
water-level  changes  resulting  from  several  senarios 
of  future  pumpage.  These  pumping  senarios  are 
designed  to  aid  water  planners  in  making  water- 
management  decisions. 

The  third  phase  was  to  describe  the  ground-water 
geochemistry  of  the  Aquia  and  Piney  Point- 
Nanjemoy  aquifers.  Methods  used  included: 

(1)  Collecting  available  chemical  data  on  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers 
from  published  sources. 

(2)  Graphing  sodium  (Na+),  calcium  (Ca2+)  and 
magnesium  (Mg2+),  and  bicarbonate  (HC03) 
concentrations  versus  distance  along  the 
flowpath.  These  plots  illustrate  water  chemis¬ 
try  changes  as  water  moves  down  the  hydro- 
logic  gradient. 

(3)  Constructing  chemical  models  based  on  simp¬ 
lified  aquifer  lithology  and  idealized  mineral 
composition. 

(4)  Testing  the  chemical  models  by  matching 
model  predictions  to  observed  data. 
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Previous  Investigations 

The  Paleocene-Eocene  aquifers  in  Maryland  were 
first  described  by  Darton  (1891),  who  considered 
them  to  be  a  thick  undivided  unit  called  the 
“Pamunkey  Formation.”  The  Pamunkey  was  later 
elevated  to  Group  status  and  subdivided  to  include 
the  Aquia  Formation,  the  Nanjemoy  Formation  and 
the  Marlboro  Clay  Formation  (Clark  and  Martin, 
1901).  The  water  resources  of  these  formations  were 
discussed  in  detail  for  each  county  in  Southern 
Maryland  by  Dryden  and  Overbeck  (1948),  Bennion 
and  Brookhart  (1949),  Cook  and  others  (1952),  and 
Martin  and  Ferguson  (1953).  The  ground-water 
resources  of  Southern  Maryland  were  described  by 
Otton  (1955),  who  first  applied  the  name  “Piney 
Point  Formation”  to  the  coarse  sands  of  Eocene  age 
that  overly  the  Nanjemoy  Formation  in  Southern 
Maryland.  Back  (1966)  described  the  regional 
ground- water  geochemistry  of  the  Aquia  and  Nan¬ 
jemoy  aquifers.  Glaser  (1968,1971)  provided  a  com¬ 
prehensive  discussion  of  the  geology  of  the  Coastal 
Plain  sediments  in  Southern  Maryland.  Weigle  and 
others  (1970)  constructed  several  hydrologic  maps 
of  Southern  Maryland.  Mack  and  others  (1971) 
discussed  the  Aquia  and  Piney  Point  aquifers  on  the 
Eastern  Shore  of  Maryland.  Hansen  (1972)  con¬ 
tributed  a  guide  to  the  aquifers  of  Maryland’s 
Coastal  Plain  with  emphasis  on  the  practical  use  of 
these  aquifers.  Hansen  (1974)  interpreted  the 
depositional  environments  of  the  Aquia  Formation. 
Kapple  and  Hansen  (1976)  constructed  a  two- 
dimensional  digital  flow  model  of  the  Aquia  Forma¬ 
tion  in  Southern  Maryland.  Mack  (1976)  collected 
geohydrologic  data  on  the  Aquia  Formation  from 
several  wells  drilled  at  Chalk  Point  in  Prince 
Georges  County.  Lucas  (1976)  collected  well  data  for 


Anne  Arundel  County.  Similar  data  were  compiled 
for  Calvert  and  St.  Marys  Counties  by  Weigle  and 
Webb  (1970)  and  Drummond  (in  preparation).  Han¬ 
sen  (1977)  catalogued  data  from  two  core  holes  in 
Prince  Georges  County  and  Queen  Annes  County. 
Woll  (1978)  provided  many  chemical  analyses  of 
water  for  Maryland.  Williams  (1979)  constructed  a 
two-dimensional  digital  fow  model  of  the  Piney 
Point  aquifer  in  Southern  Maryland  and  on  the 
Eastern  Shore  of  Maryland. 
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GEOLOGY 


Regional  Geologic  and  Stratigraphic  Framework 

The  sediments  that  make  up  the  Coastal  Plain 
physiographic  province  in  Southern  Maryland  con¬ 
sist  of  unconsolidated  clay,  silt,  sand,  and  gravel 
beds,  which  dip  gently  to  the  southeast.  These 
sediments  crop  out  in  a  concentric  band  that  lies 
parallel  to  the  Fall  Line.  The  Fall  Line  marks  the 
western  boundary  of  the  Coastal  Plain.  The  Coastal 
Plain  sediments  are  underlain  by  Precambrian  and 
Paleozoic  gneiss,  schist,  and  gabbroic  rocks,  which 
are  usually  referred  to  as  “basement”  rocks.  The 


surface  of  the  basement  rocks  underlying  Southern 
Maryland  has  been  downwarped  into  a  structure 
termed  the  “Salisbury  Embayment”  (Richards, 
1948). 

The  Salisbury  Embayment  is  part  of  a  much 
larger  regional  basement  structure  called  the 
Chesapeake-Delaware  Embayment  (Murray,  1961). 
This  structure,  which  is  prominent  in  the  basement 
rocks,  loses  form  in  the  younger  Coastal  Plain 
sediments.  Southern  Maryland  occupies  a  south- 
central  position  on  the  Chesapeake-Delaware  Em¬ 
bayment. 
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Table  1.— Generalized  stratigraphy  of  Southern  Maryland 


System 

Series 

Stratigraphic 

unit 

Thickness 
(  feet) 

Dominant  lithologic  character 

Water-bearing  properties 

Quaternary 

and 

Tertiary  (?) 

Holocene 
Pleistocene 
Pliocene  (?) 

Lowland 

and 

upland 

deposits 

0-190 

Sand,  gravel,  and  silt;  tan  to  rusty 
orange;  predominantly  quartz. 

Yields  small  to  moderate  amounts  of  water 
to  wells.  Utilized  primarily  as  a  water 
source  for  shallow,  domestic,  and  farm 
wells.  The  upper  recharging  water- 
table  aquifer  to  the  Aquia  and  Piney 
Point-Nan jemoy  aquifers  in  southern 
Maryland . 

Tertiary 

Miocene 

Chesapeake  Group 

St.  Marys 
Formation 

0-80 

Clay,  sandy,  silty;  greenish-blue  to 
yellow;  fossiliferous . 

Functions  as  a  confining  bed. 

Chop tank 

Formation 

0-60 

Clay,  silty;  olive-green  to  gray; 
fossiliferous . 

Functions  as  a  confining  bed. 

Calvert 

Formation 

0-180 

Clay,  silty;  olive-green  to  gray; 
fossiliferous;  lower  member  is 
diatomaceous  and  contains  phos- 
phatic  pebbles. 

Functions  as  a  confining  bed. 

Eocene 

Piney  Point 
Formation 

0-80 

Sand;  grayish-green  to  grayish-white; 
medium-  to  coarse-grained;  quartz 
is  most  common  mineral;  glauconitic; 
calcite-cemented  shell  beds  common. 

Important  source  of  water  in  southern 

Calvert  and  St.  Marys  Counties.  Hy¬ 
draulically  connected  to  the  upper 
sandy  portion  of  the  Nanjemoy  Formation. 

Nan jemoy 

Formation 

0-250 

Sand,  silt,  clay;  blackish-green  to 
gray;  quartz  most  common  mineral; 
glauconitic;  the  upper  portion  of 
formation  is  predominantly  sand; 
the  lower  portion  is  predominantly 
silt  and  clay. 

The  upper  sandy  portion  is  an  important 
source  of  water  in  Calvert  and  St. 

Marys  Counties  and  is  hydraulically 
connected  to  the  overlying  Piney 

Point  Formation.  The  lower  portion  of 
the  formation  functions  as  a  confining 
bed . 

Marlboro  Clay 
Formation 

0-35 

Clay,  pinkish-red  to  silvery-gray; 
very  plastic;  thin  lenses  of 
pale  gray  silt. 

Functions  as  a  confining  bed. 

Paleocene 

Aquia 

Format  ion 

0-230 

Sand,  greenish-black;  quartz  most 
common  mineral;  glauconitic; 
lenses  of  silty-clay  and  shell 
beds  common;  calcite  cementation 
common . 

A  primary  source  of  water  in  southern  Anne 
Arundel  County  and  in  St.  Marys  and 
Calvert  Counties.  An  important  source  of 
water  in  southern  Charles  and  Prince 
Georges  Counties. 

Br ightseat 
Formation 

0-40 

Silt,  clayey;  gray  to  dark  gray; 
micaceous . 

Functions  as  a  confining  bed. 

Cretaceous 

Upper 

Cretaceous 

Severn 

and 

Magothy 

Formations , 
undifferentiated 

350-1,700 

Silt,  sand,  clay,  interbedded. 

Magothy  is  a  primary  water  source  in  north¬ 
ern  Anne  Arundel  County,  but  is  not 
present  in  southern  Calvert  or  St.  Marys 
Counties.  Severn  and  Matawan  function 

as  confining  beds. 

Lower 

Cretaceous 

Patapsco, 

Arundel,  and 
Patuxent 

Formations , 
undifferentiated 

300-2,500 

Silt,  sand,  clay,  gravel  inter¬ 
bedded  . 

Yields  large  amounts  of  water  in  Charles, 
Prince  Georges,  and  northern  Anne  Arundel 
Counties.  Untested  in  southern  Calvert 
and  St.  Marys  Counties. 

Paleozoic 

and 

Precambrian 

Crystalline 

rocks 

(basement) 

Unknown 

Complex  assemblage  of  schists, 
granites,  gneisses,  and 
gabbros . 

Untested . 

The  wedge  of  sediments  that  make  up  the  Coastal 
Plain  of  Maryland  consists  of  beds  ranging  from 
Jurassic(?)  to  Holocene  in  age.  This  report  is 
directed  primarily  to  those  sediments  that  range  in 
age  from  Paleocene  to  Quaternary. 

Plates  1  and  2  show  cross  sections  of  the  idealized 
lithology,  thickness,  and  distribution  of  the  major 
Tertiary  and  Quaternary  hydrogeologic  units  in 
Southern  Maryland.  The  generalized  stratigraphy  of 
Southern  Maryland  is  shown  in  table  1. 


Topography 

Southern  Maryland  is  characterized  by  hilly,  ter¬ 
raced  upland  areas  which  contrast  markedly  with 
the  low  flat-lying  areas  of  Maryland’s  Eastern 
Shore.  Where  upland  areas  have  been  dissected  by 
streams,  valley  walls  are  often  steep  and  local  relief 
of  100  ft  is  common.  The  valleys  of  major  streams, 
such  as  the  Patuxent  and  Potomac  Rivers,  are  char¬ 
acterized  by  a  system  of  well-marked  gravel  and 
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loam-covered  terraces  which  have  been  studied  by 
Hack  (1955).  As  these  major  rivers  approach  the 
Chesapeake  Bay,  their  valleys  tend  to  become 
broader  and  there  is  a  general  lowering  of  top¬ 
ographic  relief. 

The  distinctive  upland  topography  of  Southern 
Maryland  has  developed  in  response  to  the  repeated 
steepening  and  flattening  of  stream  gradients  which 
has  accompanied  fluctuations  in  sea  level  (Glaser, 
1968).  Hack  (1957)  has  shown  that  a  river  system  ex¬ 
isted  during  Pleistocene  time  which  approached 
equilibrium  with  a  sea  level  300  to  400  ft  below  the 
present  sea  level.  At  the  close  of  Pleistocene  time, 
the  rising  sea  submerged  much  of  this  river  system 
creating  the  Chesapeake  Bay  estuary. 


Aquia  Formation 
Lithology 

The  Aquia  Formation  in  the  subsurface  consists 
of  a  medium-to  coarse-grained,  medium-to  well- 
sorted  glauconitic-quartz  sand.  Carbonate  shell 
debris  is  abundant  and,  in  places,  makes  up  approx¬ 
imately  20  percent  of  the  aquifer  material.  The  usual 
percentage  of  shell  material,  however,  is  between  1 
and  5  percent.  Quartz  sand  grains,  which  make  up 
50  to  75  percent  of  the  aquifer  material,  are  often 
angular  to  subangular  and  frequently  exhibit 
goethite  staining.  Some  grains  are  also  character¬ 
ized  by  a  high  polish.  Glauconite,  a  hydrous  iron¬ 
bearing  aluminum  silicate,  typically  makes  up  20  to 
40  percent  of  the  aquifer  material.  Glauconite  occurs 
in  sand-sized  grains  that  are  usually  fine  to  medium 
in  size.  Subspherical  grains  or  botryoidal  masses  of 
subspherical  grains  are  also  commonly  observed. 

Although  the  Aquia  is  predominantly  sand,  some 
zones  occur  which  contain  significant  quantities  of 
clay  and  silt-sized  particles.  X-ray  defraction 
analyses  of  this  fine-grained  fraction  show  that  the 
clay  mineralogy  is  mostly  of  the  mixed-layer  type 
(Hansen,  1977). 

The  heavy  mineral  suite  identified  in  the  Aquia  in¬ 
cludes  hornblende,  tourmaline,  epidote,  andalusite, 
sillimanite,  garnet,  zircon,  ilmenite,  and  magnetite. 
Ilmenite  and  magnetite  are  the  most  commonly  ob¬ 
served  heavy  minerals. 

The  lithologic  character  of  the  Aquia  Formation 
in  outcrop  differs  markedly  from  that  of  the  subsur¬ 
face.  In  outcrop,  the  Aquia  appears  rusty  or  orange- 
red  as  opposed  to  a  dark  green  in  the  subsurface. 
Shell  material  in  outcrops  is  severely  leached  and 
ledges  of  goethite-cemented  sand  are  common. 
Wolff  (1967)  presents  a  detailed  X-ray  analysis  of 


these  goethite  ledges  and  nodules  in  Anne  Arundel 
County  near  Annapolis,  Md.,  and  concludes  that 
they  result  from  the  weathering  of  glauconite. 

Stratigraphy 

The  Aquia  Formation  crops  out  in  a  band  extend¬ 
ing  from  the  Potomac  River  near  Washington,  D.C., 
through  Anne  Arundel  County  near  the  Chesapeake 
Bay  Bridge.  This  outcrop  (fig.  4),  which  is  often 
mantled  by  younger  Pleistocene  sediments,  ranges 
in  width  from  about  2  mi  near  the  Potomac  River  to 
9  mi  in  southern  Prince  Georges  County. 

The  Aquia  Formation  has  been  interpreted  by 
Glaser  (1968)  and  Hansen  (1974)  as  having  been 
deposited  in  a  regressive  marine  environment.  The 
basal  sediments  are  typically  fine-grained  sug¬ 
gesting  deposition  in  the  sublittoral  zone  (about  300 
ft  of  water).  The  upper  sections  are  considerably 
coarser  reflecting  sand  deposition  in  a  higher- 
energy,  shallower  environment. 

The  Aquia  Formation  near  the  Potomac  River  un- 
conformably  overlies  the  Upper  Cretaceous  Severn 
Formation.  To  the  east,  however,  it  overlies  the 
Brightseat  Formation  of  Paleocene  age  either  con¬ 
formably  or  with  a  minor  unconformity  (Glaser, 
1968).  The  Aquia  is  overlain  conformably  through¬ 
out  Southern  Maryland  by  the  Marlboro  Clay  For¬ 
mation. 

In  outcrop,  the  thickness  of  the  Aquia  ranges 
from  90  ft  near  the  Potomac  River  to  approximately 
150  ft  in  Anne  Arundel  County.  The  subsurface 
thickness  of  the  Aquia  as  defined  by  well  logs  shows 
a  similar  pattern  of  increasing  thickness  to  the 
northeast  and  thinning  to  the  southwest.  The 
thickness  of  the  Aquia  Formation  in  Southern 
Maryland  is  shown  in  figure  2. 

Hansen  (1974)  identified  three  sedimentary  facies 
in  the  subsurface  of  the  Aquia  on  the  basis  of  coarse 
and  medium  sand  percentage  of  the  formation. 
These  facies  are  shown  in  figure  3.  Facies  1  is  a 
coarse  sand  that  extends  from  Kent  County  to  the 
Potomac  River.  Facies  2  is  a  fine  silty  sand  which 
occurs  primarily  in  Charles  and  eastern  Prince 
Georges  County.  Facies  3  is  a  thin  predominantly 
clay-silt  facies  that  underlies  the  extreme  southern 
tip  of  St.  Marys  County  and  extends  to  the  north¬ 
west  across  the  Chesapeake  Bay.  It  is  not  clear 
which  facies  of  the  Aquia  are  present  across  the 
Potomac  River  in  the  Northern  Neck  Peninsula  of 
Virginia.  However,  the  best  available  information 
(Teifke,  1973)  suggests  that  the  Aquia  in  eastern 
part  of  Northern  Neck  is  the  equivalent  of  facies  3. 
In  the  vicinity  of  the  Aquia  Formation  type  section 
in  Virginia,  the  formation  is  a  fine,  poorly  sorted, 
silty  sand  (Drobnyk,  1965). 
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The  subsurface  relationships  of  the  Aquia  Forma¬ 
tion  to  other  Tertiary  and  Quaternary  sediments  in 
Southern  Maryland  are  shown  on  plates  1  and  2. 

Structure 

The  subsurface  structure  of  the  Aquia  Formation 
in  Maryland  is  shown  in  figure  4.  The  distinctive 
bend  along  a  northwest-trending  axis  in  southern 
Calvert  County  has  been  interpreted  by  Hansen 
(1974)  as  reflecting  the  deeper  basinal  patterns  of 
the  Chesapeake-Delaware  Embayment.  The  dip  of 
the  Aquia  ranges  from  approximately  10  ft/mi  in  St. 
Marys  County  to  20  ft/mi  in  Queen  Annes  County. 

Hansen  (1974)  has  pointed  out  an  oblique  struc¬ 
tural  relationship  between  the  outcrop  belt  and 
several  lithofacies  trends  in  the  Aquia.  He  suggests 
that  the  present  structural  strike  of  the  Aquia  For¬ 
mation  is  different  from  the  structural  strike  during 
Aquia  deposition.  This  may  indicate  some  slight 
post-depositional  tilting. 

Marlboro  Clay 

The  Marlboro  Clay,  which  overlies  the  Aquia  For¬ 
mation,  is  a  pink  to  silver-gray  plastic  clay.  X-ray 
analyses  show  that  it  is  composed  predominantly  of 
kaolinite  and  mixed-layer  clays  with  minor  amounts 
of  montmorillonite  (Hansen,  1977).  Thin  lenses  of 
silt  generally  less  than  an  inch  thick  are  common  in 
the  Marlboro.  Lignite  and  thin  lenses  of  glauconite 
sand  have  been  reported  from  some  localities 
(Glaser,  1968). 

The  Marlboro  Clay  has  been  interpreted  as 
representing  a  very  shallow,  probably  brackish 
water  environment  which  followed  the  marine 
regression  of  Aquia  time  (Glaser,  1968).  Burrows 
and  other  shallow-water  features  are  commonly 
observed  in  the  Marlboro  and  these  are  consistent 
with  that  interpretation.  Microfossil  data  from  the 
Marlboro  indicate  that  it  is  Paleocene  in  age 
(Reinhardt  and  others,  1980). 

The  Marlboro  Clay  ranges  in  thickness  from  about 
30  ft  in  Prince  Georges  and  Anne  Arundel  Counties 
to  less  than  5  ft  in  parts  of  southern  St.  Marys  Coun¬ 
ty  and  southern  Charles  County.  The  Marlboro  Clay 
has  not  been  identified  on  the  Eastern  Shore  of 
Maryland. 

Darton  (1948),  who  originally  named  the  Marlboro 
Clay,  assigned  it  as  a  member  of  the  Nanjemoy  For¬ 
mation.  Glaser  (1968)  pointed  out  that  the  Marlboro 
is  lithologically  distinct  from  overlying  and  underly¬ 
ing  sediments,  that  it  is  thick  enough  to  be  mapped, 
and  that  it  is  areally  extensive.  For  these  reasons, 
Glaser  (1968)  proposed  formational  status  for  the 


Marlboro  Clay.  Reinhardt  and  others  (1980),  in 
agreement  with  Glaser  (1968),  have  also  given  for¬ 
mational  status  to  the  Marlboro  Clay.  This  report, 
therefore,  will  follow  that  convention. 

The  subsurface  relationships  of  the  Marlboro  Clay 
Formation  with  other  Tertiary  and  Quaternary  sedi¬ 
ments  in  Southern  Maryland  are  shown  on  plates  1 
and  2. 

Nanjemoy  Formation 

The  Nanjemoy  Formation  overlies  the  Marlboro 
Clay  throughout  Southern  Maryland.  The  lower 
part  of  the  Nanjemoy  is  an  olive-green  glauconitic 
silty  clay.  Some  thin  beds  of  glauconitic  quartz  sand 
1  to  5  ft  thick  occur  in  the  lower  Nanjemoy,  but  they 
do  not  apear  to  be  laterally  extensive.  The  Nan¬ 
jemoy  Formation  tends  to  coarsen  upward.  The  up¬ 
per  part  of  the  Nanjemoy  is  commonly  a  glauconitic 
medium-sorted  fine  quartz  sand.  Glauconite  typical¬ 
ly  makes  up  30  to  50  percent  of  the  formation  in  the 
upper  sandy  portion.  The  upper  sandy  portion  of  the 
formation  tends  to  become  thinner  toward  the 
north.  In  southern  Anne  Arundel  County,  sand  beds 
cannot  be  distinguished  in  the  Nanjemoy  Formation 
by  geophysical  logs.  It  is  possible  that  this  updip 
thinning  of  the  upper  sandy  portion  is  due  to  a  par¬ 
tial  erosional  truncation  during  the  Oligocene 
withdrawal  of  the  sea. 

Hansen  (1972)  has  interpreted  the  Nanjemoy  as 
representing  regressive  marine  sedimentation  dur¬ 
ing  Eocene  time.  According  to  this  interpretation, 
the  lower  clay  and  silt  beds  were  deposited  in  deep 
water.  The  upper  sandy  portion  was  deposited  in  a 
progressively  shallower,  higher  energy  environ¬ 
ment. 

The  Nanjemoy  Formation  ranges  in  thickness 
from  about  100  ft  in  Charles  County  to  about  250  ft 
in  northern  St.  Marys  County  and  central  Calvert 
County. 

The  subsurface  relationship  of  the  Nanjemoy  For¬ 
mation  with  other  Tertiary  and  Quaternary 
sediments  is  shown  on  plates  1  and  2. 

Piney  Point  Formation 

The  Piney  Point  Formation  conformably  overlies 
the  Nanjemoy  Formation  in  much  of  Southern 
Maryland.  The  Piney  Point  is  grayish-white  medium 
to  coarse  slightly  glauconitic  quartz  sand.  The 
Piney  Point  is  readily  distinguished  from  the 
underlying  sand  beds  of  the  Nanjemoy  Formation 
by  its  relative  coarseness  and  by  its  much  lower  con¬ 
tent  of  glauconite.  Glauconite  in  the  Piney  Point 
generally  makes  up  less  than  5  percent  of  the  forma- 
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Figure  4.— Altitude  of  the  top  and  outcrop/subcrop  area  of  the  Aquia  Formation. 


11 


tion.  Beds  of  calcite-cemented  sand  and  shell  beds 
are  very  common  in  the  Piney  Point  Formation. 

Hansen  (1972)  has  interpreted  the  Piney  Point  as 
having  been  deposited  during  a  regressive  marine 
episode,  possibly  by  longshore  currents.  Upward 
coarsening  of  the  Piney  Point  has  been  noticed  by 
many  workers  and  is  probably  due  to  shoaling  dur¬ 
ing  a  seaward  shift  of  the  strandline. 

The  Piney  Point  Formation  has  been  truncated  in 
updip  areas  by  erosion  which  occurred  during  the 
Oligocene-Miocene  withdrawal  of  the  sea.  Because 
of  this  truncation,  the  Piney  Point  Formation  does 
not  crop  out  in  Southern  Maryland.  In  the  subsur¬ 
face,  the  Piney  Point  is  not  present  northwest  of  a 
line  that  runs  from  Leonardtown  in  St.  Marys  Coun¬ 
ty  through  Prince  Frederick  in  Calvert  County. 
Subsequent  to  the  Oligocene-Miocene  disconformi- 
ty,  the  Piney  Point  was  buried  by  Chesapeake 
Group  sediments  during  the  following  Miocene 
marine  transgression. 

Chesapeake  Group 

The  Chesapeake  Group  in  ascending  order  con¬ 
sists  of  the  Calvert,  Choptank,  and  St.  Marys  For¬ 
mations  of  Miocene  Age.  The  Calvert  Formation  is  a 


fossiliferous,  slightly  sandy  greenish-gray  silty  clay. 
The  lower  part  of  the  Calvert  is  a  distinctive 
diatomaceous  earth.  Dark  gray  phosphatic  pebbles 
and  quartz  sand  have  been  reported  at  its  base  at 
some  localities.  The  Calvert  disconformably  overlies 
the  Piney  Point  in  southern  Calvert  and  St.  Marys 
County.  In  northern  Calvert  County,  where  the 
Piney  Point  Formation  has  been  truncated,  the  Cal¬ 
vert  disconformably  overlies  the  Nanjemoy  Forma¬ 
tion.  (See  pis.  1  and  2.) 

The  Choptank  Formation  is  a  gray-green  clay  with 
some  yellowish-brown  sand.  The  Choptank  is  fossili¬ 
ferous  in  places  and  is  well  exposed  along  the 
Calvert  Cliffs. 

The  St.  Marys  Formation  is  a  dense,  bluish-gray 
clay  with  some  fine  clayey  sand.  This  sand,  where 
present,  is  frequently  glauconitic.  Calcite-cemented 
sandstone  beds  are  common  in  this  unit. 

The  total  thickness  of  Chesapeake  Group  sedi¬ 
ments  varies  from  about  175  ft  to  about  275  ft.  In 
the  subsurface,  it  is  difficult  to  distinguish  the  for¬ 
mations  of  the  Chesapeake  Group  from  each  other. 
For  this  reason,  these  sediments  will  be  referred  to 
collectively  as  the  Chesapeake  Group  throughout 
this  report.  The  subsurface  relationships  of  the 
Chesapeake  Group  with  other  Tertiary  sediments 
are  shown  on  plates  1  and  2. 


HYDROGEOLOGY 


Aquia  Aquifer 
Regional  Extent 

The  Aquia  Formation  is  extensively  tapped  for 
fresh  ground  water  in  Anne  Arundel,  Calvert,  and 
St.  Marys  Counties.  It  is  used  to  a  lesser  extent  in 
Charles  and  Prince  Georges  Counties.  On  the  East¬ 
ern  Shore  of  Maryland,  the  Aquia  Formation  is 
tapped  in  Queen  Annes,  Kent,  and  Talbot  Counties 
(fig.l).  Across  the  Potomac  River  in  Virginia,  the 
Aquia  is  not  widely  used  as  a  source  of  water. 

Generally,  the  Aquia  Formation  is  most  produc¬ 
tive  in  facies  1  (fig.  3).  Facies  1  is  typically  a  thick, 
medium  to  coarse  sand  and  has  the  highest  permea¬ 
bility.  Facies  2  (fig.  3)  is  also  tapped.  However,  due 
to  the  higher  silt  content  and  the  finer-grained  sand, 
it  is  not  as  productive  as  facies  1.  Facies  3,  which 
consists  predominantly  of  silt  and  clay,  is  not  an 
aquifer. 

In  this  report,  the  term  “Aquia  aquifer”  refers  to 
those  portions  of  the  Aquia  Formation  that  can  be 
utilized  to  yield  ground  water.  This  includes  facies  1 
and  2,  and  excludes  facies  3. 


Transmissivity  and  Storage  Coefficient 

Transmissivity  is  defined  as  the  rate  at  which 
water  is  transmitted  through  a  unit  width  of  aquifer 
under  a  unit  hydraulic  gradient  (Lohman  and  others, 
1972).  Transmissivity  is  closely  related  to  the 
thickness  of  an  aquifer.  Generally,  the  thickest  part 
of  an  aquifer  is  the  most  transmissive.  However,  as 
has  been  pointed  out  by  Williams  (1979),  local 
calcite  cementation  or  other  lithologic  inhomo¬ 
geneities  make  a  direct  correlation  between  aquifer 
thickness  and  transmissivity  unfeasible.  In  prac¬ 
tice,  therefore,  the  transmissivity  of  an  aquifer  is 
often  determined  by  analyzing  pumping  test  infor¬ 
mation  using  the  Theis  non-equilibrium  equation 
(Theis,  1936).  Transmissivities  obtained  by  this 
method,  however,  are  subject  to  many  assumptions 
and  simplifications  and  should  be  considered  ap¬ 
proximations  rather  than  exact. 

Figure  5  shows  the  transmissivity  distribution  of 
the  Aquia  aquifer  as  determined  by  aquifer  tests.  A 
comparison  of  this  transmissivity  distribution  (fig. 
5)  with  aquifer  thickness  (fig.  2)  shows  that  the 
thickest  part  of  the  aquifer  is  also  the  most  trans- 
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Figure  5.— Transmissivity  of  the  Aquia  aquifer. 
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Table  2.— Vertical  hydraulic  conductivities  of  confining  beds  as  determined  by  laboratory  methods 


Geologic 

U.S.G.S. 1/ 

Depth  of 

Hydraulic 

Consolidation 

County 

unit 

No. 

sample 

conductivity 

load 

Permeant 

(ft) 

(ft/s) 

( lb/  iii  ) 

St.  Marys 

Chesapeake  Group 

SM-Dd-46— ^ 

237-237.1 

2.95X10-7 

90 

Distilled  water. 

Chesapeake  Group 

3/ 

SM-Eg-28— 

217-217.5 

1 .20X10-7 

1,000 

Simulated  forma¬ 
tion  water. 

Chesapeake  Group 

2/ 

SM-Dd-49— 

236 

8 . 5X10~10 

236 

Do . 

Marlboro  Clay 

2/ 

SM-Dd-49— 

458 

-9 

1.1X10 

458 

Do. 

Chesapeake  Group 

2/ 

SM-Dd-50- 

198 

-9 

7.3X10 

200 

Do . 

Calvert 

Chesapeake  Group 

2/ 

CA-Fd-54- 

215 

6.90X10-10 

215 

Do. 

Marlboro  Clay 

2/ 

CA-Fd-54— 

516 

6.7X10-10 

516 

Do . 

Chesapeake  Group 

2/ 

CA-Cc-57- 

197 

2.3X10-9 

197 

Do . 

Prince  Georges 

Marlboro  Clay 

4  / 

PG-Df-35- 

36-38 

3. 13X10"9 

24 

Do . 

Queen  Annes 

Chesapeake  Group 

4/  • 

OA-Bg-54— ' 

45-47 

6.95X10-10 

30 

Do . 

1/  Location  of 

test  holes  shown  on 

plates  12  and 

13. 

2/  Analysis  pe 

rformed  by  TerraTek, 

Salt  Lake  City 

,  Utah. 

3/  Analysis  performed  by  Core  Laboratories,  Inc.,  Dallas,  Texas. 

4/  Analysis  performed  by  U.S.  Geological  Survey  Hydrologic  Laboratory,  Denver,  Colorado. 


missive.  Also,  the  most  transmissive  tracts  of  the 
Aquia  aquifer  correspond  to  facies  1  (fig.  3). 

Storage  coefficient  is  defined  as  the  volume  of 
water  released  from  a  unit  surface  area  of  aquifer 
per  unit  drop  in  head  (Lohman  and  others,  1972).  In 
confined  aquifers,  water  derived  from  storage  comes 
in  part  from  expansion  of  water  and  compression  of 
aquifer  material.  The  storage  coefficient  is  usually 
determined  by  analysis  of  pumping  test  data  using 
the  Theis  non-equilibrium  equation  (Theis,  1936). 
Reported  values  of  the  storage  coefficient  in  the 
Aquia  aquifer,  as  determined  by  pumping  tests, 
range  from  0.0001  to  0.0004  (Hansen,  1972). 

Upper  and  Lower  Confining  Beds 

The  clay  and  silt  formations  in  the  Coastal  Plain 
sediments  generally  have  such  low  permeability 
that  they  cannot  be  directly  tapped  for  water.  These 
beds  are  important,  however,  because  they  control 
the  vertical  leakage  of  water  between  aquifers.  In 
addition,  whereas  clay  and  silt  beds  have  low  per¬ 
meability,  they  commonly  have  high  porosity,  so 


that  they  may  contain  significant  quantities  of 
water  in  storage. 

Leakage  through  a  confining  bed  is  controlled  by 
the  thickness  of  the  bed  and  its  hydraulic  conduc¬ 
tivity.  Hydraulic  conductivity  is  defined  as  the 
volume  of  water  that  can  be  transmitted  in  unit  time 
under  a  unit  hydraulic  gradient  (Lohman,  1972).  In 
confining-bed  materials,  which  are  several  orders  of 
magnitude  less  permeable  than  the  aquifers  they 
confine,  water  movement  is  predominantly  vertical. 
In  general,  therefore,  the  vertical  hydraulic  conduc¬ 
tivity  of  a  confining  bed  is  the  important  parameter 
when  evaluating  leakage. 

The  ability  of  a  confining  bed  to  store  water  is 
described  by  its  specific  storage.  Specific  storage  is 
defined  as  the  volume  of  water  released  from  or 
taken  into  storage  per  unit  volume  of  confining  bed 
per  unit  change  in  head  (Lohman,  1972). 

There  is  little  available  information  on  the 
hydraulic  properties  of  the  Brightseat  and  Severn 
Formations  that  underlie  the  Aquia.  One  undis¬ 
turbed  core  of  material  from  the  Brightseat  is 
available  from  Prince  Georges  County  (Hansen, 
1977).  Laboratory  analyses  of  this  material  in- 


14 


76°45' 


76”  15' 


76°  30' 


15 


Table  3.— Records  of  Aquia  and  Piney  Point-Nanjemoy  observation  wells  drilled  during  the  project. 


U.S.G.S. 

Well 

No. 

State 

Permit 

No. 

Nearest 

town 

Latitude 

Longitude 

Altitude 
of  land 
surface 
(ft) 

Date 

drilled 

Drilled  depth 
(ft  below 
land  surface) 

SM-Dd-49 

SM-73-3081 

Redgate 

38°16 

16" 

115 

10/20/78 

617 

76°36 

47" 

SM-Fe-31 

SM-73-3088 

Piney  Point 

38°08 

34" 

8 

10/18/78 

598 

76°30 

35" 

SM-Dd-50 

SM-73-3082 

Leonard town 

38°18 

07" 

90 

10/26/78 

578 

76°38 

00" 

SM-Df-71 

SM-73-3431 

Lexington  Park 

38°15 

27" 

65 

7/12/79 

560 

76°28 

31" 

SM-Bb-15 

SM-73-3430 

Charlotte  Hall 

38°28 

38" 

170 

7/18/79 

460 

76°47 

01" 

SM-Bb-22 

SM-73-3787 

Charlotte  Hall 

38°28 

38" 

170 

6/27/80 

240 

76°47 

01" 

SM-Dd-62 

SM-73-3786 

Redgate 

38°16 

16" 

115 

6/30/80 

382 

76°36 

47" 

SM-Dd-63 

SM-73-3785 

Redgate 

38°16 

15" 

116 

7/07/80 

377- 

76°36 

47" 

CA-Cc-57 

CA-73-2893 

Parren 

38°36 

05" 

139 

10/04/78 

577 

76°34 

46" 

CA-Fd-54 

CA-73-2892 

Calvert  Cliffs 

38°21 

19" 

129 

10/11/78 

698 

Park 

76°25 

60" 

CA-Db-47 

CA-73-3304 

Prince  Frederick 

38°32 

39" 

142 

7/26/79 

570 

76°35 

42" 

CA-Bb-27 

CA-73-3303 

Dunkirk 

38°43 

31" 

130 

8/03/79 

320 

76°39 

52" 

CA-Bb-28 

CA-73-3721 

Dunkirk 

38°43 

31" 

130 

6/25/80 

199 

76°39 

52" 

AA-Fd-43 

AA-74-1004 

Tracys  Landing 

38°46 

46" 

150 

8/08/79 

280 

76°35 

24" 

AA-Ed-45 

AA- 74-1005 

Davidsonville 

38°54 

06" 

100 

8/15/79 

157 

76°38 

39" 

\J  Geophysical  logs:  G  =  gamma  ray;  M  =  multipoint  electric; 

S  =  single  point  electric. 

Well  locations  are  given  on  plates  12  and  13. 

dicates  a  vertical  hydraulic  conductivity  of  l.OxlO-8 

Water  Levels 

ft/s  and  a  specific  storage  of  7.4xl0_5ft_1. 

The  upper  confining  bed  of  the  Aquia  aquifer  is 

Historical  water  levels  in  the  Aquia  aquifer 

the  Marlboro  Clay  and  the  silty-clay  portion  of  the 

recorded  by  Darton  (1896),  Clark  and  others  (1918), 

Nanjemoy  Formation.  The  thickness  of  this  confin- 

Overbeck  (1948), 

Bennion  and  Brookhart  (1949), 

ing  bed  is 

shown  in  figure  6.  Several  undisturbed 

Cook  and  others  (1952),  and  Martin  and  Ferguson 

cores  of  this  material  were 

taken  as  a  part  of  this 

(1953)  make  possible  an  approximate  reconstruction 

study.  The  results  of  laboratory  tests  to  determine 

of  the  prepumping  potentiometric  surface  (pi.  3). 

the  vertical  hydraulic  conductivity  of  this  material 

This  reconstruction  is  approximate  because  of  the 

are  shown  in  table  2.  The  values  range  from  10-9  to 

lack  of  measured  levels  in  much  of  the  area.  The 

10-10  ft/s.  Specific  storage  values  are  not  available 

precision  of  the 

measurements  available  is  also 

from  these  laboratory  tests. 

uncertain;  however,  because  many  of  these  measure- 
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Construction  data 

Static  water 
level  (ft 
below  land  sur¬ 
face  measured 
on  May  1 ,  1981 ) 

Casing 
diameter 
( in. ) 

upper-lower 

Screen 
posit  ion 
(ft  below 
land  surface) 

Core 
sample 
(ft  below 
land  surface) 

Geophysical 
logs  1/ 

Driller 

Aquifer 

U.S.G.S. 

Well 

No. 

4-2 

539-617 

236-236.5 

458-458.5 

S,  G 

Shannahan  Arte¬ 
sian  Well  Co. 

143.73 

Aquia 

SM-Dd-49 

4-2 

448-458 

- 

S,  G, 

M 

Shannahan  Arte¬ 
sian  Well  Co. 

37.48 

Aquia 

SM-Fe-31 

4-  2 

503-513 

198-198.5 

S,  G 

Shannahan  Arte¬ 
sian  Well  Co. 

124.36 

Aquia 

SM-Dd-50 

4-2 

550-560 

- 

S,  G, 

M 

Calvert  Well 
Drilling  Co. 

126.66 

Aquia 

SM-Df-71 

4-2 

450-460 

- 

S,  G, 

M 

Calvert  Well 
Drilling  Co. 

165.79 

Aquia 

SM-Bb-1 5 

4-2 

208-218 

S,  G, 

M 

Calvert  Well 
Drilling  Co. 

15.78 

Piney  Point- 
Nan  jemoy 

SM-Bb-22 

4-2 

348-358 

- 

S,  G 

Calvert  Well 
Drilling  Co. 

113.60 

Piney  Point- 
Nan  jemoy 

SM-Dd-62 

4-2 

346-356 

S,  G 

Calvert  Well 
Drilling  Co. 

115.05 

Piney  Point- 
Nan  jemoy 

SM-Dd-63 

4-2 

511-521 

~ 

S,  G, 

M 

Shannahan  Arte¬ 
sian  Well  Co. 

145.23 

Aquia 

CA-Cc-57 

4-2 

638-648 

215-215.5 

516-516.5 

S,  G, 

M 

Shannahan  Arte¬ 
sian  Well  Co. 

156.31 

Aquia 

CA-Fd-54 

4-2 

560-570 

- 

S,  G, 

M 

Calvert  Well 
Drilling  Co. 

152.90 

Aquia 

CA-Db-47 

4-2 

310-320 

- 

S,  G, 

M 

Calvert  Well 
Drilling  Co. 

138.59 

Aquia 

CA-Bb-27 

4-2 

160-170 

S,  G, 

M 

Calvert  Well 
Drilling  Co. 

82.56 

Piney  Point- 
Nan  jemoy 

CA-Bb-28 

4-2 

270-280 

- 

S,  G, 

M 

Calvert  Well 
Drilling  Co. 

146.62 

Aquia 

AA-Fd-43 

4 

147-157 

- 

S,  G, 

M 

Calvert  Well 
Drilling  Co. 

66.21 

Aquia 

AA-Ed-45 

ments  are  consistent  with  each  other,  it  seems  pro¬ 
bable  that  the  measurements  are  accurate  within 
about  5  ft.  The  most  notable  features  of  plate  3  are 
the  potentiometric  highs,  which  correspond  to 
topographically  high  areas,  and  the  potentiometric 
lows,  which  correspond  to  low  areas  near  the 
Potomac  and  Patuxent  Rivers  and  the  Chesapeake 
Bay. 

Otton  (1955)  prepared  the  first  potentiometric 
map  of  the  Aquia  in  Southern  Maryland.  This  map 
(pi.  4)  shows  cones  of  depression  near  the  towns  of 
Lexington  Park  and  Leonardtown  in  St.  Marys 
County  in  approximately  1952. 

The  cone  of  depression  in  the  Lexington  Park  area 
remained  nearly  constant  from  1952  to  1967,  when 
the  potentiometric  surface  was  again  measured  by 
Weigle  (1970).  However,  water-levels  declined  from 


10  to  20  ft  near  Prince  Frederick  in  Calvert  County. 
In  St.  Marys  County,  the  water  level  declined  about 
20  ft  near  Leonardtown  because  of  increased  munic¬ 
ipal  pumpage. 

From  1967  to  1980,  there  were  general  water-level 
declines  in  the  Aquia  throughout  Southern  Mary¬ 
land.  The  potentiometric  surface  of  the  Aquia  was 
measured  May  19-23,  1980  (pi.  5).  The  most  signifi¬ 
cant  water-level  declines  in  this  period  were  in  the 
Lexington  Park  area  and  near  Cove  Point  and 
Prince  Frederick  in  Calvert  County.  The  general 
declines  in  water  levels  were  probably  caused  by 
growth  of  large  cones  of  depression  and  increasing 
domestic  and  industrial  pumpage. 

Eleven  observation  wells  were  drilled  during  this 
project  to  monitor  water  levels  in  the  Aquia  aquifer. 
Records  of  these  wells  are  listed  in  table  3  and  their 
locations  are  shown  on  plate  12. 
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Pumpage 

It  is  not  known  exactly  when  the  first  water  wells 
were  drilled  into  the  Aquia.  Darton  (1896)  lists  wells 
in  Nanjemoy  and  Oakley  drilled  into  the  Aquia 
aquifer.  As  these  were  flowing  wells,  however,  it  is 
unlikely  that  they  were  equipped  with  pumps.  Clark 
and  others  (1918)  list  many  more  wells  drilled  into 
the  Aquia  aquifer.  A  few  of  these  were  equipped 
with  suction  pumps,  but  most  relied  on  natural  arte¬ 
sian  pressure  to  produce  water.  Because  most  of 
these  wells  flowed  at  only  2  to  5  gal/min,  total  yield 
from  the  Aquia  at  this  time  was  probably  less  than 
200,000  gal/d  in  Southern  Maryland. 

Large  scale  pumping  of  the  Aquia  aquifer  began 
during  World  War  II  when  the  U.S.  Navy  built 
several  training  facilities  in  St.  Marys  County  near 
Lexington  Park  and  in  southern  Calvert  County. 
From  1941  to  1943,  pumpage  from  the  Aquia  in¬ 
creased  to  0.54  Mgal/d.  By  1945,  pumpage  increased 
to  1.8  Mgal/d.  After  World  War  II,  Aquia  pumpage 
from  the  Naval  facilities  declined  to  0.87  Mgal/d  by 
1951.  From  1951  to  1966,  pumpage  from  the  Aquia 
increased  gradually.  Much  of  this  increase  came 


from  domestic  pumpage.  From  1966  to  1980,  munic¬ 
ipal  and  domestic  use  of  Aquia  water  increased 
rapidly. 

Table  4  summarizes  the  pumpage  history  of  the 
Aquia  aquifer  from  1941  to  1980.  This  table  lists 
documented  Aquia  pumpage  from  well  fields  that 
use  more  than  5,000  gal/d.  The  locations  of  these 
well  fields  are  shown  in  figure  7.  Also  shown  are  the 
county  planning  districts  of  St.  Marys,  Calvert,  and 
Anne  Arundel  Counties  where  Aquia  pumpage  has 
been  documented.  Included  with  table  4  are  esti¬ 
mates  of  domestic  and  small  industrial  pumpage 
organized  by  county  planning  district.  The  pro¬ 
cedure  used  to  estimate  this  pumpage  was  similar  to 
the  method  employed  by  Williams  (1976).  First,  the 
approximate  number  of  domestic  wells  for  each 
county  was  determined  from  the  State  of  Maryland 
Water  Resources  Administration  fact  sheets 
WS001-WS004.  The  total  number  of  wells  in  use 
was  then  divided  among  each  county  planning  dis¬ 
trict.  Each  domestic  well  was  considered  to  pump  an 
average  of  250  gal/d.  Finally,  the  approximate 
domestic  pumpage  was  calculated  by  multiplying 
the  number  of  domestic  wells  in  each  planning 
district  by  the  average  daily  usage. 
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Table  4.— Major  users  of  the  Aquia  aquifer  and  their  pumping  rates 
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planning  district  from  the  Anne  Arundel  County  master  plan  for  water  supply  and  sewerage  systems.  July  1976.  £late 


Table  4.*— continued 


X 

o 

•H 

0) 

cl 


cd 

O 

A) 

r-» 

X  to 

A)  1 

1 

r-H 

(X  o 

vO 

B  *— < 

ON 

p 

CL  U 

r-H 

P 

60 


p  x) 

3  a) 
o  -u 
p 


4 


cl 

O 

Li 

cl 

cl 

«d 


o 

00 

I 


60 

o 

cd 

vO 

u 

1 

A) 

r—l 

> 

in 

cd 

ON 

r-H 

X 

r—l 

cd 

60 

o 

m 


<r 

ON 


s 

i 

•rH 

X 

cd 

X 


P 

60 

cd 

u 

P 

<i 


CO 

•  H  • 
0^0 

•  A)  Z 
CO  3: 


U  O  Ph 

P  *H  <J 
4-»  4-J  O 
cd  cd  w 


M  4J  • 

CL -H  o 
O  S  2 

»h  M 

a  a) 
o  cl  cl. 

< 


i  m\ 

p 

cd  cd 

T— I  d) 

cl  u 
cd 

■U  60 

P  P 

0  *H 

O 

U 


p 

0  r- 

o 

o 

•H 

u  p 

4-1 

Uh  u 

Cd 

p 

o 

•  60 

o 

O  -H 

►J 

Z  UH 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

ON  /— \ 

ON 

00 

co 

o 

/■s 

vO 

r-H 

m 

m 

/rs 

/T~S 

«r— \ 

CM 

ON 

Vs  H 

«— t 

rs 

Vs 

*— t 

vs  rH 

Vs 

Vs 

Vs 

rH 

rH 

rH 

Vs 

Vs 

O  W 

in 

00  W 

o 

oo 

o  ^ 

m 

CO 

V-/ 

V— ' 

vD 

m 

00 

r^. 

r-H 

CM 

vO 

CM 

co 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o  ^ 

ON 

in 

vO 

00 

X— \ 

r-H 

<T 

o 

X—s 

co 

H 

»-H 

•s  w-i 

Vs 

•— t 

Vs  rH 

Vs 

o 

Vs 

rH 

rH 

rH 

Vs 

o 

.O  w 

m 

v_/ 

in  w 

vO 

oo 

N-/ 

n- 

r-H 

o 

W 

W 

vO 

o 

ON 

CM 

CM 

CO 

*-H 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

O 

00 

✓“N 

/^N 

/T-S 

o 

/T-n 

/T~N 

Vs  H 

rv 

<-• 

•s  r-i 

O 

o 

rH 

o  ^ 

•s 

o 

rH 

rH 

rH 

O 

o 

o  ^ 

CM 

v-^ 

V-/ 

o 

V-/ 

V— ✓ 

v> 

o 

vO 

CM 

r-H 

m 

O 

o 

o 

O 

o 

o 

O 

/*"N 

/~s 

<r\ 

o 

/~s 

✓— N 

o  ^ 

f* 

o  - 

O 

o 

O  ^ 

o 

o 

Vs 

rH 

rH 

rH 

O 

o 

*t  H 

o 

V/ 

>—✓ 

V-/ 

m 

V_/ 

W 

O  W 

vO 

m 

r- 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

O 

o  ^ 

o 

o 

/^s 

O  /~N 

O 

o 

o 

/-s 

o 

o 

»t  H 

1 

•» 

vs 

Vs 

rH 

Vs  rH 

Vs 

Vs 

Vs 

rH 

rH 

rH 

Vs 

Vs 

O  W 

vO 

o 

CM 

O  W 

o 

m 

m 

v-/ 

V— ' 

00 

in 

m 

vO 

vO 

m 

o 

r— H 

CM 

CM 

•-H 

r-H 

CM 

r-H 

CM 

o 

O 

O 

o 

O 

o 

o 

o 

o 

o 

o 

O 

O 

o 

O 

o 

o 

o 

o 

o 

o  ^ 

O 

o 

m 

/■"S 

O 

o 

o 

o 

<rv 

/■> 

in 

o 

Vs  rH 

1 

Vs 

Vs 

rH 

Vs  rH 

Vs 

vs 

Vs 

rH 

rH 

rH 

Vs 

Vs 

in  w 

<T  ^ 

CM 

r-H 

\w/ 

m 

o 

i— H 

o 

>-✓ 

S-H' 

V-/ 

m 

m 

vO 

<t 

ON 

CM 

00 

r-H 

in 

t-H 

r-H 

r-H 

VO  *-H 

o 

*-H 

r-H  VO 

ON 

I"- 

on  00 

CM 

in 

CO 

m 

vO 

CM 

co 

1  1 

1 

CM 

■  i 

CM 

■ 

CO 

1 

co 

1 

CO  r-H 

I  1 

r-H 

| 

1 

1 

■ 

co 

CM 

■ 

CO 

X  X 

1 

o 

u 

1  1 
rO  rO 

i 

o 

1 

o 

1 

o 

1  1 

X  A) 

X 

-o 

o 

1 

o 

■ 

o 

O 

i 

-O 

O 

o 

u 

Q  Q 

o 

Q 

Q 

pH  Ph 

w 

Q 

PQ 

PQ 

PQ 

CQ 

Q 

CQ 

cd 

o 

CO 

r-H 

CM 

r-H 

CM 

CM 

CM 

r-H 

in 

CO 

00 

r-H 

Ph 

1 

(X 

1 

Pi 

1 

P-) 

1 

Pi 

1 

CL 

1 

CL 

1 

04 

1 

Oj 

1 

PH 

*3 

< 

c 

< 

3 

3 

< 

3 

3 

3 

o 

i 

o 

o 

o 

o 

o 

o 

O 

i 

o 

■ 

o 

1 

o 

1 

CM 

CM 

CM 

co 

ON 

o 

CM 

i 

o 

i 

CM 

1 

co 

1 

co 

CO 

i 

co 

I 

in 

in 

m 

i 

vO 

vO 

n- 

i 

• 

■ 

i 

cd 

1 

cd 

cd 

cd 

i 

cd 

cd 

l 

cd 

i 

cd 

■ 

cd 

cd 

■ 

P 

o 

u 

o 

o 

o 

o 

o 

o 

o 

o 

u 

< 

o 

pq 

C 

« 

< 

< 

CQ 

CJ 

PQ 

o 

■u 

CO 

CO 

14H 

u 

A) 

MH 

r-H 

A ) 

C 

H 

•iH 

cd 

4J 

to 

O 

r-H 

■u 

cd 

C 

•rH 

co 

u 

•H 

z 

cd 

4-J 

C 

CL 

CL 

Cd 

r-H 

o 

r-H 

CO 

CO 

B 

o 

o 

B 

A) 

O 

4H 

o 

o 

o 

o 

X 

4H 

o 

rC 

> 

rP 

rP 

r-H 

C 

•H 

o 

o 

u 

o 

cd 

r-H 

to 

p 

to 

cd 

CO 

CO 

P4 

4-» 

u 

4-J 

cd 

4-> 

A) 

c 

•H 

P4 

c 

PQ 

rP 

I  >> 

1 

p 

r-H 

p 

60 

-u 

o 

CO 

•H 

A ) 

vP 

o 

A) 

•rH 

CO 

tO*H 

o 

4-J 

c 

4-J 

rX 

o 

o 

vX 

PS 

>  »“< 

> 

CO 

cd 

Z 

cd 

cd 

»H 

cd 

p 

Cd  *H 

cd 

4J 

•rH 

CL 

A) 

A) 

4—1 

A) 

A) 

4-J 

O 

Z  u 

Z 

M 

4-J 

B 

t-H 

CL 

CQ 

4-J 

CL 

M 

cd 

A) 

C 

o 

A) 

cd 

A) 

P 

cd 

A) 

w 

•  pH 

• 

> 

A) 

o 

CO 

CO 

60 

> 

<P 

CO 

> 

•H 

CO 

co 

r-H 

•rH 

CO 

A) 

C 

r-H 

A) 

r-H 

u 

• 

• 

cd 

u 

cd 

vC 

O 

cd 

rC 

P 

p 

P3 

u 

CO 

o 

O 

hQ 

CJ 

o 

o 

04 

m 

vO 

r>. 

00 

ON 

o 

r-H 

CM 

CO 

m 

r-H 

r-H 

r-H 

r-H 

r-H 

r-H 

O 

O 

m  ^  ^ 


00  W  w 


o 

o 

O  '"N  ''"'N 
r".  v-x  'w' 

in 


o 

o 

m  /— n 

CM  w  w 
CO 


ON  o  H 
H  CM  CM 

I  I  I 

P  A)  P 


r-H 

I 

5! 

O 


CO 

I 

cd 

o 


< 


p 

o 


cd 

u 

o 

cl 

u 

o 

u 


•0 

p 

r-H 

o 

o 


o 

o 

O  '“N 

ON  w 

m 


o  o  o 

o  o  o 

o  o  o 

m  o  o 

oo  o  o 

«-H  *~H  in 


o 

o 

o  ^ 


o 

o 

o  ^ 


ON  W 


m  w 


o  o  o 

o  o  o 

o  o  o 

o  o  o 

m  oo  o 

co 


o  o 
<f  m 


o  o 
o  o 
o  o 

o  o 
co  m 


o 

o 

o 

o 

o 


o 

o 

o 

o 

m 


o 

o 

o  ^ 

m  w 
vO 

oo 


o 

o 

O 


o 

o 

O 


VO  'W / 

t— H 

CM 


O  W 

o 

vO 


•— 1  <t 
I  I 

rd  -O 

Q  Q 


co 

CM  CM 
I  I 
X  X 
W  W 


I 


m 

I 

3 

O 

I 

■vf 

I 

cd 

o 


o 

r-H 

I 

% 

o 

ON 

vO 

I 

cd 

u 


CQ 


< 


<  pq  o 


V* 

o 

•rH 

u 

A) 

X 

P 

u 

04 

p 

o 

p 

•rH 

U 

CH 


o 

•  rH 

w 

u 

UH 

4-1 

mh 

o 

•rH  ^ 

p 

r-H  4-J 

r-H 

O  P 

w 

P 

4J  r-H 

X 

U  Oh 

p 

A) 

p 

>  u 

^H  P 

CO 

P  3 

P 

o  o 

e> 

^  04 

p 

to  u 

Jh 

p  p 

o 

p  p 

B 

CL  *-H 

•H 

B  u 

4-J 

o  0 

r-H 

o  z 

P 

CQ 

CO 

P 

X 

p 


cd 

B 

CO 

X 

§ 


co 

A) 

B 

o 

Q 


m 


20 


_i/  Pumpage  included  in  the  above  value. 

3/  Population  planning  area  from  the  Calvert  County  comprehensive  water  and  sewerage  plan.  1978  update,  Figure 
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jJ  Pumpage  included  in  the  above  value. 

4/  Sanitary  districts  from  the  St.  Marys  County  comprehensive  water  and  sewerage  plan;  1977  update,  Map 


Piney  Point-Nanjemoy  Aquifer 
Nomenclature 

The  Piney  Point  Formation  in  Southern  Maryland 
is  lithologically  and  paleontologically  distinct  from 
the  underlying  Nanjemoy  Formation  (Otton,  1955). 
However,  as  has  been  pointed  out  by  Weigle  (1970) 
and  Williams  (1979),  the  upper  sands  of  the  Nan¬ 
jemoy  Formation  are  directly  overlain  by  Piney 
Point  Formation  sands  in  many  places.  This  sug¬ 
gests  that  the  two  formations  are  hydraulically  con¬ 
nected  and  function  as  a  single  aquifer.  In  view  of 
this,  Weigle  (1970)  proposed  the  name  “Piney  Point- 
Nanjemoy  hydraulic  unit”  for  use  in  Southern 
Maryland.  This  report  will  follow  that  convention. 
Where  it  is  used,  the  term  “Piney  Point-Nanjemoy 
aquifer”  is  synonymous  with  the  term  “Piney  Point- 
Nanjemoy  hydraulic  unit.” 

Regional  Extent 

The  Piney  Point-Nanjemoy  aquifer  is  used  as  a 
source  of  ground  water  in  Calvert  and  St.  Marys 
Counties.  On  the  Eastern  Shore  of  Maryland,  the 
Nanjemoy  Formation  is  predominantly  clay  and  the 
term  “Piney  Point-Nanjemoy  aquifer”  does  not  apply. 

Figure  8  shows  the  altitude  of  the  top  of  the  Piney 
Point-Nanjemoy  aquifer.  Its  thickness  is  shown  in 
figure  9.  In  northern  Calvert  and  northern  St. 
Marys  Counties,  the  aquifer  is  relatively  thin  and 
unproductive.  The  aquifer  thickens  rapidly  to  the 
southeast  and  becomes  correspondingly  more  pro¬ 
ductive.  Southeast  of  Leonardtown  in  St.  Marys 
County,  wells  screened  in  the  Piney  Point- 
Nanjemoy  commonly  yield  more  than  60  gal/min. 

Transmissivity  and  Storage  Coefficient 

Because  the  Piney  Point-Nanjemoy  aquifer  is 
utilized  mainly  by  homeowners  (Hansen,  1972), 
aquifer  tests  are  not  commonly  performed  on  newly 
drilled  wells.  For  small-capacity  wells,  drillers  are 
usually  concerned  only  to  produce  a  short-term  yield 
of  about  20  gal/min  and  pumping  tests  to  determine 
the  long-term  yield  potential  are  deemed  unneces¬ 
sary.  Faced  with  these  difficulties,  Williams  (1979) 
constructed  a  transmissivity  map  of  the  Piney 
Point-Nanjemoy  aquifer  in  Southern  Maryland  based 
upon  the  few  available  single-well  aquifer  tests  and 
specific  capacities  of  wells  as  reported  by  drillers. 
The  portion  of  the  transmissivity  map  which  covers 
the  project  area  of  this  report  is  shown  in  figure  10. 

Only  a  few  storage  coefficient  values  have  been 
determined  for  the  Piney  Point-Nanjemoy  aquifer. 


These  have  been  tabulated  by  Hansen  (1972)  and 
range  from  0.0003  to  0.0004.  A  multiple-well  pump¬ 
ing  test  of  the  Piney  Point-Nanjemoy  aquifer  near 
Leonardtown  in  St.  Marys  County  indicated  a 
storage  coefficient  of  0.0004. 

Upper  and  Lower  Confining  Beds 

The  Piney  Point-Nanjemoy  aquifer  is  overlain 
throughout  Southern  Maryland  by  silt  and  clay  of 
the  Chesapeake  Group.  These  sediments  act  as  an 
upper  confining  bed  to  the  Piney  Point-Nanjemoy 
aquifer.  The  thickness  of  this  confing  bed  is  shown 
in  figure  11.  In  general,  the  Chesapeake  Group 
sediments  are  thickest  in  the  upland  areas  of  St. 
Marys,  Charles,  and  Calvert  Counties  and  are  thin¬ 
nest  in  the  valley  of  the  Patuxent  River.  This  sug¬ 
gests  some  post-deposition  erosional  influence  on 
their  thickness. 

The  vertical  hydraulic  conductivity  of  the  Ches¬ 
apeake  Group  has  been  determined  in  a  few  loca¬ 
tions  by  laboratory  methods.  These  values  are 
shown  in  table  2. 

The  specific  storage  of  the  Chesapeake  Group 
sediments  has  been  determined  in  only  one  location 
in  Southern  Maryland.  This  is  at  the  Old  Leeland 
Road  test  well  in  Prince  Georges  County  (Hansen, 
1977)  where  specific  storage  is  10~5  ft.-1. 

The  lower  confining  beds  of  the  Piney  Point- 
Nanjemoy  aquifer  are  the  silty  clay  beds  of  the  Nan¬ 
jemoy  Formation  and  the  Marlboro  Clay.  The  hy¬ 
draulic  properties  of  these  beds  are  discussed  in  the 
section  on  the  upper  confining  bed  of  the  Aquia 
aquifer. 

Water  Levels 

Water  levels  in  the  Piney  Point-Nanjemoy  aquifer 
show  a  similar  history  to  the  water  levels  of  the 
Aquia  aquifer.  The  prepumping  potentiometric  map 
(pi.  6),  as  determined  from  historical  water  levels 
(Darton,  1896;  Clark  and  others,  1918;  Overbeck, 
1948;  and  Otton,  1952),  shows  potentiometric  highs 
in  upland  areas  and  potentiometric  lows  conform 
with  the  topographic  lows  in  river  valleys. 

The  Piney  Point-Nanjemoy  was  developed  in  the 
early  1900’s.  Its  relative  shallowness,  combined 
with  the  fact  that  many  of  the  early  wells  flowed, 
made  the  Piney  Point-Nanjemoy  aquifer  attractive 
to  many  users.  Use  of  this  aquifer  increased  rapidly 
during  World  War  II.  Water  levels  in  the  Lexington 
Park  and  Solomons  area  declined,  and  many  arte¬ 
sian  wells  ceased  flowing.  During  the  late  1940’s, 
the  Navy  reduced  its  pumpage  of  the  Piney  Point- 
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Figure  7.— Locations  of  well  fields  pumping  more  than  5,000  gal/d  from  the  Aquia 
aquifer  and  boundaries  of  county  planning  districts. 
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Figure  8. — Altitude  of  the  top  of  the  Piney  Point-Nanjemoy  aquifer. 
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Figure  9.— Thickness  of  the  Piney  Point-Nanjemoy  aquifer. 
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Figure  10.—Transmissivity  of  the  Piney  Point-Nanjemoy  aquifer. 
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76”  30' 


76”  15' 


76”45' 


Figure  12.— Location  of  well  fields  pumping  more  than  5,000  gal/d  from 
the  Piney  Point-Nanjemoy  aquifer  and  boundaries  of 
county  planning  districts. 
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Nanjemoy  aquifer  in  response  to  citizens’  com¬ 
plaints  about  declining  water  levels.  Some  pumping 
was  continued,  however,  and  by  1952  a  cone  of 
depression  about  20  ft  below  sea  level  had  developed 
around  Patuxent  Naval  Air  Training  Station 
(PNATS).  The  1952  potentiometric  surface  is  shown 
on  plate  7. 

From  1952  to  1980,  pumping  from  the  Piney 
Point-Nanjemoy  aquifer  increased.  Much  of  this 
pumpage  came  from  homes  and  other  small-capacity 


wells.  The  1980  potentiometric  surface  of  the  Piney 
Point-Nanjemoy  aquifer  is  shown  on  plate  8.  The 
pumpage  history  of  the  Piney  Point-Nanjemoy 
aquifer  from  1941  to  1980  is  summarized  in  table  5. 
The  locations  of  well  fields  producing  more  than 
5,000  gal/d  are  shown  in  figure  12. 

Four  Piney  Point-Nanjemoy  observation  wells 
were  drilled  during  this  project  to  observe  water 
levels.  Records  of  these  wells  are  listed  in  table  3 
and  the  well  locations  are  shown  on  plate  13. 
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Table  5.— Major  users  of  the  Piney  Point-Nanjemoy  aquifer  and  their  pumping  rates 
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Table  5.— continued 
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Hydrologic  Effects  of  Pleistocene  Erosion 
in  Southern  Maryland 

Pleistocene  Erosion 

The  worldwide  lowering  of  sea  level  during  Pleis¬ 
tocene  time  and  the  consequent  development  of  a 
river  system  approaching  equilibrium  with  that  sea 
level  has  had  a  strong  influence  on  the  present 
topography  and  drainage  system  of  Southern  Mary¬ 
land.  Hack  (1955,  1957)  has  shown  that  the  Chesa¬ 
peake  Bay  Estuary  is  a  drowned  Pleistocene  river 
system  and  that  this  river  system  approached  equi¬ 
librium  with  a  sea  level  300  to  400  ft  below  present 
sea  level.  In  the  process  of  approaching  this  dynam¬ 
ic  equilibrium  with  the  lower  sea  level,  there  was 
considerable  downward  erosion  of  river  channels.  In 
many  areas  of  Southern  Maryland,  these  river  chan¬ 
nels  have  eroded  through  sediments  of  Miocene, 
Eocene,  and  Paleocene  age. 

The  effects  of  this  erosion  can  clearly  be  observed 
at  some  locations.  One  such  location  is  at  the  Ches¬ 
apeake  Bay  Bridge.  Figure  13  shows  a  geologic 
cross  section  at  this  site  that  was  constructed  from 
bore  hole  data  gathered  during  bridge  construction. 
This  figure  shows  a  buried  channel  that  has  a  base 
about  200  ft  below  sea  level.  Hack  (1957)  has  inter¬ 
preted  this  channel  as  representing  the  Pleistocene 
Susquehanna  River.  The  channel  of  this  ancient 
river  has  completely  cut  through  the  Aquia  aquifer 
at  this  location. 

Other  examples  of  the  effects  of  Pleistocene  ero¬ 
sion  can  be  observed  near  Benedict,  Md.  (fig.  14), 
and  near  Solomons,  Md.  (fig.  15).  Bore  hole  data 
gathered  at  these  sites  show  that  Chesapeake  Group 
sediments  have  been  eroded  and  buried  by  sedi¬ 
ments  of  Pleistocene  age. 

Areal  Distribution  of  Pleistocene  Erosional  Channels 

To  evaluate  the  effects  that  Pleistocene  erosional 
channels  have  on  the  hydrogeology  of  the  Aquia  and 
Piney  Point-Nanjemoy  aquifers,  it  is  necessary  to 
determine  where  these  aquifers  or  their  confining 
beds  have  been  disturbed.  Hack  (1957)  estimates  the 
longitudinal  profiles  of  the  Pleistocene  Susquehan¬ 
na  River  and  some  of  its  major  tributaries.  Figue  16 
shows  the  approximate  locations  and  approximate 


depths  of  these  Pleistocene  channels,  which  are  as¬ 
sumed  to  coincide  with  the  present  day  deep  channel 
of  the  Chesapeake  Bay.  The  information  of  figure  16 
combined  with  a  knowledge  of  the  areal  distribu¬ 
tion,  altitude  of  top,  and  thickness  of  geologic  for¬ 
mations  makes  possible  an  estimate  of  where  these  for¬ 
mations  have  been  eroded  by  Pleistocene  channels. 

Schubel  and  Zabawa  (1972)  have  identified  at 
least  one  other  channel  in  the  Chesapeake  Bay 
region  using  continuous  seismic  reflection  profiling. 
It  is  possible,  therefore,  that  the  truncation  of  Ter¬ 
tiary  sediments  is  more  extensive  in  the  Bay  region 
than  can  be  inferred  from  figure  16,  which  considers 
only  those  channels  identified  by  Hack  (1957). 

Hydrogeologic  Effects  of  Pleistocene 
Erosional  Channels 

Disruption  of  aquifers  and  confining  beds  by 
Pleistocene  erosional  channels  appears  to  have  a 
significant  effect  on  head  distributions  in  the  Aquia 
and  Piney  Point-Nanjemoy  aquifers.  Where  aquifers 
have  been  breached  by  erosional  channels,  aquifer 
water  levels  will  tend  to  approach  the  river  water 
level  by  leakage  into  or  out  of  the  aquifer.  Similarly, 
where  confining  beds  have  been  breached,  increased 
leakage  into  or  out  of  the  aquifer  will  cause  a  head 
distribution  adjustment.  This  reasoning  is  consis¬ 
tent  with  the  observed  prepumping  potentiometric 
surfaces  of  the  Aquia  and  Piney  Point-Nanjemoy 
aquifers.  A  comparison  of  figure  16  with  plates  3 
and  6  shows  that  potentiometric  lows  correspond  to 
areas  in  river  channels  where  erosion  has  disturbed 
either  the  aquifers  or  their  confining  beds.  This  pro¬ 
cess  appears  to  be  one  of  the  major  controls  on  the 
natural  head  distribution  of  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers  in  Southern  Maryland. 

In  some  areas  close  to  the  Chesapeake  Bay, 
Pleistocene  erosion  has  increased  the  potential  for 
salt-water  intrusion  into  aquifers.  For  example, 
figure  13  shows  that  at  the  Chesapeake  Bay  Bridge, 
Pleistocene  river  sands  have  been  deposited  on  an 
erosional  surface  which  truncates  the  Aquia  aquifer. 
If  the  potentiometric  surface  of  the  Aquia  aquifer 
drops  significantly  below  sea  level  near  this  area,  in¬ 
trusion  of  brackish  water  into  the  Aquia  is  possible. 
Future  water-supply  decisions  in  these  areas  should 
consider  this  potential  for  contamination. 
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Figure  13.— Geologic  cross  section  at  the  Chesapeake  Bay  Bridge. 
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Figure  14.— Geologic  cross  section  at  the  Patuxent  River  Bridge. 
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Figure  15.— Geologic  cross  section  at  the  Thomas  Johnson  Bridge. 
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DIGITAL  MODEL 


Theory  of  Modeling  Gound-Water  Flow 

The  flow  of  ground  water  in  a  heterogeneous  and 
anisotropic  medium  can  be  described  mathematically 
as 


J_  (T  Jh)  +  _!  (T  _3h) 

3x  xx  3x;  3y  yy  dyJ 


J-  (bK  ) 

3z  zz  dz 


S1-^  +  bW(x*y,z,t) 

O  L 


in  which 


T  = 


T  .  = 


yy 


K  = 


h 

S' 

b 

t 

x,  y,  z 
W(x,  y,  z,  t) 


principal  transmissivity  tensor  in 
the  x  direction  (L2T-1); 
principal  transmissivity  tensor  in 
the  y  direction  (L2T-1); 
hydraulic  conductivity  tensor  in 
the  z  direction  (LT-1); 
hydraulic  head  (L); 
storage  coefficient  (dimensionless); 
thickness  of  hydraulic  unit  (L); 
time  (T); 

space  coordinates  (L); 
volumetric  flux  per  unit  volume 
(T-1)  (Trescott,  1975). 

The  volumetric  flux  term  W(x,  y,  z,  t)  allows  the 
model  to  adjust  calculated  head  values  to  different 
rates  of  pumpage  and  leakage  across  confining  beds. 

Equation  1  can  be  solved  by  subdividing  the 
region  into  blocks  where  the  medium  properties  are 
assumed  to  be  uniform.  The  continuous  derivatives 
are  replaced  by  finite-difference  approximations  at  a 
point  in  the  center  of  each  block.  The  result  is  N 
algebraic  equations  in  N  unknowns,  where  N  is  the 
number  of  blocks  representing  the  porous  medium. 
The  finite-difference  approximation  of  equation  1 
can  be  written  as 


A  t  =  time  increment  (T); 

i  =  index  in  the  y  dimension; 
j  =  index  in  the  x  dimension; 
k  =  index  in  the  z  dimension; 
h  =  hydraulic  head  (L); 

S'  =  storage  coefficient  (dimensionless); 
b  =  thickness  of  hydraulic  unit  (L); 

Txx  =  principal  transmissivity  tensor  in 
the  x  direction  (L2T-1); 

Tyy  =  principal  transmissivity  tensor  in 
the  y  direction  (L2T-1); 

Kzz  =  hydraulic  conductivity  tensor  in 
the  z  direction  (LT-1). 

This  set  of  algebraic  equations  can  be  solved  with  a 
digital  computer. 

The  digital  model  program  used  in  this  study  was 
developed  by  Trescott  and  Larson  (1976).  In  this 
program,  the  Strongly  Implicit  Procedure  (SIP) 
numerical  technique  is  used  to  solve  equation  2.  The 
derivation  of  the  SIP  algorithm  is  given  in  Wein¬ 
stein  and  others  (1969).  The  program  incorporates  a 
transient  leakage  calculation,  developed  and  docu¬ 
mented  by  Posson  and  others  (1980),  to  simulate  the 
effects  of  confining-bed  storage.  For  this  study,  the 
quasi  three-dimensional  approach  was  used.  This  ap¬ 
proach  incorporates  the  effects  of  vertical  flow 
through  confining  beds  in  the  vertical  components 
of  hydraulic  conductivity  of  adjacent  aquifers  (Tres¬ 
cott,  1975).  This  is  justified  if  the  ratio  of  aquifer 
hydraulic  conductivity  to  the  confining  bed’s  hy¬ 
draulic  conductivity  is  greater  than  100:1.  In  the 
aquifer  system  of  Southern  Maryland,  it  is  esti¬ 
mated  that  this  ratio  is  greater  than  1,000:1. 

Description  of  the  Model 


(h-  i.,  -  h,  . 
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in  which 

A  Xj  =  space  increment  in  the  x  direc¬ 
tion  for  column  j  (L); 

A  y;  =  space  increment  in  the  y  direc¬ 
tion  for  row  i  (L); 

A  zk  =  space  increment  in  the  z  direc¬ 
tion  for  layer  k  (L); 


Conceptual  Model 

A  digital  ground-water  flow  model  is  a  math¬ 
ematical  representation  of  an  aquifer  system  and  is 
constructed  by  (1)  quantifying  the  physical  charac¬ 
teristics  of  an  aquifer  system,  and  (2)  translating 
these  characteristics  into  a  form  which  can  be 
manipulated  by  a  computer.  Most  natural  ground- 
water  flow  systems,  however,  are  far  too  complex  to 
be  precisely  represented  in  this  manner.  As  a  prac¬ 
tical  matter,  therefore,  digital  flow  models  are  con¬ 
structed  from  idealized  and  simplified  represen¬ 
tations  of  aquifer  systems.  These  simplified  ver¬ 
sions  of  the  natural  systems  are  called  conceptual 
models. 

Figure  17  shows  schematically  the  conceptual 
model  which  was  used  to  construct  the  digital  model 
described  in  this  report  as  interpreted  from  geologic 
and  hydrologic  information  previously  discussed. 
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channels 
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Figure  17.— Schematic  diagram  showing  the  conceptual  model  used  in  constructing  the  digital  model. 


The  project  area  is  conceived  to  be  covered  by  a 
mantle  of  Pleistocene  sediments  which  acts  as  a 
water-table  aquifer.  This  water-table  aquifer  is 
recharged  directly  from  precipitation.  Most  of  the 
water  in  this  aquifer  is  ultimately  discharged  into 
streams  or  by  evapotranspiration.  A  small  amount 
of  water,  however,  percolates  vertically  and  re¬ 
charges  the  Aquia  and  Piney  Point-Nanjemoy 
aquifers.  Where  the  Aquia  aquifer  is  directly  over- 
lain  by  Pleistocene  sediments  (for  example,  the 
outcrop-subcrop  area),  this  vertical  leakage  is 
greater  than  where  the  Aquia  is  overlain  by  con- 
fining-bed  material.  Because  of  the  subsurface  trun¬ 
cation  of  the  Piney  Point  Formation  and  the  upper 
sands  of  the  Nanjemoy  Formation,  the  Piney  Point- 
Nanjemoy  aquifer  does  not  have  an  outcrop  area. 
Discharge  from  the  Aquia  and  Piney  Point-Nan¬ 
jemoy  aquifers  occurs  as  upward  leakage  through 
confining  beds  in  downgradient  areas.  Additional 
discharge  occurs  in  major  river  valleys,  such  as  the 
Patuxent  River,  where  confining-bed  material  has 
been  eroded  by  Pleistocene  channels. 

This  conceptual  model  is  consistent  with  the  main 
features  of  the  natural  flow  system  as  shown  on 
plates  3  and  6.  It  is  also  consistent  with  the  distribu¬ 
tion  and  lithology  of  the  Tertiary  sediments  as 
shown  on  plates  1  and  2.  It  should  be  remembered, 
however,  that  this  conceptual  model  represents  a 
simplification  of  the  natural  system.  It  therefore  in¬ 
troduces  the  possibility  of  error  in  areas  where  the 
simplifications  may  not  be  good  approximations. 
Future  studies  with  more  complete  data  will  possi¬ 
bly  add  refinements  to  the  conceptual  model  used  in 
this  report. 

Grid  Design 

The  model  area  was  divided  into  a  rectangular 
grid  having  3  layers,  28  rows,  and  55  columns.  A 
variable  grid  spacing  was  used  so  that  the  smallest 
grid  blocks  coincide  with  areas  where  greater  ac¬ 
curacy  was  needed.  A  multiplication  factor  of  1.5 
was  used  to  expand  the  block  size  outside  the  area  of 
interest.  The  smallest  block  in  plan  view  was  1  mi^ 
and  the  largest  was  25  mi^.  By  convention,  the 
point  at  the  center  of  each  block  is  called  the  node. 
The  layer-row-column  system  (I,J,K)  was  used  to 
label  each  node.  For  example,  the  index  for  the  node 
(1,5,8)  refers  to  the  center  of  the  block  which  is 
located  in  layer  1,  row  5,  column  8.  For  the  model  in 
this  report,  layer  1  represents  the  Aquia  aquifer, 
layer  2  represents  the  Piney  Point-Nanjemoy 
aquifer,  and  layer  3  represents  the  Pleistocene  water 
table  aquifer.  Each  input  value  (transmissivity,  ver¬ 
tical  hydraulic  conductivity,  and  so  forth)  assigned 
to  a  node  is  considered  to  be  the  average  value  over 


the  entire  block.  Similarly,  each  output  value  (hy¬ 
draulic  head,  drawdown )  is  also  an  average  value  for 
that  block.  The  grid  design  used  for  this  study  is 
shown  in  figure  18. 

Model  Boundaries 

The  first  step  in  translating  the  conceptual  model 
of  the  Aquia  and  Piney  Point-Nanjemoy  aquifer 
system  into  a  digital  model  was  to  select  ap¬ 
propriate  boundaries.  These  boundaries  were  placed 
so  that  they  were  consistent  with  the  hydrogeologic 
framework  of  the  aquifer  system.  Because  natural 
aquifer  boundaries  in  Southern  Maryland  do  not 
coincide  with  the  study  area,  the  area  modeled  was 
considerably  larger  than  the  project  area. 

The  Aquia  aquifer  has  several  natural  boundaries 
which  were  utilized  in  the  modeling  process.  The 
eastern  boundary  of  this  aquifer  coincides  with  the 
change  from  facies  I  to  facies  III  (fig.  4)  and  was 
modeled  as  no-flow.  The  western  boundary  coincides 
with  the  Aquia  outcrop-subcrop  area  and  was 
modeled  as  a  no-flow  boundary.  The  boundary  to  the 
north,  where  the  Aquia  aquifer  extends  through 
Queen  Annes  County  to  Delaware,  was  assumed  to 
be  no-flow  and  was  placed  far  enough  away  from  the 
study  area  to  have  minimal  effect  on  model  results. 
Although  the  Aquia  Formation  extends  to  the  south 
into  Virginia,  its  lithology  changes  to  a  poorly 
sorted  fine  silty  sand  (Drobnyk,  1965)  that  is  not 
generally  utilized  as  an  aquifer.  No  data  on  the  yield 
characteristics  of  the  Aquia  are  available  for  this 
area  of  Virginia.  However,  based  on  the  lithology,  it 
is  assumed  that  the  transmissivity  is  relatively  low. 
It  was  judged  that  this  area  of  low  transmissivity 
could  be  best  simulated  by  treating  it  as  a  no-flow 
boundary.  The  location  of  simulated  no-flow  bound¬ 
aries  for  the  Aquia  aquifer  are  shown  on  figure  18. 

The  Piney  Point-Nanjemoy  aquifer  also  has  natu¬ 
ral  boundaries  which  were  utilized  in  the  modeling 
process.  The  western  boundary,  which  was  modeled 
as  no-flow,  approximates  the  subsurface  truncation 
of  the  Piney  Point-Nanjemoy  aquifer.  The  southern 
boundary  borders  an  area  of  low  transmissivity  and 
is  also  modeled  as  no-flow  (Williams,  1979).  The 
northeastern  and  eastern  boundaries  are  artificial 
since  the  Piney  Point  Formation  extends  in  this 
direction  up  to  New  Jersey.  However,  no-flow 
boundaries  in  this  direction  were  placed  far  enough 
away  so  that  they  would  have  little  or  no  effect  on 
model  results  in  the  project  area.  The  location  of 
simulated  no-flow  boundaries  for  the  Piney  Point- 
Nanjemoy  aquifer  are  shown  in  figure  18. 

The  lower  boundary,  which  separates  the  base  of 
the  Aquia  aquifer  from  the  underlying  Cretaceous 
sediments,  was  modeled  as  a  no-flow  boundary. 
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Figure  18.— Finite  difference  grid  and  simulated  no-flow  boundaries  of  the  digital  model. 
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Figure  19.— Hydrographs  of  observation  wells  in  the  Pleistocene  water-table  aquifer. 
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Some  vertical  leakage  might  exist  between  the 
Aquia  aquifer  and  the  underlying  Magothy  aquifer. 
However,  the  thick  clay  layers  of  the  Severn  Forma¬ 
tion  that  separates  the  Aquia  and  Magothy  aquifers 
are  believed  to  restrict  leakage  through  the  lower 
boundary. 

The  upper  boundary  of  the  aquifer  system  was 
modeled  as  a  constant-head  boundary.  Conceptual¬ 
ly,  the  upper  boundary  simulates  the  Pleistocene 
water-table  aquifer  in  Southern  Maryland.  The 
water-table  aquifer  recharges  the  underlying  con¬ 
fined  aquifer  system  by  vertical  leakage.  The  rate  of 
vertical  leakage  depends  on  the  head  gradient 
across,  the  thickness,  and  vertical  hydraulic  conduc¬ 
tivity  of  the  confining  beds.  An  assumption  implicit 
in  this  treatment  of  the  upper  boundary  is  that 
heads  in  the  Pleistocene  water-table  aquifer  remain 
relatively  constant.  Figure  19  shows  three  hydro¬ 
graphs  of  water-table  wells  in  Southern  Maryland. 
The  hydrographs  indicate  that  although  short-term 
variations  occur,  the  long-term  water-level  trend  is 
relatively  constant.  This  long-term  trend  indicates 
that  the  constant-head  assumption  for  the  upper 
boundary  of  the  model  is  reasonable. 

Figure  20  summarizes  the  boundary  conditions 
used  in  the  digital  model.  The  upper  Pleistocene 


aquifer  (layer  3)  is  modeled  as  a  constant-head 
boundary.  The  subsurface  truncation  of  the  Piney 
Point-Nanjemoy  aquifer  (layer  2)  is  modeled  as  a  no¬ 
flow  boundary.  Areas  of  layer  2  overlying  the  Aquia 
aquifer  (layer  1)  northwest  of  the  subsurface  trunca¬ 
tion  were  modeled  as  a  constant-head  boundary  in 
order  to  simulate  the  Pleistocene  water-table  aquifer 
in  areas  where  the  Piney  Point-Nanjemoy  aquifer  is 
absent.  In  areas  where  Pleistocene  sands  directly 
overlie  Aquia  sands,  a  thin  (1  ft)  confining  bed  was 
considered  to  separate  layers  1  and  2  of  the  model. 
The  subsurface  facies  change  of  the  Aquia  aquifer 
was  modeled  as  a  no-flow  boundary. 

Transmissivity  and  Storage  Coefficient 

Values  of  transmissivity  and  storage  coefficient 
for  the  Aquia  and  Piney  Point-Nanjemoy  aquifers 
are  discussed  in  the  Hydrogeology  section.  The 
transmissivity  distributions  of  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers  are  shown  in  figures  5  and 
10,  respectively. 

Average  transmissivity  values  for  the  Aquia  and 
Piney  Point-Nanjemoy  aquifers  were  estimated 
from  figures  5  and  10  for  each  block  of  the  model 
grid.  These  values  were  then  entered  into  the  model 


Subsurface  truncation  of  the 
Piney  Point  -Nanjemoy  aquifer 


Figure  20.— Schematic  diagram  showing  boundary  conditions  used  in  the  digital  model. 
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as  a  matrix.  The  value  of  lxl(M  is  a  representative 
value  of  storage  coefficient  for  the  confined  portions 
of  the  Aquia  aquifer  and  was  used  throughout  the 
modeled  area.  A  value  of  3x1  CM  is  a  representative 
storage  coefficient  value  for  the  Piney  Point- 
Nanjemoy  aquifer  and  was  used  throughout  the 
modeled  area. 

Starting  Heads 

The  prepumping  potentiometric  surfaces  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers  were  es¬ 
timated  on  the  basis  of  historical  water  levels  (pis.  3 
and  6).  Average  head  values  for  each  grid  block  were 
estimated  from  these  figures  and  used  as  starting 
head  values  for  model  simulations. 

The  Pleistocene  water-table  aquifer  provides 
recharge  to  the  Aquia-Piney  Point-Nanjemoy  sys¬ 
tem  in  Southern  Maryland.  The  potentiometric  sur¬ 
face  of  this  unit  was  estimated  from  the  altitude  of 
perennial  streams  that  intersect  the  water  table  (pi. 
9).  Average  head  values  for  each  grid  block  were 
entered  into  the  model  as  a  constant-head  matrix. 

Vertical  Hydraulic  Conductivities  of 

Confining  Beds 

The  upper  confining  beds  of  the  Aquia  aquifer  are 
the  clays  and  silts  of  the  Marlboro  Clay  and  Nan- 
jemoy  Formations.  The  Piney  Point-Nanjemoy  is 
confined  above  by  the  silts  and  clays  of  the  Ches¬ 
apeake  Group.  Some  values  of  vertical  hydraulic 
conductivity  for  these  confining  beds  are  available 
from  laboratory  core  analyses  (table  2).  For  the  most 
part,  however,  vertical  hydraulic  conductivities 
were  determined  by  calibration  of  the  model.  In  all 
cases  the  vertical  hydraulic  conductivities  as  deter¬ 
mined  from  model  calibration  fall  within  the  range  of 
experimentally  determined  values  (10-9  to  1(M°  ft/s). 

Specific  Storage 

Very  few  data  on  the  specific  storage  of  the  confin¬ 
ing  beds  in  Southern  Maryland  are  available. 
Hansen  (1977)  presented  the  results  of  consolidation 
tests  on  the  Marlboro  Clay  which  shows  a  specific 
storage  on  the  order  of  1(M  ft-1.  In  the  absence  of 
other  data,  this  value  was  used  in  the  model. 

Thickness  of  Confining  Beds 

The  thickness  of  the  upper  confining  beds  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers  was  esti¬ 
mated  for  each  grid  block  from  figures  6  and  11. 
These  values  were  then  entered  into  the  model  as  a 
matrix. 


Model  Calibration 

Before  a  model  is  used  to  predict  future  water 
levels,  it  is  necessary  to  verify  the  ability  of  the 
model  to  reproduce  historical  water  levels  and 
water-level  changes  due  to  pumpage.  Although 
model  input  represents  the  best  available  data, 
these  parameters  frequently  have  to  be  adjusted 
(calibrated)  in  order  for  the  model  to  reproduce 
historical  water  levels.  In  calibrating  the  Aquia  and 
Piney  Point-Nanjemoy  model,  parameters  that  were 
judged  to  be  relatively  well  known  were  not  varied 
during  calibration.  The  model  parameters  not  varied 
during  calibration  were  starting  head  of  the  Pleis¬ 
tocene  water-table  aquifer,  thickness  of  confining 
beds,  aquifer  transmissivity,  storage  coefficient, 
and  documented  pumpage  by  major  users.  The  ver¬ 
tical  hydraulic  conductivities  of  confining  beds  were 
varied  to  calibrate  the  model. 

The  method  of  calibrating  a  digital  flow  model  is, 
in  large  measure,  determined  by  the  kinds  of  data 
available.  Ideally,  long-term  hydrographs  of  obser¬ 
vation  wells  and  documented  pumpage  would  be 
available  for  each  grid  block  in  the  model.  Such  com¬ 
plete  data,  however,  are  rarely  available.  In 
Southern  Maryland,  enough  water  levels  were  re¬ 
corded  early  in  the  history  of  aquifer  development  to 
estimate  the  prepumping  potentiometric  surfaces  of 
the  Aquia  and  Piney  Point-Nanjemoy  aquifers  (pis. 
3  and  6).  In  addition  to  this  data,  Otton  (1955) 
published  the  potentiometric  surfaces  of  the  Aquia 
and  Piney  Point-Nanjemoy  aquifers  as  they  existed 
in  approximately  1952  (pis.  4  and  7).  Finally,  the 
potentiometric  surfaces  of  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers  were  measured  in  1980  as 
part  of  this  study  (pis.  5  and  8).  A  few  hydrographs 
of  water  levels  in  the  aquifers  of  Southern  Maryland 
are  available  (Weigle,  1970;  Williams,  1979).  These 
hydrographs,  however,  cover  only  a  small  portion  of 
the  aquifer-development  history  and  are  available 
for  only  a  few  locations.  It  was  determined,  there¬ 
fore,  that  the  prepumping,  1952,  and  1980  potentio¬ 
metric  surfaces  were  the  best  data  available  for 
calibrating  the  model. 

The  calibration  of  the  Aquia  and  Piney  Point- 
Nanjemoy  model  was  accomplished  in  three  steps. 
First,  the  model  was  adjusted  so  that  it  could 
reproduce  the  estimated  steady-state  (prepumping) 
potentiometric  surfaces  of  the  aquifers.  This  phase 
of  calibration  was  done  on  a  trial-and-error  basis  by 
adjusting  confining-bed  vertical  hydraulic  conduc¬ 
tivities.  It  was  found  that  relatively  high  values  of 
vertical  hydraulic  conductivity  had  to  be  assigned 
to  the  confining  beds  overlying  the  Piney  Point- 
Nanjemoy  aquifer  in  major  river  valleys.  In  these 
areas  it  had  been  documented  that  Pleistocene 


44 


Table  6. — Pumpage  rates  used  during  the  calibration  period  1890-1980 
[Average  pumpage  in  million  gallons  per  day] 


Pumping  period 

County 

Aquifer 

1 

2 

3 

4 

5 

6 

7 

1890-1942 

1942-46 

1946-52 

1952-66 

1966-76 

1976-79 

1979-80 

St.  Marys  — 

Piney  Point- 
Nanjemoy 

0.14 

0.46 

0.77 

0.69 

0 .79 

0.92 

0.99 

Aquia 

.05 

1.73 

.80 

1.25 

2.00 

2.43 

2.62 

Calvert  - 

Piney  Point- 
Nanjemoy 

- 

.30 

.15 

.37 

.64 

.83 

1.01 

Aquia 

- 

.15 

.08 

.13 

.61 

1.44 

1.81 

Anne  Arundel  — ^ 

Piney  Point- 
Nanjemoy 

- 

- 

- 

- 

- 

- 

- 

Aquia 

- 

- 

.019 

.032 

.10 

.49 

.68 

Charles 

Piney  Point- 
Nanjemoy 

- 

- 

- 

- 

- 

- 

- 

Aquia 

- 

- 

- 

.05 

.10 

.27 

.28 

Prince  Georges 

Piney  Point- 
Nanjemoy 

- 

- 

- 

- 

- 

- 

- 

Aquia 

- 

- 

- 

.05 

.10 

.32 

.34 

Queen  Annes 

Piney  point- 
Nanjemoy 

- 

- 

- 

- 

- 

- 

- 

Aquia 

- 

- 

- 

- 

.1 

.34 

.34 

Talbot 

Piney  Point- 
Nanjemoy 

- 

- 

- 

- 

- 

- 

Aquia 

- 

- 

- 

- 

.80 

.80 

.80 

Dorchester 

Piney  Point- 
Nanjemoy 

.20 

3.91 

3.91 

3.91 

3.91 

3.91 

3.91 

Aquia 

- 

- 

- 

- 

- 

- 

- 

JV  Location  of  pumping  centers  shown  in  figures  11  and  12. 


streams  had  eroded  channels  into  the  confining-bed 
material  and  deposited  more  permeable  silty  sands 
and  gravels.  Because  of  this  geologic  evidence,  it 
was  judged  that  the  higher  values  of  vertical  hy¬ 
draulic  conductivity  in  these  areas  were  justified. 
This  phase  of  calibration  was  considered  complete 
when  the  model  could  reproduce  the  shape  of  the 
prepumping  potentiometric  contours  (pis.  3  and  6) 
and  calculated  water  levels  agreed  with  measured 
water  levels  within  5  ft. 

In  order  to  simulate  the  pumping  history  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers,  pump¬ 
age  was  divided  into  seven  pumping  periods  as 
listed  in  table  6.  The  year  1890  was  chosen  to  start 
the  simulation  period  because  little  pumpage  had 
taken  place  in  these  aquifers  at  that  time  and  it  is 
reasonable  to  assume  near  steady-state  conditions. 
The  divisions  of  pumping  periods  shown  on  table  6 
were  chosen  in  order  to  best  simulate  significant 
changes  in  pumping  patterns.  For  example  the 


heavy  pumping  during  World  War  II,  for  the  period 
1942-46,  was  simulated  as  one  pumping  period. 

The  second  phase  of  calibration  was  to  match  the 
computed  1952  potentiometric  surfaces  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers  with  the 
potentiometric  maps  published  by  Otton  (1952).  It 
was  not  necessary  to  change  model  parameters  dur¬ 
ing  this  phase  of  calibration  in  order  to  obtain  a 
match  that  was  considered  acceptable.  In  all  cases, 
the  calculated  water  levels  agreed  with  measured 
water  levels  within  5  ft. 

The  last  stage  of  calibration  was  to  match  the 
computed  potentiometric  surfaces  for  1980  to  the 
potentiometric  surfaces  of  the  aquifers  as  measured 
in  1980.  Some  changes  in  the  vertical  hydraulic  con¬ 
ductivity  of  the  Piney  Point-Nanjemoy  upper  con¬ 
fining  bed  were  necessary  in  order  to  match  the  com¬ 
puted  and  observed  water  levels  within  a  margin  of 
5  ft.  Therefore,  after  the  changes  were  made  for  the 
last  calibration  period,  a  final  calibration  run  was 
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made  to  insure  that  the  changes  did  not  significant¬ 
ly  alter  the  matches  of  the  first  two  calibration 
periods.  The  final  matches  between  the  computed 
prepumping,  1952,  and  1980  potentiometric  sur¬ 
faces  and  available  water-level  measurements  for 
those  periods  are  shown  on  plates  10  and  11. 

The  level  of  confidence  placed  on  model  predic¬ 
tions  should  be  tempered  by  an  evaluation  of  the 
model  calibration.  In  general,  the  greatest  amount 
of  confidence  in  the  model  is  warranted  where  water- 
level  data  are  plentiful.  Similarly,  less  confidence  in 
the  model  calibration  is  warranted  where  water-level 
data  are  scarce.  Plates  10  and  11  show  the  water- 
level  control  points  which  were  used  for  each  phase 
of  model  calibration.  The  area  where  the  greatest 
amount  of  data  were  available  for  calibration  is  in 
the  Lexington  Park— PN ATS,  Leonardtown,  and 
Piney  Point  areas  of  St.  Marys  County.  The  model 
calibration  in  these  areas,  therefore,  should  be  con¬ 
sidered  the  most  reliable.  Adequate  data  were  also 
available  near  the  Chalk  Point,  Benedict,  Prince 
Frederick,  North  Beach,  and  Shady  Side  areas.  In 
the  areas  near  the  subcrop  zone  of  the  Aquia,  few 
water  levels  were  available  and  the  calibration 
should  be  given  less  confidence.  Also,  areas  on  the 
Eastern  Shore  of  Maryland  had  few  data  available 
and  the  calibration  should  not  be  given  a  high  level 
of  confidence. 

Water  Budget 

The  water  budget  calculated  by  the  digital  model 
during  the  calibration  period  is  of  interest  for  two 
reasons.  First,  it  provides  a  rough  check  on  the 


model  solution  by  indicating  if  the  quantity  of  water 
flowing  into  and  out  of  the  aquifer  system  is  reason¬ 
able.  Secondly,  it  can  provide  estimates  of  the  total 
recharge  to  the  aquifer  system.  These  recharge 
estimates  are  useful  to  water  planners  because  they 
may  reflect  the  yield  potential  of  the  aquifer  system. 

Table  7  summarizes  the  water  budget  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers  (con¬ 
sidered  together)  for  the  90-year  calibration  period 
as  calculated  by  the  model.  The  main  source  of 
water  (recharge)  to  the  aquifer  system  is  leakage 
from  the  overlying  Pleistocene  water-table  aquifer. 
Additional  water  is  derived  from  storage  in  the 
aquifers  and  confining  beds.  Storage  is  treated  as  a 
source  by  this  model  (Trescott,  1975).  Water  is 
discharged  from  the  aquifer  system  by  upward 
leakage  to  the  overlying  Pleistocene  aquifer  and  by 
pumpage. 

The  data  in  table  7  show  that  the  rate  of  recharge 
to  the  aquifer  system  varied  over  the  calibration 
period  from  0.22  to  0.57  in./yr.  These  recharge  rates 
are  consistent  with  rates  to  be  expected  of  confined 
aquifer  systems  (Walton,  1970). 

The  yield  potential  of  an  aquifer  system  is,  in  part, 
limited  by  the  rate  that  it  is  recharged.  It  is  impor¬ 
tant  to  realize,  however,  that  factors  such  as  aquifer 
transmissivity,  proximity  of  aquifer  boundaries, 
and  well-field  spacing  can  also  limit  yield  potential. 
Aquifer  transmissivity  is  particularly  important  in 
this  respect.  If  transmissivity  is  low,  then  draw¬ 
downs  resulting  from  pumpage  can  be  more  of  a 
limiting  factor  to  yield  than  recharge  rates.  Another 
factor  which  complicates  this  type  of  analysis  is 
that  aquifer  recharge  rates  are  not  necessarily  con- 


Table  7.— Model-calculated  water  budget  of  the  Aquia  and  Piney  Point-Nanjemoy  aquifer  system 

for  the  calibration  period 


Pumping 

period 

Sources  — 

Discharges  — 

Percent  difference 
between  sources 
and  discharges 

T  V  2/ 

Leakage— 

Storage 

T  L.  3/‘ 

Leakage— 

Pumpage 

( in. /yr) 

(ft3/s) 

(gal/d) 

(ft3/s) 

(in./yr) 

(ft3/s) 

(gal/d) 

( f  t 3 /s  ) 

1 

0.22 

9.72 

2,160 

3.33X10-3 

0.21 

9.14 

390,000 

0.60 

<  1 

2 

.41 

18.1 

253,600 

0.390 

.19 

8.40 

6,500,000 

9.98 

<  1 

3 

.39 

17.2 

33,700 

.052 

.19 

8.40 

5,800,000 

8.93 

<  1 

4 

.40 

17.7 

33,300 

.051 

.17 

7.51 

6,500,000 

10.0 

<  1 

5 

.45 

19.9 

70,400 

.109 

.16 

7.06 

9,140,000 

14.1 

<  1 

6 

.52 

22.9 

189,000 

.292 

.13 

5.74 

11,500,000 

17.7 

<  1 

7 

.57 

25.2 

230,136 

.355 

.12 

5.30 

12,960,000 

20.0 

<  1 

1 /  Sources  and  discharges  are  considered  £or  the  aquifer  system  as  a  whole. 
2 J  Figures  assume  recharge  area  of  600  mi  ^ 

3/  Figures  assume  discharge  area  of  600  mi  . 
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stant.  If  more  potential  recharge  (precipitation)  is 
available  to  an  aquifer  system  than  is  utilized,  then 
recharge  rates  can  increase  if  pumping  stress  is  in¬ 
creased.  Conversely,  the  recharge  rate  can  decrease 
if  pumping  stress  is  decreased. 

The  complex  relationship  between  aquifer  re¬ 
charge  rates  and  aquifer  yield  can  be  illustrated 
with  the  model-calculated  water  budget  shown  in 
table  7.  During  the  first  pumping  period,  recharge  to 
the  aquifer  system  (leakage)  was  0.22  in./yr  (6.3 
Mgal/d).  The  pumpage  during  this  period  was  only 
0.39  Mgal/d  and,  therefore,  recharge  greatly  ex¬ 
ceeded  pumpage.  In  contrast,  recharge  to  the 
aquifer  system  during  the  last  pumping  period  was 
0.57  in./yr  (16.3  Mgal/d),  while  pumpage  increased 
to  about  13  Mgal/d.  Therefore,  in  the  last  pumping 
period  recharge  still  exceeded  pumpage,  but  the  dif¬ 
ference  was  much  less  than  in  the  first  pumping 
period.  It  is  significant  that  during  the  last  pumping 
period,  large  water-level  declines  occurred  in  some 
places.  These  water-level  declines,  which  are  a  func¬ 
tion  of  aquifer  transmissivity  and  proximity  of 
aquifer  boundaries,  impose  a  practical  limit  on 
aquifer  yield  in  certain  places  even  though  aquifer 
recharge  remains  higher  than  pumpage.  This  sug¬ 
gests  that  aquifer  properties  at  individual  pumping 
locations  are  more  limiting  to  water  production  than 
recharge  rates  in  this  aquifer  system.  From  a  water 
planner’s  point  of  view,  the  amount  of  water-level 
drawdown  considered  acceptable  will  ultimately 
limit  aquifer  yields  rather  than  recharge  rates. 

Simulated  Future  Water  Levels 

The  water  budget  analysis  provided  by  the  model 
suggests  that  aquifer  yield  properties  and  draw¬ 
downs  resulting  from  pumpage  rather  than  the 
recharge  rate  to  the  aquifer  system  are  the  factors 
that  will  limit  the  quantity  of  available  water.  The 
next  step  in  the  model  analysis,  therefore,  was  to 
estimate  aquifer  drawdowns  based  upon  alternative 
scenarios  of  future  pumpage.  The  simulations  con¬ 
sist  of  the  following: 

1.  A  series  of  simulations  using  estimates  of  future 
pumpage  in  Southern  Maryland  for  the  intervals 
(a)  1980-1985,  (b)  1980-1990,  and  (c)  1980-2000. 

2.  A  series  of  simulations  to  predict  water  levels 
in  the  year  2000  assuming  the  Aquia  aquifer  is 
pumped  an  average  of  1.5,  2.0,  2.5,  and  3.0  Mgal/d 
in  the  Lexington  Park— PN ATS  area.  These  sim¬ 
ulations  include  1980  pumpage. 

3.  A  series  of  simulations  designed  to  show  the  ef¬ 
fects  of  withdrawing  1  Mgal/d  from  northern  St. 
Marys  County  for  10  years.  Two  pumping 
schemes  were  simulated:  (a)  1  Mgal/d  of  pumpage 


from  a  single  well  field;  and  (b)  1  Mgal/d  of  pump¬ 
age  spaced  among  four  well  fields. 

4.  A  series  of  simulations  designed  to  show  draw¬ 
down  in  the  Piney  Point-Nanjemoy  aquifer  due  to 
pumpage  in  the  Aquia  aquifer.  Two  simulations 
are  presented  that  show  (a)  drawdown  in  the 
Piney  Point-Nanjemoy  due  to  a  1  Mgal/d  with¬ 
drawal  from  the  Aquia  in  the  Lexington  Park- 
PNATS  area,  and  (b)  drawdown  in  the  Piney 
Point-Nanjemoy  due  to  a  2  Mgal/d  withdrawal 
from  the  Aquia  in  the  Lexington  Park-PNATS 
area. 

5.  A  series  of  simulations  designed  to  show  draw¬ 
down  effects  of  a  hypothetical  0.5-Mgal/d  well 
field  in  southern  Anne  Arundel  County.  Three 
different  sites  were  simulated  at  this  pumping 
rate  for  a  period  of  10  years.  These  sites  were 
located  near  (a)  Shady  Side,  (b)  Fairhaven,  and  (c) 
Bristol. 

6.  A  simulation  of  a  0.5  Mgal/d  withdrawal  from  the 
Aquia  aquifer  at  the  Chalk  Point  Power  Plant  in 
Prince  Georges  County. 

7.  A  hypothetical  simulation  which  assumes  that 
all  appropriated  users  of  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers  daily  pump  their  max¬ 
imum  Groundwater  Appropriation  Permit  (GAP) 
allocations  for  a  period  of  10  years. 

8.  A  simulation  which  assumes  a  pumpage  of  0.5 
Mgal/d  from  the  Aquia  aquifer  in  Chesapeake 
Beach. 

Each  of  these  pumping  scenarios  was  selected  and 
simulated  in  response  to  information  requests  from 
the  planning  departments  of  St.  Marys,  Calvert,  and 
Anne  Arundel  Counties.  They  are  designed  to  il¬ 
lustrate  the  yield  potential  of  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers  in  areas  where  pumpage  is 
expected  to  increase  in  the  future.  Because  of  the 
many  simplifications  made  of  the  natural  aquifer 
system  during  construction  of  the  digital  model,  the 
drawdown  and  water-level  predictions  made  should 
not  be  considered  exact.  Instead,  they  should  be 
considered  the  best  estimates  which  can  be  made, 
given  the  current  available  data  and  technology. 

Projected  Water  Levels  Based  on  Estimates 
of  Future  Withdrawal 

A  series  of  transient  simulations  of  the  Aquia  and 
Piney  Point-Nanjemoy  aquifers  were  made  using 
estimates  of  future  withdrawals.  These  estimates 
(table  8)  assume  that  future  pumpage  will  be  direct¬ 
ly  proportional  to  the  projected  population  growth. 
For  example,  if  the  projected  population  increase  is 
10  percent,  pumpage  is  assumed  to  increase  10  per¬ 
cent.  These  estimates  also  assumed  that  all  in- 
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Table  8.— Projected  rates  of  future  pumpage  from  the  Aquia  and  Piney  Point-Nanjemoy  aquifers 

based  on  county  population  growth  estimates 


County 

County  planning  area 

Projected  population  — ^ 

Aquifer 

2/ 

Projected  pumping  rates  —  (gal/d) 

1985 

1990 

2000 

1985 

1990 

2000 

St.  Marys 

Pine  Hill  Run 

31,026 

33,527 

39,194 

Aquia 

1,740,000 

1,890,000 

1,970,000 

Piney  Point-Nanjemoy 

700,000 

750,000 

980,000 

Leonard  town 

6,093 

6,926 

8,817 

Aquia 

380,000 

440,000 

550,000 

Piney  Point-Nanjemoy 

0 

0 

0 

Luck land  Run 

6,620 

8,445 

12,582 

Aquia 

310,000 

390,000 

590,000 

Piney  Point-Nanjemoy 

0 

0 

0 

Flood  Creek 

1,733 

1,974 

2,522 

Aquia 

10,000 

11,000 

16,000 

Piney  Point-Nanjemoy 

20,000 

23,000 

28,000 

Piney  Point 

3,100 

3,286 

3,708 

Aquia 

180,000 

190,000 

210,000 

Piney  Point-Nanjemoy 

62,000 

66,000 

74,000 

Lake  Conoy 

2,056 

2,223 

2,603 

Aquia 

0 

0 

0 

Piney  Point-Nanjemoy 

130,000 

140,000 

116,000 

Carrol  Pond 

1,426 

1,593 

1,973 

Aquia 

0 

0 

0 

Piney  Point-Nanjemoy 

160,000 

180,000 

220,000 

Manor  Run 

2,530 

3,030 

4,165 

Aquia 

32,000 

39,000 

53,000 

Piney  Point-Nanjemoy 

16,000 

19,000 

27,000 

Indian  Creek 

6,965 

8,800 

12,937 

Aquia 

280,000 

360,000 

520,000 

Piney  Point-Nanjemoy 

0 

0 

0 

Dukeharts  Creek 

7,871 

8,880 

11,170 

Aquia 

37,000 

42,000 

52,000 

Piney  Point-Nanjemoy 

0 

0 

0 

Calvert 

A 

13,000 

15,500 

20,500 

Aquia 

684,800 

729,600 

854,400 

Piney  Point-Nanjemoy 

374,400 

425,600 

582,400 

B 

11,150 

13,100 

17,300 

Aquia 

599,200 

638,400 

747,600 

Piney  Point-Nanjemoy 

327,600 

372,400 

509,600 

C 

16,450 

19,800 

26,200 

Aquia 

856,000 

912,000 

1,068,000 

Piney  Point-Nanjemoy 

468,000 

532,000 

728,000 

Anne  Arunde 1 

8  (Rural) 

2/ 

2/ 

2/ 

Aquia 

616,000 

805,000 

1,043,000 

Piney  Point-Nanjemoy 

0 

0 

0 

6  (Annapolis) 

2/ 

2/ 

2/ 

Aquia 

264,000 

345,000 

447,000 

Piney  Point-Nanjemoy 

0 

0 

0 

Char les 

3/ 

2/ 

2/ 

y 

Aquia 

310,000 

340,000 

370,000 

Piney  Point-Nanjemoy 

0 

0 

0 

Prince  Georges 

3/ 

2/ 

2/ 

y 

Aquia 

410,000 

490,000 

590,000 

Piney  Point-Nanjemoy 

0 

0 

0 

_1/  Population  projections  from  Calvert  County  Comprehensive  Water  and  Sewerage  Plan,  1978  update,  and 
St.  Marys  County  Comprehensive  Water  and  Sewerage  Plan,  1977  update. 

2, 1  Projected  pumping  rate  assumes  that  future  pumpage  is  directly  proportional  to  population  growth. 

3/  Population  projections  by  planning  district  not  available. 


creases  of  future  pumpage  in  Southern  Maryland 
would  come  from  the  Aquia  and  Piney  Point- 
Nanjemoy  aquifers.  Because  the  Maryland  State 
Water  Resources  Administration  has  mandated 
that  the  Piney  Point-Nanjemoy  aquifer  be  reserved 
for  domestic  and  small  industrial  users,  all  increases 


in  municipal  and  large  industrial  users  were  as¬ 
sumed  to  come  from  the  Aquia  aquifer.  The  pump¬ 
ing  centers  used  in  this  simulation  are  those  located 
in  figure  7  and  12. 

Figures  21  and  22  show  the  predicted  poten- 
tiometric  surface  of  the  Aquia  and  Piney  Point- 
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Figure  22.— Simulated  1985  potentiometric  surface  of  the  Piney  Point-Nanjemoy  aquifer  and 
simulated  drawdown  from  1980  based  on  estimates  of  future  pumpage. 


Nanjemoy  aquifers  in  1985.  Also  shown  are  the 
drawdowns  from  the  1980  potentiometric  surface. 
The  maximum  head  declines  for  the  Aquia  aquifer 
are  predicted  to  occur  in  the  Lexington 
Park— PNATS  area,  in  northern  St.  Marys  County, 
and  near  Prince  Frederick  in  Calvert  County.  In  the 
Lexington  Park— PNATS  area,  declines  of  8  to  10  ft 
are  predicted.  In  northern  St.  Marys  County,  where 
presently  rapid  population  growth  is  taking  place, 
10  to  15  ft  of  head  decline  is  predicted.  Near  Prince 
Frederick,  approximately  8  to  10  ft  of  head  decline  is 
predicted.  The  maximum  head  declines  for  the  Piney 
Point-Nanjemoy  aquifer  are  predicted  to  occur  in 
the  Lexington  Park-PNATS  area.  Declines  in  this 
area  are  predicted  to  be  about  6  ft  for  the  simulated 
period. 

Figures  23  and  24  show  the  Aquia  and  Piney 
Point-Nanjemoy  potentiometric  surfaces  predicted 
for  1990  based  on  estimated  future  withdrawals. 
Also  shown  are  the  predicted  water-level  declines 
from  the  1980  potentiometric  surfaces.  For  the 
Aquia  aquifer,  the  greatest  drawdowns  are  pre¬ 
dicted  to  occur  in  the  Lexington  Park-PNATS  area, 
northern  St.  Marys  County,  and  Prince  Frederick. 
For  the  Piney  Point-Nanjemoy  aquifer,  the  greatest 
drawdowns  are  predicted  to  occur  in  the  Lexington 
Park-PNATS  area,  near  Prince  Frederick,  and  in  the 
North  Beach  area. 

Figures  25  and  26  show  the  predicted  potentio¬ 
metric  surfaces  of  the  Aquia  and  Piney  Point- 
Nanjemoy  aquifers  in  the  year  2000.  Also  shown  are 
drawdowns  from  the  1980  potentiometric  surfaces. 

This  series  of  simulations  suggests  that  water 
levels  in  the  Aquia  and  Piney  Point-Nanjemoy 
aquifers  will  decline  in  the  future  in  many  areas.  The 
predicted  maximum  rate  of  decline  for  the  Aquia 
aquifer  is  approximately  1.5  ft/yr  near  Lexington 
Park.  The  predicted  maximum  rate  of  decline  for  the 
Piney  Point-Nanjemoy  aquifer  is  approximately  1 
ft/yr  near  Lexington  Park.  The  assumption  that 
pumpage  will  be  proportional  to  population  growth 
is  probably  conservative  and  higher  pumpages  than 
the  projections  used  in  these  simulations  are  possi¬ 
ble.  In  this  case,  water-level  declines  greater  than 
those  predicted  will  occur. 

This  series  of  simulations  predicts  that  by  the 
year  2000,  the  potentiometric  surface  of  the  Aquia 
aquifer  near  Lexington  Park  will  be  more  than  100  ft 
below  sea  level.  The  actual  pumping  levels  of  pro¬ 
duction  wells  will  be  much  lower,  probably  on  the 
order  of  200  to  300  ft  below  sea  level.  These 
predicted  low  water  levels  may  make  production  of 
water  from  some  domestic  wells  in  this  area  (which 
typically  are  designed  to  lift  water  250  to  300  ft)  im¬ 
practical. 


Projected  Water  Levels  Based  on  Hypothetical 

Withdrawal  Rates  near  Lexington  Park-PNATS 

The  Lexington  Park-PNATS  area  of  St.  Marys 
County  is  presently  the  heaviest  user  of  ground 
water  from  the  Aquia  aquifer  in  Southern  Maryland. 
It  is  important  for  planners  to  be  able  to  estimate 
the  amount  of  water  which  is  available  from  the 
Aquia  aquifer  in  this  area.  To  aid  planners  in  this 
decision-making  process,  a  series  of  transient  simu¬ 
lations  was  done  for  the  Lexington  Park-PNATS 
area  that  simulates  different  pumping  rates  for  the 
period  1980-2000.  Included  in  these  simulations  is 
the  average  1980  pumpage  (tables  4  and  5).  The 
pumping  center  simulated  is  location  26  in  figure  7. 

The  first  simulation  specified  a  total  withdrawal 
rate  of  1.5  Mgal/d  in  the  Lexington  Park-PNATS 
area  for  20  years.  The  approximate  pumpage  rate  in 
this  area  for  May  1980  was  1.5  Mgal/d,  so  that  no 
net  increase  in  pumpage  is  simulated.  Figure  27 
shows  the  predicted  drawdown  from  1980  to  the 
year  2000.  According  to  this  simulation,  1.5  Mgal/d 
of  pumpage  will  result  in  about  10  ft  of  drawdown  in 
the  Lexington  Park-PNATS  area.  The  drawdown  in 
northern  St.  Marys  County  is  due  to  the  1980  rate  of 
domestic  pumpage. 

The  second  simulation  assumed  a  total  with¬ 
drawal  in  the  Lexington  Park-PNATS  area  of  2.0 
Mgal/d  for  20  years,  which  is  a  net  increase  of  0.5 
Mgal/d  from  1980.  Figure  28  shows  the  predicted 
drawdown  from  the  1980  potentiometric  surface. 
The  maximum  predicted  drawdown  for  the  year 
2000  is  approximately  40  ft. 

The  third  simulation  in  this  series  assumed  a  total 
withdrawal  in  the  Lexington  Park-PNATS  area  of 
2.5  Mgal/d  for  20  years,  which  is  a  net  increase  of  1.0 
Mgal/d  from  1980.  Figure  29  shows  the  predicted 
drawdown  from  the  1980  potentiometric  surface  in 
the  Aquia  aquifer.  This  drawdown  is  predicted  to  be 
approximately  60  ft  during  the  simulation  period  for 
the  center  of  the  cone  of  depression. 

The  last  simulation  in  this  series  specified  an 
average  withdrawal  rate  of  3.0  Mgal/d  from  the 
Aquia  in  the  Lexington  Park-PNATS  area,  which 
is  a  net  increase  of  1.5  Mgal/d  from  1980.  The  pre¬ 
dicted  drawdown  from  1980  is  shown  in  figure  30. 
This  drawdown  is  predicted  to  be  approximately  90 
ft  during  the  simulation  period  for  the  center  of  the 
cone  of  depression. 

This  series  of  simulations  provides  water  planners 
with  drawdown  estimates  in  the  Aquia  aquifer  for 
various  rates  of  pumpage  in  the  Lexington  Park- 
PNATS  area.  These  simulations  are  designed  to 
help  planners  estimate  the  approximate  quantity  of 
water  available  in  this  area.  The  quantity  of  water 
available  ultimately  depends  on  how  much  draw- 
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Figure  23.— Simulated  1990  potentiometric  surface  of  the  Aquia  aquifer  and  simulated  drawdown 
from  1980  based  on  estimates  of  future  pumpage. 
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Figure  24.— Simulated  1990  potentiometric  surface  of  the  Piney  Point-Nanjemoy  aquifer  and 
simulated  drawdown  from  1980  based  on  estimates  of  future  pumpage. 
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Figure  26.— Simulated  2000  potentiometric  surface  of  the  Piney  Point-Nanjemoy  aquifer  and 
simulated  drawdown  from  1980  based  on  estimates  of  future  pumpage. 
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Figure  27.— Simulated  drawdown  from  the  1980  potentiometric  surface  in  the 
Aquia  aquifer  based  on  20  years  of  1.5  Mgal/d  pumpage  from 
the  Lexington  Park— PNATS  area. 
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Figure  28.— Simulated  drawdown  from  the  1980  potentiometric  surface  in  the 
Aquia  aquifer  based  on  20  years  of  2.0  Mgal/d  pumpage  in  the 
Lexington  Park-PNATS  area. 
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Aquia  aquifer  based  on  20  years  of  2.5  Mgal/d  pumpage  in 
the  Lexington  Park-PNATS  area. 
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Aquia  aquifer  based  on  20  years  of  3.0  Mgal/d  pumpage  in 
the  Lexington  Park-PNATS  area. 
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down  is  judged  to  be  acceptable.  If,  for  example,  it  is 
decided  that  an  additional  average  drawdown  (draw¬ 
down  in  a  grid  block)  of  more  than  60  ft  would  be 
undesirable,  the  effective  limit  of  aquifer  yield 
would  be  approximately  2.0  to  2.5  Mgal/d  in  this 
area. 

It  should  be  reemphasized  at  this  point  that  the 
model  calculates  average  drawdowns  for  each  grid 
block  (1  mi2  in  the  Lexington  Park-PNATS  area). 
Pumping  levels  in  individual  wells  within  this  grid 
block  are  not  calculated  by  the  model.  These  in¬ 
dividual  pumping  levels,  which  are  dependent  on 
pumping  rates,  well  diameter,  and  well  efficiency 
will  be  greater  than  the  average  drawdowns  given 
by  the  model. 

Simulation  of  1  Mgal/d  Withdrawal  from 
Northern  St.  Marys  County 

Northern  St.  Marys  County  is  presently  a  rapid- 
growth  area.  It  has  been  estimated  by  county  plan¬ 
ning  officials  that  in  the  next  20  years  the  popula¬ 
tion  of  this  area  will  more  than  double  (table  8).  This 
increasing  population  will  increase  demands  for 
water.  The  Aquia  and  Piney  Point-Nanjemoy  model 
was  used  to  simulate  two  pumping  scenarios  for  this 
region.  For  these  simulations,  it  was  assumed  that 
growth  in  this  area  would  result  in  a  demand  for  an 
additional  1  Mgal/d.  The  first  simulation  placed  a  1 
Mgal/d  well  field  at  the  present  location  of  the  St. 
Marys  County  Metropolitan  Commission  (METCOM) 
Rolling  Acres  (table  4,  no.  41)  well  field  in  northern 
St.  Marys  County.  The  second  scenario  again 
simulated  withdrawal  of  1  Mgal/d,  but  spread  the 
pumpage  between  the  Birch  Manor,  Country  Lakes, 
King  and  Kennedy,  and  Rolling  Acres  (table  4,  nos. 
19,  21,  36  and  41)  well  fields  in  northern  St.  Marys 
County.  Average  1980  pumpage  was  also  used  for 
these  simulations  so  that  the  1  Mgal/d  is  an  increase 
above  1980  pumpage. 

Figure  31  shows  the  predicted  drawdown  of  the 
Aquia  aquifer  based  upon  a  1 -Mgal/d  withdrawal 
from  the  METCOM  Rolling  Acres  wells  for  a  period 
of  10  years.  At  the  end  of  the  simulation  period,  the 
drawdown  in  the  pumping  node  had  declined  slight¬ 
ly  more  than  100  ft. 

Figure  32  shows  the  predicted  drawdown  after  10 
years  of  pumping  the  Rolling  Acres,  Birch  Manor, 
Country  Lakes,  and  King  and  Kennedy  wells  at  a 
combined  rate  of  1  Mgal/d.  In  this  case,  the  largest 
predicted  drawdown  in  one  node  is  about  45  ft. 

A  comparison  of  figures  31  and  32  shows  the  ef¬ 
fects  of  spreading  pumpage  between  well  fields  as 
opposed  to  withdrawing  all  the  required  water  from 
one  well  field.  The  impacts  as  far  as  decline  of  water 


level  is  concerned  are  less  if  pumpage  is  spread  out. 
Careful  spacing  of  well  fields  can  minimize  draw¬ 
down  impacts  of  pumpage. 

Simulation  of  Drawdowns  in  the  Piney  Point- 
Nanjemoy  Aquifer  due  to  Aquia  Pumpage 

In  Southern  Maryland,  the  Maryland  State  Water 
Resources  Administration  has  mandated  that  the 
Piney  Point-Nanjemoy  aquifer  be  reserved  for  use 
by  homeowners  and  small  industries.  Large  in¬ 
dustry  and  municipal  wells  therefore  have  been 
developed  almost  exclusively  in  the  Aquia  aquifer. 
One  question  that  has  arisen  over  this  practice  is 
what  effects  do  drawdowns  in  the  Aquia  aquifer 
have  on  heads  in  the  Piney  Point-Nanjemoy  aquifer. 
It  has  been  suggested  that  drawdowns  in  the  Aquia 
could  produce  drawdowns  in  the  Piney  Point- 
Nanjemoy. 

It  is  possible  to  address  this  type  of  question  with 
the  digital  model  developed  for  this  study.  Tran¬ 
sient  leakage  between  aquifers  is  considered  by  this 
model  and  it  is  therefore  possible  to  simulate  pump¬ 
age  in  one  aquifer  and  predict  resulting  drawdowns 
in  adjacent  aquifers. 

A  model  run  designed  to  address  this  question 
simulated  a  1 -Mgal/d  withdrawal  from  the  Aquia 
aquifer  in  the  Lexington  Park-PNATS  area.  The 
simulation  was  started  with  steady-state  heads  in 
both  the  Aquia  and  Piney  Point-Nanjemoy  aquifers 
and  no  other  pumpage  was  simulated.  These  start¬ 
ing  conditions  were  specified  to  eliminate  drawdown 
effects  of  other  pumpage  and  to  isolate  the  effects  of 
vertical  leakage  on  the  Piney  Point-Nanjemoy 
aquifer.  No  pumpage  was  simulated  in  the  Piney 
Point-Nanjemoy  aquifer. 

Figure  33  shows  the  calculated  drawdowns  in  the 
Piney  Point-Nanjemoy  aquifer  due  to  withdrawal  of 
1  Mgal/d  from  the  Aquia  for  a  period  of  10  years. 
This  pumpage  resulted  in  a  drawdown  of  about  90  ft 
in  the  Aquia  for  this  time  interval.  The  correspond¬ 
ing  drawdown  in  the  Piney  Point-Nanjemoy  aquifer 
is  about  4  ft  for  the  same  time  interval. 

An  additional  simulation  was  made  that  specified 
a  2-Mgal/d  withdrawal  from  the  Aquia  in  the  Lex¬ 
ington  Park-PNATS  area.  Again  the  simulation 
was  started  under  initial  steady-state  conditions 
and  only  withdrawal  from  the  Aquia  was  specified. 
The  simulation  period  was  again  10  years.  Figure  34 
shows  the  drawdown  in  the  Piney  Point-Nanjemoy 
aquifer  which  results  from  the  scenario.  The  Aquia 
shows  a  drawdown  of  about  180  ft.  The  correspond¬ 
ing  drawdown  in  the  Piney  Point-Nanjemoy  over 
the  same  time  period  is  about  7  ft. 
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Figure  31. — Simulated  drawdown  from  1980  in  the  Aquia  aquifer  based  on  10  years  of  1.0-Mgal/d 
pumpage  in  the  northern  area  of  St.  Marys  County  from  one  well  field. 


61 


76*45' 


76°  30 1 


76°  15' 


Figure  32.— Simulated  drawdown  from  1980  in  the  Aquia  aquifer  based  on  10  years  of  1.0-Mgal/d 
pumpage  in  northern  St.  Marys  County  from  four  well  fields. 
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Figure  33.— Simulated  drawdown  from  a  steady-state  potentiometric  surface  in  the  Piney  Point-Nanjemoy 
aquifer  based  on  10  years  of  1.0  Mgal/d  pumpage  from  the  Aquia  aquifer  near  Lexington  Park 
-PNATS. 
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Figure  34.— Simulated  drawdown  from  a  steady-state  potentiometric  surface  in  the  Piney  Point-Nanjemoy 
aquifer  based  on  10  years  of  2.0  Mgal/d  pumpage  from  the  Aquia  aquifer  near  Lexington  Park 


-PNATS. 
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These  simulations  show  that  major  drawdowns  in 
the  Aquia  can  have  a  drawdown  impact  on  the  Piney 
Point-Nanjemoy  aquifer.  These  drawdowns  in  the 
Piney  Point-Nanjemoy  are,  however,  rather  small. 
It  is  not  likely,  therefore,  that  significant  draw¬ 
downs  in  the  Piney  Point-Nanjemoy  have  been  caused 
by  Aquia  pumpage.  It  should  be  noted  that  these 
simulations  assume  an  undisturbed  confining  bed 
between  the  Aquia  and  the  Piney  Point-Nanjemoy 
aquifers  and  that  there  is  no  artificial  connection 
between  the  two  aquifers.  It  is  possible  for  im¬ 
properly  grouted  wells  to  provide  a  conduit  of  water 
between  aquifers.  Significant  artificial  connection  of 
the  Aquia  and  Piney  Point-Nanjemoy  aquifers  caused 
by  poorly  grouted  wells,  or  multiple-screened  wells, 
could  increase  the  drawdown  impacts  of  one  aquifer 
on  the  other. 

Simulation  of  a  Hypothetical  0.5-Mgal/d  Well  Field 
in  Southern  Anne  Arundel  County 

Most  of  the  present  water  usage  in  southern  Anne 
Arundel  County  is  by  homeowners'  wells  and  by 
small  industries.  The  Aquia  aquifer  in  southern 
Anne  Arundel  County  is  relatively  shallow  and  is  a 
good  producer  of  ground  water.  Because  of  this,  the 
majority  of  wells  in  this  area  are  screened  in  the 
Aquia. 

At  the  present  time,  there  is  no  municipal  water 
service  planned  for  southern  Anne  Arundel  County. 
However,  such  a  system  may  become  desirable 
sometime  in  the  future  near  one  or  more  of  the 
population  centers  such  as  Shady  Side,  Fairhaven, 
or  Bristol.  For  this  reason,  three  simulations  were 
made  of  hypothetical  well  fields  near  each  of  these 
towns. 

Figure  35  shows  the  predicted  drawdown  for  a  0.5 
Mgal/d  well  field  near  Shady  Side.  Domestic  and 
municipal  pumpage  for  1980  is  included  in  this  sim¬ 
ulation.  The  calculated  heads  in  this  simulation 
show  a  drawdown  of  about  20  ft  for  the  10-year  sim¬ 
ulation  period  from  the  1980  potentiometric  surface. 

Figure  36  shows  the  predicted  drawdown  from  the 
1980  potentiometric  surface  for  a  0.5  Mgal/d  well 
field  near  Bristol  in  southern  Anne  Arundel  County. 
Domestic  and  municipal  pumpage  for  1980  is  in¬ 
cluded  in  this  simulation  and  the  period  simulated  is 
10  years.  Calculated  drawdowns  due  to  the  well  field 
are  about  26  ft  for  the  simulation  period. 

The  final  simulation  of  pumpage  in  southern  Anne 
Arundel  County  was  near  Fairhaven.  Again,  a  0.5 
Mgal/d  well  field  was  simulated  for  a  period  of  10 
years.  The  results  of  this  simulation  are  shown  in 
figure  37.  The  maximum  drawdowns  from  the  1980 
potentiometric  surface  are  about  28  ft. 


The  moderate  drawdowns  produced  by  the 
simulated  pumpage  near  these  towns  in  Anne 
Arundel  County  suggest  that  municipal  well  fields 
in  these  locations  are  feasible.  It  should  be  noted, 
however,  that  the  predicted  cones  of  depression  for 
the  Shady  Side  and  Fairhaven  simulations  may  in¬ 
tersect  the  subsurface  truncation  of  the  Aquia 
aquifer.  This  introduces  the  possibility  that  chlo¬ 
rides  may  enter  the  Aquia  aquifer  due  to  leakage 
from  the  Chesapeake  Bay.  Monitoring  chloride  con¬ 
centrations  of  Aquia  water  in  these  areas  would  be  a 
prudent  precaution.  It  is  interesting  to  note  that  the 
simulated  0.5  Mgal/d  simulation  near  Shady  Side 
produces  less  drawdown  than  the  0.5  Mgal/d  simula¬ 
tion  near  Fairhaven.  This  reflects  the  increasing 
transmissivity  of  the  Aquia  aquifer  north  of 
Fairhaven. 

Simulation  of  a  0.5  Mgal/d  Withdrawal  at  the 
Chalk  Point  Power  Plant,  Prince  Georges  County 

The  Chalk  Point  Power  Plant,  which  is  fired  by 
fossil  fuels,  currently  is  a  heavy  user  of  Cretaceous 
aquifers  (Mack,  1976).  This  heavy  usage  has  led  to 
large  drawdowns  in  the  Magothy  aquifer  in  par¬ 
ticular.  The  Aquia  aquifer  (but  not  the  Piney  Point- 
Nanjemoy  aquifer)  is  present  at  Chalk  Point  and  its 
ability  to  supply  water  to  the  power  plant  is  of  in¬ 
terest  to  water  planners.  For  this  reason,  a  model 
run  which  simulated  the  withdrawal  of  0.5  Mgal/d 
from  the  Aquia  at  Chalk  Point  for  a  period  of  10 
years  was  made.  The  pumpage  for  1980  is  included 
in  this  simulation.  Figure  38  shows  a  simulated 
drawdown  of  45  ft  from  the  1980  potentiometric  sur¬ 
face  after  10  years  of  pumping. 

This  simulation  indicates  that  0.5  Mgal/d  of 
Aquia  pumpage  from  the  Chalk  Point  Power  Plant 
would  have  a  relatively  large  impact  on  water  levels 
in  the  area.  This  is  partly  because  of  the  lower 
transmissivity  in  the  Aquia  aquifer  west  of  Calvert 
County,  and  partly  because  of  domestic  and  small 
industrial  use  in  the  same  area. 

Simulation  of  Maximum  GAP  Pumpage 

The  Maryland  State  Water  Resources  Adminis¬ 
tration  presently  has  a  policy  that  all  users  of 
ground  water  who  pump  more  than  10,000  gal/d 
must  apply  for  a  Ground- Water  Appropriation  Per¬ 
mit  (GAP).  These  GAP  permits  specify  an  average 
amount  and  a  maximum  amount  of  water  that  each 
user  may  consume.  The  average  is  defined  to  be  the 
daily  pumpage  averaged  over  a  year,  while  the  maxi¬ 
mum  is  the  daily  average  for  the  highest  use  month. 
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In  practice,  ground-water  users  seldom  pump  the 
maximum  amount  of  water  allocated  to  them.  How¬ 
ever,  to  illustrate  a  high  stress  situation  it  is  in¬ 
structive  to  estimate  the  effects  on  water  levels 
assuming  each  user  actually  did  pump  their  maxi¬ 
mum  appropriation  daily  over  a  10-year  simulation 
period.  This  is  a  hypothetical  situation  and  does  not 
simulate  the  effects  of  presently  allocated  ground- 
water  in  the  study  area. 

Figures  39  and  40  show  the  predicted  drawdowns 
from  the  1980  potentiometric  surfaces  of  the  Aquia 
and  Piney  Point-Nanjemoy  aquifers  assuming  each 
holder  of  a  1980  GAP  permit  (tables  4  and  5)  pumped 
their  maximum  allotment  of  ground  water  for  a 
period  of  10  years.  Results  suggest  that  large  draw¬ 
downs  would  probably  occur  in  the  Aquia  aquifier 
near  Lexington  Park-PNATS,  Leonardtown,  Piney 
Point,  the  Calvert  Industrial  Park,  and  Cove  Point. 
In  the  Piney  Point-Nanjemoy  aquifer,  there  are 
small  drawdowns  caused  in  St.  Mary’s  County.  This 
simulation  suggests  that  near  Piney  Point,  Leonard¬ 
town,  and  Lexington  Park,  each  additional  appropri¬ 
ation  of  Aquia  ground  water  should  be  carefully 
considered.  In  other  areas  of  Southern  Maryland, 
the  Aquia  aquifer  does  not  appear  overly  stressed. 
The  simulation  of  maximum  GAP  withdrawals  in 
the  Piney  Point-Nanjemoy  aquifer  shows  that  no 
drawdowns  greater  than  about  20  ft  would  occur. 
This  indicates  that  the  Piney  Point-Nanjemoy  aqui¬ 
fer  is  not  presently  over-allocated. 

Simulation  of  a  0.5  Mgal/d  Withdrawal  at 
Chesapeake  Beach,  Calvert  County 

The  community  of  Chesapeake  Beach  in  Calvert 
County  is  presently  evaluating  the  feasibility  of  con¬ 


structing  a  municipal  water  system  pumping  from 
the  Aquia  aquifer.  In  view  of  this,  the  Calvert  County 
Planning  Department  requested  that  a  simulation 
be  performed  specifying  a  pumpage  of  0.5  Mgal/d 
combined  with  estimated  future  pumpage  (table  8) 
for  a  period  of  5  years  at  Chesapeake  Beach. 

Figure  41  shows  the  drawdowns  which  resulted 
from  this  simulation.  Approximately  20  ft  of  draw¬ 
down  is  predicted  for  the  simulation  period.  It 
should  be  noted  that  transient  flow  was  occurring  at 
the  end  of  the  simulation  period.  This  indicates  that 
the  water  levels  at  Chesapeake  Beach  were  still 
dropping  at  the  end  of  the  simulation  period.  Thus, 
at  the  end  of  10  years,  the  drawdown  would  be 
greater  than  those  shown  in  figure  38,  which  was  for 
a  period  of  5  years.  The  simulation  performed  for 
Fairhaven  in  Anne  Arundel  County  (5  miles  from 
Chesapeake  Beach)  indicated  a  drawdown  of  28  ft 
for  a  0.5  Mgal/d  withdrawal  for  10  years.  It  is  likely 
that  after  10  years  the  drawdown  at  Chesapeake 
Beach  would  also  be  about  28  ft. 

This  simulation  indicates  that  a  municipal  water 
system  pumping  from  the  Aquia  aquifer  at  Chesa¬ 
peake  Beach  would  produce  moderate  drawdowns. 
This  simulation  also  indicates  that  the  edge  of  the 
drawdown  cone  may  possibly  intersect  the  subsur¬ 
face  truncation  of  the  Aquia  aquifer  under  the  Ches¬ 
apeake  Bay.  There  is  a  possibility  that  this  could 
result  in  chlorides  moving  into  the  aquifer  from  the 
Bay.  The  likelihood  of  this  possibility  cannot  be 
evaluated  within  the  scope  of  this  project.  However, 
it  would  be  a  reasonable  precaution  to  monitor  the 
Chesapeake  Beach  municipal  system  for  possible 
chloride  contamination. 
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Figure  35. — Simulated  drawdown  from  the  1980  potentiometric  surface  in  the  Aquia  aquifer  based  on 
0.5  Mgal/d  pumpage  for  10  years  near  Shady  Side,  Anne  Arundel  County. 
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Figure  36.— Simulated  drawdown  from  the  1980  potentiometric  surface  in  the  Aquia  aquifer  based  on 
0.5  Mgal/d  pumpage  for  10  years  near  Bristol,  Anne  Arundel  County. 
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Figure  37.— Simulated  drawdown  from  the  1980  potentiometric  surface  in  the  Aquia  aquifer  based  on 
0.5  Mgal/d  pumpage  for  10  years  near  Fairhaven,  Anne  Arundel  County. 
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Figure  38.— Simulated  drawdown  from  the  1980  potentiometric  surface  of  the  Aquia  aquifer  based  on 
0.5  Mgal/d  pumpage  for  10  years  at  the  Chalk  Point  Power  Plant. 
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Figure  39.— Simulated  drawdown  from  the  1980  potentiometric  surface  of  the  Aquia  aquifer  assuming  all  holders 
of  1980  Ground-Water  Appropriation  Permits  pump  their  maximum  allocations  for  10  years. 
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Figure  40.— Simulated  drawdown  from  the  1980  potentiometric  surface  of  the  Piney  Point-Nanjemoy  aquifer  as 
suming  all  holders  of  1980  Ground-Water  Appropriation  Permits  pump  their  maximum  allocations  fo 
10  years. 
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Figure  41.— Simulated  drawdown  from  the  1980  potentiometric  surface  in  the  Aquia  aquifer  based  on  0.5  Mgal/d 
pumpage  for  5  years  in  the  Chesapeake  Beach  area,  Calvert  County. 
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GROUND  WATER  GEOCHEMISTRY 


Introduction 

The  chemical  quality  of  water  produced  by  an 
aquifer  system  is  an  important  consideration  in  de¬ 
termining  its  potential  uses.  In  many  communities 
ground  water  is  utilized  by  a  wide  spectrum  of  citi¬ 
zens,  some  of  whom  may  have  widely  different 
water-quality  requirements.  It  is  appropriate,  there¬ 
fore,  to  include  a  discussion  on  ground-water  quality 
in  a  description  of  an  aquifer  system.  It  is  also  ap¬ 
propriate  to  discuss  the  chemical  and  physical  pro¬ 
cesses  which  control  the  chemical  character  of  the 
water.  For  the  Aquia  and  Piney  Point-Nanjemoy 
aquifers,  this  is  of  particular  interest  because  the 


chemical  character  of  their  water  changes  markedly 
from  place  to  place.  Aquia  water  in  Anne  Arundel 
County,  for  example,  has  a  chemical  character  which 
is  significantly  different  from  Aquia  water  in  St. 
Mary’s  County.  This  change  in  chemical  character 
can  be  directly  related  to  chemical  reactions  which 
occur  between  ground-water  and  aquifer  material. 

The  discussion  of  ground-water  chemistry  will  be 
in  two  parts.  The  first  part  is  a  brief  and  general 
description  of  water  quality  in  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers.  It  is  meant  primarily  to 
describe  the  water  quality  of  the  two  aquifers  and 
show  how  this  quality  changes  from  place  to  place. 
The  second  part  is  a  more  detailed  discussion  of 


Table  9.— Representative  analyses  of  Aquia  water 
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the  chemical  evolution  of  Aquia  and  Piney  Point- 
Nanjemoy  water.  For  this  discussion,  principles  of 
equilibrium  chemistry  are  utilized  to  identify  the 
major  chemical  processes  controlling  the  chemical 
character  of  the  water. 

Water  Quality 

Aquia  Aquifer 

Water  produced  from  the  Aquia  aquifer  generally 
has  good  chemical  quality.  Near  the  outcrop/subcrop 
area  of  the  Aquia  the  water  is  a  calcium  and  mag¬ 
nesium  bicarbonate  type  water  with  relatively  low 
(6. 5-7. 5)  pH.  Downgradient  from  the  outcrop/sub¬ 
crop  area  Aquia  water  changes  to  a  sodium  bicar¬ 
bonate  type  water  with  relatively  high  (7. 5-8. 5)  pH. 
Iron  concentrations  tend  to  be  high  (0.1-0. 5  mg/L) 
near  the  outcrop  area,  but  downgradient  are  very 


low  ((  0.1  mg/L).  Water  in  the  Aquia  aquifer  has 
generally  low  (8-12  mg/L)  concentrations  of  sulfate, 
which  do  not  change  significantly  along  the  flow- 
path.  Aquia  water  normally  contains  about  2  to  3 
mg/L  of  chloride,  which  is  close  to  the  concentration 
of  chloride  in  precipitation  in  Southern  Maryland 
(Junge  and  Werby,  1958).  Aquia  water  is  naturally 
fluoridated  and  typically  contains  0.1  to  0.3  mg/L  of 
fluoride.  In  addition  to  these  dissolved  chemical 
species,  Aquia  water  commonly  contains  dissolved 
silica,  manganese,  potassium,  and  nitrate  in  minor 
amounts.  Representative  analyses  of  Aquia  water 
are  shown  in  Table  9. 


Piney  Point-Nanjemoy  Aquifer 

In  many  respects,  ground  water  produced  from 
the  Piney  Point-Nanjemoy  aquifer  is  very  similar  in 
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chemical  quality  to  water  produced  from  the  Aquia 
aquifer.  This  is  primarily  because  the  lithologies  of 
the  Piney  Point-Nanjemoy  and  Aquia  aquifers  are 
very  similar.  In  upgradient  areas,  Piney  Point- 
Nanjemoy  water  is  a  calcium  magnesium  bicarbon¬ 
ate  type  water  with  relatively  low  (6. 0-7.0)  pH. 
Downgradient,  the  water  changes  to  a  sodium  bi¬ 
carbonate  type  water  with  relatively  high  (7. 0-8. 2) 
pH.  Concentrations  of  other  dissolved  ions  such  as 
iron  (0-0.5  mg/L),  sulfate  (10  mg/L),  chloride  (2-3 
mg/L),  and  fluoride  (0. 1-0.5  mg/L),  also  tend  to  be 


very  similar  to  Aquia  water.  Two  dissolved  species 
which  are  notably  different  from  Aquia  water  are 
potassium  and  silica.  Piney  Point-Nanjemoy  water 
commonly  has  about  12-15  mg/L  potassium,  which 
is  about  double  the  average  potassium  concentra¬ 
tions  in  Aquia  water.  Similarly,  silica  concentra¬ 
tions  in  Piney  Point-Nanjemoy  water  average  be¬ 
tween  50  and  60  mg/L,  which  is  about  4  or  5  times 
the  average  composition  of  silica  in  Aquia  water. 
Representative  analyses  of  Piney  Point-Nanjemoy 
water  are  shown  in  table  10. 


Table  10.— Representative  analyses  of  Piney  Point-Nanjemoy  water 
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Chemical  Evolution  of  Ground  Water 
Background 

In  the  past  30  years,  the  application  of  equilibrium 
chemistry  to  surface-  and  ground-water  systems  has 
provided  insight  into  many  types  of  geologic  prob¬ 
lems.  Foster  (1950),  Garrels  (1967),  and  Garrels  and 
MacKenzie  (1967)  have  applied  the  principles  of 
equilibrium  chemistry  to  the  composition  of  ground 
water  in  different  lithologic  terrains.  Jones  (1966), 
and  Eugster  and  Jones  (1979)  have  shown  how  the 
composition  of  brine  lakes  evolve  in  a  systematic 


fashion  which  depends  on  the  chemistry  of  the  en¬ 
closed  basin.  Similarly,  Cleaves  and  others  (1970) 
applied  the  principles  of  equilibrium  chemistry  to 
the  weathering  of  metamorphic  silicate  rocks  in 
Maryland  and  pointed  out  the  geomorphic  implica¬ 
tions  of  such  weathering.  Clearly  the  techniques  of 
ground-water  geochemistry  as  they  have  developed 
are  proving  to  be  excellent  tools  in  geologic  investi¬ 
gations.  In  this  report,  these  principles  are  utilized 
to  provide  a  framework  within  which  the  chemical 
composition  of  ground  water  from  the  Aquia  and 
Piney  Point-Nanjemoy  aquifers  in  Southern  Mary¬ 
land  may  be  understood. 
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Procedures 


Aquia  Aquifer 


All  chemical  analyses  of  ground  water  used  in  sub¬ 
sequent  calculations  have  been  taken  from  the 
WATSTORE  data  base  maintained  by  the  U.S. 
Geological  Survey.  Only  analyses  that  had  less  than 
a  5-percent  charge-balance  error  were  used.  Activi¬ 
ties  and  molalities  of  chemical  species  were  obtained 
by  the  use  of  the  computer  program  WATEQF 
(Plummer  and  others,  1976).  Additionally,  WATEQF 
was  used  to  calculate  the  calcite  saturation  indices 
of  Aquia  water.  Data  on  the  distribution  of  calcite 
cementation  in  the  Aquia  Formation  were  obtained 
from  drilling  logs  published  by  Bennion  and  Brook- 
hart  (1949),  Cook  and  others  (1952),  Martin  and  Fer¬ 
guson  (1953),  and  Otton  (1955). 

The  major  ions  in  Aquia  water  are  calcium,  mag¬ 
nesium,  sodium,  and  bicarbonate,  which  collectively 
account  for  about  95  percent  of  all  dissolved  constit¬ 
uents.  The  reaction  models  developed  in  this  paper 
are  designed  to  explain  the  distributions  of  these 
four  ionic  species. 

The  general  procedures  used  to  develop  the  chemi¬ 
cal  reaction  models  described  in  this  report  are: 

1.  Document  the  distribution  of  dissolved  cal¬ 
cium  (Ca2+),  magnesium  (Mg2+),  sodium  (Na+), 
and  bicarbonate  (HC03 )  and  demonstrate  how 
these  distributions  relate  to  aquifer  mineral¬ 
ogy  and  flow  patterns. 

2.  Develop  a  set  of  reaction  models  to  describe 
the  observed  changes  in  water  chemistry  along 
the  flowpath. 

3.  Test  the  reaction  models  by  comparing  ob¬ 
served  molar  ratios  of  dissolved  ions  with 
molar  ratios  predicted  by  the  reaction  models. 

The  procedure  for  evaluating  the  origin  of  calcite 
cementation  in  the  Aquia  Formation  will  be  to: 

1.  Perform  mass  balance  calculations  with  the 
verified  chemical  reaction  models. 

2.  Compare  observed  zones  of  calcite  cementation 
to  zones  predicted  by  the  mass  balance  cal¬ 
culations. 

3.  Show  that  the  zones  of  calcite  cementation  and 
zones  that  lack  cementation  as  predicted  by 
mass  balance  calculations  are  consistent  with 
the  calcite  saturation  state  of  Aquia  water. 

Analyses  of  Aquia  glauconite  for  cation  exchange 
properties  were  performed  at  the  U.S.  Geological 
Survey  Denver  Central  Laboratory  using  standard 
techniques.  Glauconite  was  separated  from  other 
aquifer  materials  for  these  analyses  using  a  mag¬ 
netic  separator. 


Mineralogy 

The  mineralogy  of  the  Aquia  aquifer  is  dominated 
by  quartz  sand,  glauconite,  and  shell  debris  which 
together  make  up  about  98  percent  of  the  aquifer 
material.  Quartz  is  the  most  abundant  mineral 
(50-65  percent).  Glauconite  is  the  second  most  abun¬ 
dant  mineral  (20-45  percent).  Broken  shell  debris 
commonly  makes  up  between  1  and  10  percent  of  the 
aquifer  material.  In  addition  to  these  major  constit¬ 
uents,  the  Aquia  contains  traces  of  magnetite, 
garnet,  kaolinite,  ilmenite,  hornblende,  tourmaline, 
epidote,  augite,  andalusite,  sillimenite,  kyanite, 
staurolite,  potassium  feldspar,  plagioclase,  pyrite, 
and  lignite. 

All  of  the  minerals  in  the  Aquia  react  to  some  de¬ 
gree  with  ground  water.  The  complexity  generated 
by  this  mineralogy  makes  it  impossible  to  construct 
a  unique  reaction  model  describing  the  water  chem¬ 
istry  (Plummer  and  others,  1983).  However,  the  ef¬ 
fects  of  slow-reacting  minerals,  such  as  quartz,  are 
small  compared  to  the  effects  of  fast-reacting  miner¬ 
als,  such  as  carbonates  (Garrels,  1975).  This  will  be 
particularly  true  if  the  faster-reacting  minerals  are 
much  more  abundant  than  the  slower-reacting  min¬ 
erals.  The  approach  taken  in  this  paper  is  to  assume 
that  fast-reacting  minerals  have  a  much  greater  im¬ 
pact  on  water  chemistry  than  minerals  that  react 
relatively  slowly.  If  this  assumption  results  in  an  ac¬ 
ceptable  approximation,  then  the  reaction  models 
which  are  constructed  will  be  non-unique  but  never¬ 
theless  useful. 

Under  ambient  temperature,  pressure,  and  pH 
conditions  in  the  Aquia  aquifer,  the  oxide  and  sili¬ 
cate  trace  minerals  can  be  expected  to  react  rela¬ 
tively  slowly  with  water.  For  the  purpose  of  this 
paper,  they  will  be  considered  nonreactive.  Similar¬ 
ly,  quartz  can  be  expected  to  react  relatively  slowly 
and  will  also  be  considered  nonreactive. 

Glauconite  is  an  iron-bearing,  hydrous  aluminum 
silicate  of  the  illite  family  and  is  an  efficient  cation 
exchange  medium.  Cation  exchange  reactions  are 
relatively  rapid  (Foster,  1950),  and  will  probably 
have  a  significant  impact  on  the  chemistry  of  Aquia 
water. 

The  shell  material  in  the  Aquia  consists  of  calcium 
carbonate  present  both  as  calcite  and  aragonite.  Ad¬ 
ditionally,  shell  material  commonly  contains  up  to 
18  mole  percent  magnesium  calcite  (Wollast  and 
others,  1971).  Under  most  conditions  encountered  in 
fresh-water  aquifers,  chemical  reactions  involving 
carbonates  are  relatively  rapid  (Reardon,  1981).  For 
this  reason,  the  carbonate  shell  material  present  in 
the  Aquia  can  be  expected  to  have  a  major  influence 
on  the  water  chemistry. 
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Carbon  dioxide  (C02)  contributes  to  the  chemical 
evolution  of  Aquia  water.  Carbon  dioxide  is  present 
in  soil  gas  and  can  react  with  percolating  water  in 
the  unsaturated  zone  (Pearson  and  Friedman,  1970). 
Carbon  dioxide  can  also  be  generated  in  the  aquifer 
by  slow  oxidation  of  lignite  (Foster,  1950).  In  view 
of  these  sources,  it  is  probably  reasonable  to  assume 
that  the  Aquia  aquifer  is  an  open  system  with  re¬ 
spect  to  carbon  dioxide. 

Ion  Distributions 

A  convenient  method  of  showing  the  distribution 
of  dissolved  ionic  species  in  an  aquifer  is  to  graph 
ion  concentrations  as  a  function  of  distance  along 
flowpaths.  This  method  has  the  advantage  of  show¬ 
ing  directly  how  water  chemistry  evolves  as  it 
moves  down  the  hydrologic  gradient.  In  the  Aquia 
aquifer,  the  distance  that  water  has  moved  along 
flowpaths  can  be  estimated  from  the  prepumping 
potentiometric  surface  (pi.  3).  As  a  matter  of  conven¬ 
tion  the  northern  limit  of  the  outcrop  area  can  be 
chosen  as  zero  distance. 

The  concentrations  of  Ca2++Mg2+,  Na+,  and 
HC03  plotted  as  a  function  of  distance  along  the 
flowpath  are  shown  in  figure  42.  The  curves  drawn 
through  the  data  points  are  to  illustrate  major 
trends  and  do  not  represent  a  rigorous  statistical  fit. 
The  data  of  figure  42  suggest  that  the  Aquia  aquifer 
can  be  divided  into  three  regions  based  on  the 
changes  in  water  chemistry: 

Region  /—This  region  is  about  26  miles  wide.  It  is 
parallel  to  and  includes  the  outcrop/subcrop  area 
of  the  aquifer.  Ca2++Mg2+  and  HC03  concentra¬ 
tions  increase  rapidly  and  level  off.  Concentra¬ 
tions  of  Na+  remain  relatively  low. 

Region  II— This  region  is  immediately  downgrad- 
ient  of  Region  I  and  is  about  22  miles  wide.  The 
water  chemistry  is  characterized  by  relatively  con¬ 
stant  HC03  concentrations,  rising  Na+  concentra¬ 
tions,  and  declining  Ca2++Mg2+  concentrations. 

Region  III— This  region  is  downgradient  from 
Region  II  and  continues  to  the  downdip  extent  of 
the  aquifer.  The  water  chemistry  is  characterized 
by  low  Ca2++Mg2+  concentrations  and  rapidly  in¬ 
creasing  Na+  and  HC03  concentrations. 

In  addition  to  the  changes  in  water  chemistry 
shown  in  figure  42,  the  pH  of  Aquia  water  consis¬ 
tently  increases  downgradient.  In  Region  I,  pH 
ranges  from  about  6.0  to  7.5.  In  Region  II,  pH 
ranges  from  about  7.5  to  8.2,  and  in  Region  III,  pH 
ranges  from  8.0  to  8.9.  Figure  43  shows  Regions  I, 
II,  and  III  in  plan  view,  and  representative  water 
analyses  for  each  region. 


EXPLANATION 


Na  + 

HCO  3 


Figure  42.— Concentrations  of  dissolved  sodium, 
calcium  and  magnesium,  and 
bicarbonate  versus  distance  along 
the  flowpath  for  Aquia  water. 

The  presence  of  different  water-chemistry  pat¬ 
terns  in  different  parts  of  the  Aquia  aquifer  implies 
that  the  chemical  processes  also  differ.  If  this  is  the 
case,  a  single  chemical  model  will  not  be  sufficient  to 
explain  the  changes  in  water  chemistry.  The  ap¬ 
proach  taken  will  be  to  build  separate  (but  related) 
chemical  models  for  each  region  of  the  Aquia  aquifer. 

Development  of  Working  Equations 

The  lithology  of  the  Aquia  aquifer  suggests  that 
shell  material,  carbon  dioxide,  and  glauconite  are 
the  aquifer  materials  that  most  affect  the  chemistry 
of  Aquia  water.  In  order  to  quantify  these  effects, 
chemical  equations  that  describe  the  reaction  of 
each  phase  with  water  must  be  developed. 

Shell  material  is  commonly  composed  of  at  least 
three  phases— calcite,  aragonite,  and  magnesium 
calcite  (Wollast  and  others,  1971).  The  reaction  of 
calcite  and  aragonite  with  water  can  be  represented 
by  the  equation 

CaC031sl+ H20(L)  ^  Ca^+HCC^+OH-,,  (3) 

Calcite  containing  a  significant  percentage  of  mag¬ 
nesium  can  be  treated  similarly: 

(Cal_IMgI)C03l,1+H20lu^ 

(l-x^aa+HCO^+OH,;,,  (4) 

The  reaction  of  C02  with  water  can  be  written 
C02Is1+H20(L1  ^H^+HCO*.,,  (5) 
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The  cation  exchange  reaction  with  glauconite  acting 
as  the  exchange  medium  can  be  written 

Xft+Na,  •  Glau,.d|  ^  2  NaJ^+X  •  Glau„dl  (6) 

where  X  can  be  Ca2+  or  Mg2+. 

These  equations  are  based  upon  idealized  mineral 
compositions  and,  therefore,  approximate  the  natu¬ 
ral  system. 


Development  of  Reaction  Models 

In  the  simplified  system  that  was  considered  in 
developing  the  working  equations,  the  source  of 
Ca2+  and  Mg2+  ions  is  dissolution  of  calcite,  arago¬ 
nite,  and  magnesium  calcite  (eqs.  3  and  4).  The  only 
important  source  of  Na+  is  the  glauconite  exchange 
reaction  (eq.  6).  HC03  ions  are  produced  by  dissolu¬ 
tion  of  calcite,  aragonite,  and  magnesium  calcite 
(eqs.  3  and  4)  and  by  reaction  of  C02  gas  with  water 
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Figure  43.— Regions  I,  II,  and  III  of  the  Aquia  aquifer  and  representative  water  analyses  from  each  region. 
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(eq.  5).  Various  combinations  of  these  working  equa¬ 
tions  can  be  utilized  to  account  for  the  water  chemis¬ 
try  changes  in  Region  I,  II,  and  III. 

In  Region  I,  sources  of  Ca2+,  Mg2+,  and  HCO~ions 
are  required.  Because  Na+  remains  essentially  con¬ 
stant  in  this  region,  a  Na+  source  is  not  required. 
Therefore,  the  chemistry  of  Region  I  can  be  simu¬ 


lated  with  the  working  equations: 

CaC03(g)+H20(L)  Cat^)-t-HC03(a(l)+0H(aq)  (3) 

<Ca1_IMgI)C03M+H20,LI  - 

<1 ' -’t)Caa+xMga+HC03';iql+OHM  (4) 

COw+H20(L)^H(tq)+HCO-.„  (5) 

Adding  these  equations  algebraically  and  balancing 
the  coefficients  yields: 

(Ca1_IMgI)C031Ij+C02(g|+H20(L|  - 
(1  -x)Ca^+xMg»+2HCO;>q)  |7) 


Given  the  simplifying  assumptions  made  in  devel¬ 
oping  the  working  equations,  equation  7  is  a  stoichi¬ 
ometric  model  for  the  ground-water  chemistry  of 
Region  I  and  will  be  referred  to  as  Model  I.  The  coef¬ 
ficients  of  Model  I  assume  that  the  pH  of  the  sys¬ 
tem  is  constant.  However,  because  H+  and  OH  -  con¬ 
centrations  are  low  ( ~  10-8  M  and  ~  10-6  M,  re¬ 
spectively)  compared  to  the  concentrations  of  the 
major  ions  (  ~  10-3  M  to  10-4),  the  stoichiometry  is 
valid  for  the  range  of  pH  found  in  Region  I. 

In  Region  II,  HC03  concentrations  remain  ap¬ 
proximately  constant.  Concentrations  of 
Ca2++Mg2+  decrease  and  Na+  concentrations  in¬ 
crease.  Because  no  source  of  HC03  is  required,  equa¬ 
tions  3,  4,  and  5  are  not  required.  In  this  region  only 
equation  6,  which  is  a  sink  for  Ca2++Mg2+  and  a 
source  of  Na+,  applies. 

X(2J+Na2  •  Glau(ad)^X  •  Glau(ad)+2Na+q)  (6) 
where 

X=Ca2+  or  Mg2+. 

This  equation  describes  the  major  changes  in  water 
chemistry  that  occur  in  Region  II  and  will  hereafter 
be  referred  to  as  Model  II. 

In  Region  III,  Ca2++Mg2+  concentrations  remain 
constant  whereas  Na+  and  HC03  concentrations  in¬ 
crease  rapidly.  Because  a  source  of  HC03  is  needed 
in  this  region,  equations  3,  4,  and  5  apply.  Because  a 
source  of  Na+  and  a  sink  for  Ca2+  and  Mg2+  ions  is 
needed,  equation  6  also  applies.  The  working  equa¬ 
tions  for  this  system  will  be 

CaC03(S)+H20,L)  ^  Ca^+HC03(aq)+0H(aq)  (3) 
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Figure  44.— Calcium  and  magnesium  concentra¬ 
tion  versus  bicarbonate  concentration 
for  water  analyses  in  Region  I  of  the 
Aquia  aquifer. 

(Cai-xMgx)C03(s)+H20(L)^ 

(1  -x)Ca2;+xMg2a;+HC03+0H(;q)  (4) 

GG2(g)^”H20(Lij  T—  H(aq)  +  HC03(aq)  ^ 

X(2at)  +  Na2  *  GlaU(ad)  ^  2  Na(tq)  +  X  *  GlaU<ad)  (6) 

where 

X=Ca2+  or  Mg2+. 

Adding  these  equations  algebraically  and  balancing 
the  coefficient  yields 

(Ca1_„Mg2)C03l,l+C02(g)+H20(LI+Na2  •  Glaulld)  = 

2Na,:q)+2HCO-.„+Ca1_>MgI  •  Glau(1(i)  (8) 

This  equation  describes  the  water-chemistry  changes 
in  Region  III  and  will  hereafter  be  referred  to  as 
Model  III.  This  equation  is  similar  to  one  derived  by 
Thorstenson  and  others  (1979). 

Model  Predictions 

Models  I,  II,  and  III  must  be  tested  to  determine 
if  they  accurately  predict  water-chemistry  changes 
in  the  regions  for  which  they  were  designed.  One 
method  of  testing  these  models  is  to  graphically 
compare  observed  molar  ratios  of  dissolved  species 
with  those  predicted  by  the  models.  Good  matches 
between  predicted  and  observed  curves  will  support 
the  applicability  of  the  models. 
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Model  I  states: 

(CBi-xMgJCOg^  +  CO^  +  HgO^,  T— 

( 1  “  X>Ca  fat)  +  xMgfaJ  +  2HCO  “aq). 

Inspection  of  this  equation  shows  that  as  carbonate 
shell  material  dissolves,  the  molar  proportion  of 
Ca2++Mg2+  to  HC03  should  be  linear  with  a  slope  of 
+  2.0.  Because  it  is  difficult  to  quantify  the  relative 
proportions  of  Ca2+  and  Mg2+  in  shell  material,  they 
have  been  considered  together.  Figure  44  shows 
Ca2++Mg2+  plotted  as  a  function  of  HC03  for  analy¬ 
ses  of  water  in  Region  I.  The  linear  least  squares 
best  fit  for  this  plot  has  a  slope  of  1.87,  which  is 
close  to  the  predicted  slope  of  2.0.  The  r-square  cor¬ 
relation  for  this  plot  is  0.67,  indicating  a  good  fit  be¬ 
tween  the  model  and  the  data.  The  statistical  signifi¬ 
cance  of  this  correlation  relative  to  sample  size  can 
be  tested  using  the  Fisher  Z  transformation  (Miller 
and  Freund,  1977,  p.  325).  Use  of  this  test  statistic 
shows  that  the  correlation  of  Ca2++Mg2+  to  HC03 
is  significant  at  the  95-percent  confidence  level. 

Model  II  states: 

xfaJ+Na2  *  Glau(ad)^2Na+q)+x  •  Glau(ad) 
X=Ca2+  or  Mg2+. 

This  equation  predicts  that  Na+  ions  will  enter  solu¬ 
tion  at  twice  the  rate  that  Ca2++Mg2+  ions  are 
removed  from  solution.  Therefore,  the  molar  rela¬ 
tionship  between  Ca2++Mg2+  and  Na+  in  Region  II 
should  be  a  straight  line  with  slope  —2.0.  Figure  45 
shows  Ca2++Mg2+  plotted  as  a  function  of  Na+  for 
water  analyses  taken  in  Region  II.  The  linear  least 
squares  best  fit  for  this  plot  has  a  slope  of  —1.75, 
which  is  close  to  the  predicted  slope  of  —2.0.  Also, 
the  r-square  correlation  for  this  plot  is  0.86  indicat¬ 
ing  a  good  fit  between  the  model  and  the  data.  This 
correlation  is  significant  at  the  95-percent  confi¬ 
dence  level. 

Model  III  states: 

(Caj_x  Mgx)C03(s)+C02(g)+H20(L)+Na2  •  Glau(ad) 

^  2Na2J+2HCO-aq)+Cai_x,  Mgx)  •  Glau(ad, 

This  equation  predicts  that  the  molar  proportions  of 
Na+  and  HC03  will  be  related  linearly  with  a  slope 
of  +1.0.  Figure  46  shows  the  plot  of  Na+  as  a  func¬ 
tion  of  HC03  for  all  available  water  analyses  in 
Region  III.  The  least  squares  best  fit  for  this  data  is 
a  line  with  slope  0.99,  which  is  very  close  to  the  pre¬ 
dicted  slope  of  1.0.  The  linear  r-square  for  this  plot  is 
0.98,  which  shows  excellent  correlation  between  the 
model  and  the  data.  This  correlation  is  significant  at 
the  95-percent  confidence  level. 
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Figure  45.— Calcium  and  magnesium  concentra¬ 
tion  versus  sodium  concentration  for 
water  analyses  in  Region  II  of  the 
Aquia  aquifer. 
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Figure  46.— Sodium  concentration  versus  bicar¬ 
bonate  concentration  for  water  analy¬ 
ses  in  Region  III  of  the  Aquia  aquifer. 

The  generally  close  matches  between  predicted 
water  chemistry  and  measured  water  chemistry  sup¬ 
ports  the  validity  of  Models  I,  II,  and  III.  These 
models  were  constructed  by  assuming  relatively 
simple  aquifer  mineralogy  and  by  assuming  ideal¬ 
ized  mineral  compositions,  and  therefore  these  are 
not  unique  models.  However,  for  the  system  speci¬ 
fied,  they  work  well. 
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EXPLANATION 

▲ 

CONTROL  POINT  FOR  CATION 
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• 

CONTROL  POINT  FOR  CON¬ 
CENTRATION  OF  EXCHANGE¬ 
ABLE  SODIUM  ON  AQUIA 
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IN  MILLIMOLES  PER  GRAM 

Figure  47.— Cation  exchange  capacity  and  concentration  of  exchangeable  sodium  of  glauconite  samples  from 
the  Aquia  aquifer  plotted  versus  distance  along  the  flowpath. 
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Discussion 

It  has  been  shown  that  concentration  changes  of 
Ca2+,  Mg2+,  Na+,  and  HC03  in  Aquia  water  can  be 
described  by  three  models.  Model  I  applies  to  the 
region  including  and  closest  to  the  outcrop  area. 
Model  II  applies  to  the  region  immediately  down- 
gradient  from  this.  Model  III  applies  to  the  most 
southeastern  area  of  the  Aquia  aquifer.  These  re¬ 
gions  are  shown  in  figure  43.  Although  these  models 
can  predict  water-chemistry  changes  along  the  flow- 
path,  they  do  not  explain  why  certain  reactions  oc¬ 
cur  in  some  regions  of  the  aquifer,  but  not  in  other 
regions.  Two  questions  in  particular  should  be  ad¬ 
dressed: 

1.  Why  is  cation  exchange  an  important  reaction 
in  Regions  II  and  III,  but  not  in  Region  I? 

2.  Why  is  shell  material  dissolution  an  important 
reaction  in  Regions  I  and  III,  but  not  in  Region 
II? 

The  first  question  was  addressed  by  determining 
the  cation  exchange  properties  of  some  Aquia 
glauconite  samples.  Six  samples  of  Aquia  aquifer 
material  were  obtained.  Five  of  these  samples  were 
from  drill  cuttings  and  one  was  from  an  Aquia  out¬ 
crop.  The  glauconite  was  separated  from  the  sam¬ 
ples  and  analyzed  for  cation  exchange  capacity  and 
concentrations  of  exchangeable  cations.  The  cation 
exchange  capacity  and  concentration  of  exchange¬ 
able  sodium  plotted  as  a  function  of  distance  along 
the  flowpath  are  shown  in  figure  47.  Also  shown  are 
the  locations  of  the  six  sample  sites.  The  cation  ex¬ 
change  capacity  of  these  samples  is  relatively  con¬ 
stant.  The  concentrations  of  exchangeable  sodium, 
however,  are  low  near  the  outcrop  area  but  increase 
substantially  at  about  30  miles  downgradient.  A 
comparison  of  figure  47  with  figure  42  shows  that 
the  increase  of  exchangeable  sodium  on  the  glauco¬ 
nite  coincides  with  the  increase  of  dissolved  Na+  in 
Aquia  water.  It  can  be  concluded,  therefore,  that  ca¬ 
tion  exchange  is  not  an  important  reaction  in  Region 
I  because  the  glauconite  has  low  concentrations  of 
exchangeable  sodium.  Cation  exchange  occurs  in 
Regions  II  and  III  because  the  glauconite  has 
higher  concentrations  of  exchangeable  sodium. 

The  lower  exchangeable  sodium  content  of  up- 
gradient  glauconite,  as  opposed  to  that  of  down- 
gradient  glauconite,  has  probably  developed  over 
the  history  of  the  aquifer.  In  upgradient  areas,  Ca2+ 
and  Mg2+  produced  by  shell  material  dissolution 
would  quickly  react  with  glauconite  and  thus 
deplete  exchangeable  sodium.  Water  from  upgrad¬ 
ient  areas  (now  low  in  Ca2+  and  Mg2+)  would  subse¬ 
quently  move  downgradient.  Glauconite  in  down- 
gradient  areas  would  be  in  contact  with  water  con¬ 


taining  little  Ca2+  and  Mg2+  and  exchangeable 
sodium  would  not  be  depleted  as  quickly. 

The  second  question  can  be  addressed  by  evaluat¬ 
ing  the  equilibrium  of  calcite  in  the  Aquia  aquifer. 
The  equilibrium  of  a  mineral-water  reaction  is  usu¬ 
ally  defined  in  terms  of  the  Saturation  Index  (SI).  SI 
is  defined  as 


IAP=ion  activity  product  of  the  mineral-water 
reaction. 

K=the  equilibrium  constant  for  the  mineral- 
water  reaction. 

If  the  SI  for  calcite  is  zero  or  close  to  zero  for  a 
particular  water  analysis,  then  that  water  is  in 
equilibrium  with  calcite.  If  SI  is  negative,  the  water 
is  subsaturated.  Langmuir  (1971)  suggests  that  in¬ 
dexes  within  ±0.1  of  zero  should  be  considered 
saturated.  Therefore,  if  SI  is  zero  ±0.1  or  positive, 
dissolution  of  calcite  should  not  occur.  If  SI  is  more 
negative  than  —0.1,  dissolution  of  calcite  should 
occur. 

The  calcite  Si’s  as  calculated  by  WATEQF  for 
available  analyses  of  Aquia  water  in  Southern  Mary¬ 
land  are  shown  in  figure  48.  In  Region  I,  the  Aquia 
aquifer  is  recharged  by  meteroic  water  that  is  rela¬ 
tively  unmineralized  and  is  subsaturated  with  re¬ 
spect  to  calcite.  When  this  water  enters  the  aquifer 
it  will  rapidly  dissolve  calcite  until  equilibrium  is 
reached  (Reardon,  1981).  Figure  48  shows  that  water 
analyses  in  Region  I  are  either  close  to  saturation  or 
supersaturated  (average  SI  =  ±0.07)  indicating  that 
calcite  dissolution  has  occurred  and  equilibrium  is 
being  approached.  In  Region  II  where  water  has 
moved  downgradient  from  Region  I,  the  analyses 
show  that  the  water  remains  close  (average 
SI  =  —0.07)  to  calcite  saturation.  Therefore,  signifi¬ 
cant  calcite  dissolution  in  Region  II  is  not  to  be  ex¬ 
pected.  In  Region  III,  practically  all  analyses  are 
subsaturated  (average  SI  =—0.31)  and  calcite  disso¬ 
lution  should  occur.  This  subsaturation  in  Region 
III  is  probably  due  to  the  loss  of  Ca2+  ions  to  the  ca¬ 
tion  exchange  reaction.  This  lowers  the  IAP  for  the 
calcite  equilibrium  and  results  in  subsaturation. 

Calcite  Cementation 

The  Aquia  aquifer  in  places  is  characterized  by 
considerable  thickness  of  sandstone  ledges  that  are 
cemented  by  calcite.  These  ledges,  or  “hardbeds”  as 
they  are  often  called,  are  distinctive  because  their 
relative  hardness  makes  them  difficult  to  drill.  Indi¬ 
vidual  hardbeds  range  in  thickness  from  a  few 


84 


76°45' 


76°  30' 


76°  15' 


85 


76°45' 


76°  30'  76°  15' 


Figure  49.— Approximate  thickness  of  calcite  cementation  in  the  Aquia  aquifer. 


86 


inches  to  several  feet  and  are  often  associated  with 
lenses  of  abundant  shell  material.  There  may  be  any¬ 
where  from  one  to  several  dozen  separate  hardbeds 
in  a  given  vertical  section  of  the  Aquia.  The  spatial 
distribution  of  this  calcite  cementation  follows  a 
pattern.  Localities  in  and  near  the  outcrop  area  are 
characterized  by  abundant  calcite  cementation, 
whereas  areas  further  downgradient  exhibit  little 
cementation.  This  pattern  is  illustrated  in  figure  49, 
which  shows  thicknesses  of  cementation  as  recorded 
in  published  drilling  logs.  In  Prince  Georges  and 
Anne  Arundel  Counties,  which  are  near  the  outcrop 
area,  total  thicknesses  of  30  feet  or  greater  are  com¬ 
mon.  Further  downgradient  in  St.  Mary’s  County, 
well  logs  record  little  or  no  cementation.  This  ob¬ 
served  difference  is  not  related  to  a  change  in  avail¬ 
able  shell  material.  Shell  material  is  as  abundant  in 
areas  lacking  cementation  as  it  is  in  areas  with 
abundant  cementation. 

The  calcite  cementation  may  be  a  primary  sedi¬ 
mentary  feature  of  the  Aquia  as  has  been  proposed 
by  Hansen  (1974)  and  Wright  and  Huffman  (1979). 
It  is  also  possible  that  the  cementation  is  secondary 
in  origin  and  has  resulted  from  the  reaction  of 
ground  water  with  aquifer  material.  If  the  calcite 
cementation  is  secondary,  the  chemical  reaction 
models  developed  in  this  paper  should  be  consistent 
with  the  observed  spatial  distribution  of  the  cemen¬ 
tation.  The  procedure  now  will  be  to  determine  if  the 
chemical  reaction  models  are  consistent  with  the 
distribution  of  calcite  cementation. 

Mass  Balance 

In  order  to  test  if  the  reaction  models  are  consis¬ 
tent  with  the  distribution  of  calcite  cementation,  the 
mass  balance  technique  (Plummer  and  Back,  1980) 
was  used.  This  technique  has  the  advantage  of  quan¬ 
titatively  considering  solid  phase  products  of  reac¬ 
tions  and  can  take  into  account  the  production  or 
nonproduction  of  secondary  calcite. 

For  a  specified  system,  mass  balance  can  be  de¬ 
scribed  by  N  equations  of  the  form 

&  «,  0cj=  AMc  _liK  (9) 

where 

0  =  total  number  of  minerals  and  gases 
considered; 

N= minimum  number  of  constituents  needed  to 
define  the  composition  of  the  chosen 
minerals  and  gases; 

a  ^stoichiometric  coefficient  of  the  mineral  or 
gas  in  moles  per  Kg; 


d  .= stoichiometric  coefficent  of  the  cth 
constituent  in  the  jth  mineral  or  gas; 

A  Me = the  change  in  moles  per  Kg  H20  of  the  Cth 
constituent  in  the  aqueous  phase  along  the 
reaction  path. 

When  A  Me  is  known  from  chemical  analyses  and 
N  =  <P ,  a  mass  balance  of  the  chosen  minerals  and 
gases  can  be  obtained  by  simultaneous  solution  of 
the  N  equations.  This  algebraic  problem  can  be  con¬ 
veniently  solved  by  computers.  One  such  computer 
program,  BALANCE  (Plummer,  written  commun., 
1980),  was  used  to  perform  the  calculations  in  this 
paper. 

Mass  balance  calculations  were  performed  for  the 
water  chemistries  of  Regions  I,  II,  and  III  using  the 
reactions  of  Models  I,  II,  and  III.  In  Region  I,  the 
number  of  minerals  and  gases  ( 0 )  considered  by 
Model  I  was  equal  to  the  number  of  active  dissolved 
constituents  (N)  so  that  a  mass  balance  could  be 
found  directly.  For  Regions  II  and  III,  however,  it 
was  necessary  to  transform  the  equations  of  Models 
II  and  III  so  that  the  requirement  0  =N  was  satis¬ 
fied.  The  cation  exchange  reaction 

X,24,+Na2  •  Glau(adl^2Na,+q|+X  •  Glauladl  (6) 

X2+=Ca2+  or  Mg2+ 

was  written  as  two  reactions 

(Ca?++Na2  •  Glau(ad}  ^2Na(aql+Ca  •  Glaulld)  (101 

Mgf.;+Na2  •  Glau(Idl  =  2Na+q(+Mg  •  Glauladl(ll) 

Also,  the  calcite  and  Mg-calcite  dissolution-precipi¬ 
tation  reactions 

CaC0„+H20)LI=  Ca2;+HC03„+0HM  (3) 

and 

(Cai_xMgx)C03(s)  +  H20(L)^- 

(l-x)Ca,24l+xMg24+HC03M+0H-aql  (4| 

were  combined  with  the  carbon  dioxide  dissolution 
reaction 

+  ^  H^  +  HCOa^q)  (5) 

to  yield 

(Cai_xMgx)C03(s)+C02(g)+H20(L)  t— 

(1— x)Caa+xMg^q)+2HCO-aq)  (7) 

and 

CaC03(.,+C02(g)+H20(L)  -  Ca»+2HCO^,  (12) 

This  approach  eliminated  the  direct  consideration  of 
C02  gas  as  a  phase  in  BALANCE. 
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Table  11.— Results  of  mass  balance  calculations 
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Wollast  and  others  (1971)  show  that  natural  car¬ 
bonate  material  commonly  contains  between  1  and 
18  percent  magnesium  carbonate.  On  this  basis  it 
was  judged  that  10  mole  percent  magnesium  car¬ 
bonate  would  be  a  reasonable  estimate  of  shell  mate¬ 
rial  composition  in  the  Aquia.  Under  this  assump¬ 
tion,  equation  10  becomes 

(Ca9Mgi)C03(sl+C02(gl+H20(L1^ 

•9Ca^q, + .  1  Mg(aq| + 2HC03laqr  (13) 

In  Region  II,  it  was  necessary  to  specify  equations  7 
and  13  in  addition  to  equations  10  and  11  so  that 
</>=  N.  Equation  13  was  used  to  represent  the  dis¬ 
solution  reaction  of  magnesium  calcite  in  Region 
III. 

The  results  of  the  mass  balance  calculations  are 
shown  in  table  11.  These  calculations  use  the  water 
analyses  of  wells  CH-Cgl,  SM-Ddl,  and  SM-Ff35 
(fig.  43)  as  the  water  compositions  of  Regions  I,  II, 
and  III,  respectively.  The  transformed  equations  of 
Models  I,  II,  and  III  are  shown  in  table  11  together 
with  the  chemical  processes  that  the  equations  rep¬ 
resent.  Table  11  is  designed  to  show  the  composi¬ 
tional  changes  of  water  as  it  moves  along  the  flow- 
path.  The  average  composition  of  rainfall  in  south¬ 
ern  Maryland  (Junge  and  Werby,  1958)  is  acted  on 
by  the  Model  I  reactions  to  produce  the  water  com¬ 
position  of  Region  I.  Region  I  water  is  then  acted  on 
by  the  Model  II  reactions  to  produce  the  composi¬ 
tion  of  Region  II  water.  Finally,  Region  II  water  is 
acted  on  by  the  reactions  of  Model  III  to  produce 
the  composition  of  Region  III  water. 

The  calculations  summarized  in  table  11  show 
that  dissolution  of  magnesium  calcite  and  simultan¬ 
eous  precipitation  of  calcite  is  required  to  obtain 
mass  balance  in  Region  I.  In  Region  II,  the  calcula¬ 
tions  indicate  that  cation  exchange  reactions  are  the 
most  important  chemical  processes.  Magnesium  cal¬ 
cite  dissolution  and  calcite  precipitation  account  for 
only  a  small  portion  of  the  mass  balance.  This  is  con¬ 
sistent  with  the  earlier  conclusion  that  Region  II 
water  chemistry  is  dominated  by  cation  exchange 
reactions.  In  Region  III,  cation  exchange,  dissolu¬ 
tion  of  magnesium  calcite,  and  dissolution  of  calcite 
is  required  for  mass  balance. 

A  comparison  of  figures  43  and  49  shows  that  in 
Region  I,  where  mass  balance  requires  precipitation 
of  calcite,  the  greatest  thicknesses  of  calcite  ce¬ 
mentation  are  observed.  Region  II,  which  does  not 
require  significant  precipitation  of  calcite  for  mass 
balance,  shows  much  less  cementation.  Region  III, 
where  dissolution  of  magnesium  calcite  and  calcite 
is  needed  for  mass  balance,  calcite  cementation  is 
generally  absent.  This  mass  balance  is  also  consis¬ 
tent  with  the  equilibrium  chemistry  of  the  system. 


In  Region  I,  where  calcite  cementation  precipitation 
is  needed  for  mass  balance,  the  water  is  generally 
either  near  saturation  or  supersaturated  with  re¬ 
spect  to  calcite  (fig.  48).  Similarly  in  Region  III, 
where  calcite  dissolution  is  required,  the  water  is 
subsaturated  with  respect  to  calcite. 

While  the  preceding  analysis  presents  an  internal¬ 
ly  consistent  picture,  it  should  be  reemphasized  that 
the  chemical  models  used  are  not  unique.  Therefore, 
the  mass  balance  obtained  also  is  not  unique.  Spe¬ 
cifically,  the  exact  composition  of  magnesium  cal¬ 
cite  in  the  Aquia  is  not  known  and,  therefore,  this 
phase  cannot  be  evaluated  exactly.  However,  addi¬ 
tional  mass  balance  calculations  that  specified  5  and 
15  mole  percent  magnesium  calcite  resulted  in  exact¬ 
ly  the  same  patterns  of  dissolution-precipitation 
shown  in  table  11.  Only  the  relative  amounts  of 
magnesium  calcite  dissolving  and  calcite  precipitat¬ 
ing  were  affected.  Also,  because  aragonite  has  the 
same  stoichiometry  as  calcite,  its  relative  impor¬ 
tance  in  the  dissolution  and  precipitation  of  calcite 
cannot  be  evaluated  with  the  mass  balance  tech¬ 
nique.  Because  aragonite  is  less  stable  than  calcite, 
it  is  likely  that  dissolution  of  aragonite  and  dissolu¬ 
tion  of  magnesium-rich  calcite  contribute  to  the  up- 
gradient  supersaturation  of  water  with  respect  to 
calcite.  These  uncertainties,  however,  do  not  alter 
the  basic  picture  presented  in  this  paper.  The  chemi¬ 
cal  reaction  models  that  were  designed  to  predict 
changes  in  water  chemistry  along  the  flowpath  are 
also  consistent  with  the  observed  distribution  of  cal¬ 
cite  cementation.  This  strongly  suggests  that  the 
observed  calcite  cementation  in  the  Aquia  is  a  post 
depositional  feature  that  has  resulted  from  the  reac¬ 
tion  of  ground  water  with  aquifer  material. 

Piney  Point  Nanjemoy  Aquifer 
Background 

The  chemical  quality  of  Piney  Point-Nanjemoy 
water  is  very  similar  to  the  chemical  quality  of 
Aquia  water.  Both  are  calcium  magnesium  bicar¬ 
bonate  type  water  near  their  outcrop  areas,  and  both 
change  to  a  sodium  bicarbonate  type  water  down- 
gradient.  This  is  not  surprising  because  the  mineral¬ 
ogy  and  ground-water  flow  patterns  of  the  two  aqui¬ 
fers  are  very  similar.  Given  this,  it  seems  reasonable 
that  the  techniques  which  were  used  to  investigate 
Aquia  water  chemistry  should  be  applicable  to 
Piney  Point-Nanjemoy  water  chemistry.  For  the 
purposes  of  this  report,  the  procedures  used  to  eval¬ 
uate  the  chemical  evolution  of  Piney  Point-Nanje¬ 
moy  water  will  be  the  same  as  those  used  for  the 
Aquia  aquifer. 
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Mineralogy 

The  Piney  Point-Nanjemoy  aquifer  consists  of 
olive-green  to  greenish-gray  glauconitic  quartz 
sand.  Sands  near  the  top  of  the  unit  have  less  glau¬ 
conite  (approximately  5  to  10-percent  glauconite) 
than  sands  near  the  bottom  of  the  unit  (40  to 
60-percent  glauconite).  Shell  material  is  common  in 
the  unit  and  makes  up  1  to  5  percent  of  the  aquifer 
material.  Staurolite,  sillimanite,  garnet,  tourmaline, 
and  zircon  are  the  most  commonly  identified  heavy 
minerals  present  in  the  Piney  Point-Nanjemoy 
(Glaser,  1971). 

For  the  purposes  of  this  report,  quartz  is  consid¬ 
ered  to  be  a  non-reactive  phase.  It  is  also  assumed 
that  the  minor  minerals  have  negligible  impact  on 
water  chemistry.  Given  these  assumptions,  shell 
material,  glauconite,  and  carbon  dioxide  gas  are  the 
phases  which  will  have  the  most  effect  on  water 
chemistry. 

Ion  Distributions 

Figure  50  shows  the  molar  concentrations  of 
calcium  and  magnesium,  sodium,  and  bicarbonate 
plotted  as  functions  of  distance  along  the  flowpath. 
The  curves  drawn  through  the  data  are  to  illustrate 
major  trends  and  do  not  represent  statistical  signifi¬ 
cance.  Figure  50  is  strikingly  similar  to  figure  42 
which  shows  calcium  and  magnesium,  sodium,  and 
bicarbonate  concentrations  versus  distance  for 
Aquia  water.  This  similarity  implies  that  the 
chemical  processes  which  control  the  composition  of 
water  from  the  two  aquifers  are  similar.  Inspection 
of  figure  50  suggests  that  the  Piney  Point-Nan¬ 
jemoy  aquifer  can  be  divided  into  three  regions 
based  on  observed  water-chemistry  changes: 
Region  7— This  region  is  about  35  mi  wide  and  is 
parallel  to  and  includes  the  outcrop  area.  Calcium 
and  magnesium  and  bicarbonate  concentrations 
are  high  and  constant  throughout  the  region,  and 
sodium  concentrations  are  low  and  constant. 
Region  II— This  region  is  immediately  downgradi- 
ent  from  Region  I,  and  is  about  8  mi  wide.  It  is 
characterized  by  decreasing  calcium  and  magne¬ 
sium,  increasing  sodium,  and  constant  bicar¬ 
bonate  concentrations. 

Region  III— This  region  is  downgradient  from 
Region  II  and  is  about  10  mi  wide.  It  is  character¬ 
ized  by  rapidly  increasing  sodium  and  bicarbonate 
concentrations,  and  constant  low  concentrations 
of  calcium  and  magnesium. 

Figure  51  shows  Regions  I,  II,  and  III  of  the 
Piney  Point-Nanjemoy  aquifer  in  plan  view.  Also 
shown  are  representative  analyses  of  water  from 
each  region.  The  different  trends  in  each  region  im- 


Figure  50.— Concentrations  of  dissolved  sodium, 
calcium  and  magnesium,  and 
bicarbonate  versus  distance  along 
the  flowpath  for  Piney  Point- 
Nanjemoy  water. 

ply  that  the  chemical  processes  occurring  in  each 
region  are  also  different.  The  approach  taken  now 
will  be  to  build  separate  chemical  models  for  each 
region  of  the  Piney  Point-Nanjemoy  aquifer. 

Working  Equations  and  Reaction  Models 

The  bulk  mineralogy  of  the  Piney  Point-Nanjemoy 
aquifer  is  very  similar  to  that  of  the  Aquia  aquifer. 
Because  of  this  similarity,  the  development  of  work¬ 
ing  equations  which  describe  the  dominant  chemical 
reactions  is  identical  to  the  development  used  for 
the  Aquia  aquifer.  Furthermore,  the  patterns  of  ma¬ 
jor  ion  concentration  changes  in  the  Piney  Point- 
Nanjemoy  and  Aquia  aquifers  are  very  similar. 
Because  of  this,  the  reaction  models  which  describe 
the  chemistry  of  the  Piney  Point-Nanjemoy  will  be 
identical  to  those  which  describe  the  chemistry  of 
the  Aquia.  The  only  difference  is  that  Regions  I,  II, 
and  III  of  the  aquifers  are  not  the  same.  The  de¬ 
tailed  procedures  for  development  of  working  equa¬ 
tions  and  reaction  models  have  been  given  previ¬ 
ously  in  the  section  on  the  chemical  evolution  of 
Aquia  water  and  will  not  be  repeated  here.  Follow¬ 
ing  those  procedures,  Model  I  (eq.  7)  will  be  a  chemi¬ 
cal  model  for  Region  I  of  the  Piney  Point-Nanjemoy 
aquifer.  Model  II  (eq.  6)  will  be  a  chemical  model  for 
Region  II,  and  Model  III  (eq.  8)  will  be  a  chemical 
model  for  Region  III  of  the  Piney  Point-Nanjemoy 
aquifer. 
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Model  Predictions 


The  chemical  models  can  be  tested  by  comparing 
observed  molar  relationships  of  major  ions  deter¬ 
mined  by  linear  regression  to  molar  relationships 
predicted  by  the  models.  Good  matches  between 
predicted  and  measured  curves  will  support  the  ac¬ 
curacy  of  the  models. 

Model  I  states 


(CV.MgJ  co,(s1+co2W+h2oiL- 

(1— x)Caj;i+xMg|.q)+2HCOw  (7) 

Inspection  of  this  equation  shows  that  in  Region  I 
the  molar  relationship  of  calcium  and  magnesium  to 
bicarbonate  should  be  a  straight  line  with  slope 
+2.0.  Figure  52  shows  calcium  and  magnesium 
plotted  versus  bicarbonate  for  available  water  anal¬ 
yses  in  Region  I.  The  least  square  best  fit  for  this 
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Figure  51.— Regions  I,  II,  III  of  the  Piney  Point-Nanjemoy  aquifer  and  representative  water  analyses  from  each 
region. 
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plot  has  a  slope  of  1.95,  which  is  very  close  to  the 
predicted  slope.  The  r-square  correlation  for  this 
plot  is  0.75  indicating  a  good  fit  between  the  model 
and  the  data.  This  correlation  is  significant  at  the 
99-percent  confidence  level. 

Model  II  states 

XfaJ  +  Na  *  GlaU(ad,^  2Na(aq)  +  X  *  GlaU(ad) 

where 

X=Ca2+  or  Mg2+. 

This  equation  predicts  that  the  molar  relationship 
between  calcium  and  magnesium  and  sodium  should 
be  a  straight  line  with  slope  —2.0.  Figure  53  shows 
calcium  and  magnesium  plotted  versus  sodium  for 
water  analyses  taken  in  Region  II.  The  least  squares 
best  fit  for  this  plot  has  a  slope  of  —1.20  and  the 
r-square  correlation  of  the  plot  is  0.68.  This  correla¬ 
tion  is  significant  at  the  99-percent  confidence  level. 

It  is  significant  that  the  slope  of  the  regression 
line  (  —  1.20)  differs  substantially  from  the  slope  pre¬ 
dicted  by  Model  II  (  —  2.0).  This  implies  that  the 
assumptions  used  to  build  Model  II  were  not  suffi¬ 
ciently  accurate  to  describe  the  chemistry  of  Region 
II  within  an  acceptable  error  margin.  Model  II 
assumes  that  cation  exchange  is  the  dominant 
chemical  process  in  Region  II.  The  negative  slope  of 
the  calcium  and  magnesium  versus  sodium  plot  sug¬ 
gests  that  cation  exchange  is  indeed  occurring. 
However,  because  the  slope  of  the  regression  is  not 
closer  to  —2.0,  it  can  be  concluded  that  chemical 
processes  other  than  cation  exchange  are  modifying 
the  composition  of  Region  II  water.  It  is  possible 
that  these  additional  chemical  processes  include 
dissolution-precipitation  reactions  of  carbonate 
shell  material. 

Model  III  states 

(CVxMgx)  C03(s)+H20(L)+Na2  Glau^ 

2Na(tq)+2HC03(aq)  +  (Cai-x.Mgx)  *  GlaU(ad)- 

This  equation  predicts  that  sodium  and  bicar¬ 
bonate  will  increase  with  a  1:1  molar  ratio.  Figure  54 
shows  the  plot  of  sodium  versus  bicarbonate  for 
available  water  analyses  in  Region  III.  The  least 
squares  best  fit  for  this  data  is  a  line  with  slope  0.99, 
which  is  very  close  to  the  predicted  slope  of  1.0.  The 
linear  r-square  for  this  plot  is  0.69  indicating  a  good 
fit  between  the  model  and  the  data.  This  correlation 
is  significant  at  the  99-percent  confidence  level. 

Discussion 

It  has  been  shown  that  concentrations  of  dis¬ 
solved  calcium,  magnesium,  sodium,  and  bicar¬ 
bonate  in  the  Piney  Point-Nanjemoy  aquifer  change 


Ca  ♦  Mg 


Figure  52.— Calcium  and  magnesium  concentra¬ 
tion  versus  bicarbonate  concentration 
for  water  analyses  in  Region  I  of  the 
Piney  Point-Nanjemoy  aquifer. 


Figure  53.— Sodium  concentration  versus  calcium 
and  magnesium  concentration  for 
water  analyses  in  Region  II  of  the 
Piney  Point-Nanjemoy  aquifer. 

systematically  along  the  flowpath.  Based  on  these 
changes  in  chemistry,  the  aquifer  was  divided  into 
three  regions.  Because  water  chemistry  of  these  re¬ 
gions  is  similar  to  that  of  the  Aquia  aquifer,  the 
chemical  models  derived  for  the  Aquia  aquifer  were 
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applied  to  the  chemistry  of  the  Piney  Point-Nanje- 
moy  aquifer.  Good  matches  between  predicted  and 
measured  water  chemistry  were  found  in  Regions  I 
and  III  of  the  Piney  Point-Nanjemoy  aquifer.  A 
poor  match  was  found  for  Region  II.  This  analysis 
suggests  that  the  chemical  processes  which  control 
the  water  chemistry  of  the  Piney  Point-Nanjemoy 
aquifer  are  similar,  but  not  identical  to  the  processes 
which  control  the  water  chemistry  of  the  Aquia 
aquifer. 

Calcite  Cementation 

There  is  a  significant  difference  in  the  patterns  of 
calcite  cementation  between  the  Piney  Point- 
Nanjemoy  aquifer  and  the  Aquia  aquifer.  In  the 
Aquia  aquifer,  calcite  cementation  is  present  near 
the  outcrop  area,  but  generally  absent  in  down- 
gradient  areas.  In  contrast,  calcite  cementation  is 
present  in  both  upgradient  and  downgradient  areas 
of  the  Piney  Point-Nanjemoy  aquifer.  This  pattern 
of  calcite  cementation  is  consistent  with  the  equi¬ 
librium  chemistry  of  the  Piney  Point-Nanjemoy 
aquifer.  Figure  55  shows  calculated  saturation  in¬ 
dexes  of  Piney  Point-Nanjemoy  water  analyses 
(table  10).  Figure  55  shows  that  the  water  is  super¬ 
saturated  with  respect  to  calcite  along  the  entire 
length  of  the  aquifer. 

It  is  likely  that  the  supersaturation  of  Piney 
Point-Nanjemoy  water  has  resulted  from  dissolu¬ 
tion  of  magnesium  calcite  and  aragonite.  Both  of 
these  phases  are  less  stable  than  stoichiometric  cal¬ 
cite  in  aqueous  systems.  It  is  possible  that  the  disso¬ 
lution  of  magnesium  calcite  and  aragonite  and 
simultaneous  precipitation  of  calcite  is  a  continuing 
process  in  the  Piney  Point-Nanjemoy  aquifer.  Wol- 
last  (1971)  has  proposed  such  a  mechanism  on  theo¬ 
retical  grounds.  Wollast  (1971)  writes: 

“If  we  accept  that  aragonite  is  more  soluble  than  cal¬ 
cite,  then  we  must  expect  to  see  a  gradual  dissolution 
of  the  aragonite  and  a  preference  for  the  nucleation  of 
calcite.  .  . ” 

A  process  such  as  that  proposed  by  Wollast  (1971) 
could  possibly  explain  the  lack  of  correlation  be¬ 
tween  the  water  chemistry  of  Region  II  and  model 
II.  It  should  also  be  noted  that  this  process  is  not 
inconsistent  with  the  good  match  of  water  chemis¬ 
try  in  Region  III  with  model  III.  Model  III  requires 
the  dissolution  of  calcium  carbonate.  This  calcium 
carbonate  can  be  calcite  or  aragonite.  The  same 
argument  applies  to  Region  I. 


Figure  54.— Sodium  concentration  versus 

bicarbonate  concentration  for  water 
analyses  in  Region  III  of  the  Piney 
Point-Nanjemoy  aquifer. 


Log  a  SiQ, 


Figure  56.— Composition  of  Piney  Point-Nanjemoy 
water  plotted  as  a  function  of 
Na  + ,  H  + ,  and  Si02  in  the  system 
Nao0  -  A1,0, -Si0o-  H,0. 


Dissolved  Silica 

solved  silica.  Piney  Point-Nanjemoy  water  typically 
An  interesting  feature  of  Piney  Point-Nanjemoy  contains  45  to  65  mg/L  silica.  In  contrast,  Aquia 

water  is  the  relatively  high  concentration  of  dis-  water  commonly  contains  only  10  to  15  mg/L  silica. 
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Figure  55.— Calcite  saturation  indices  calculated  from  the  chemical  analyses  of  the  Piney  Point-Nanjemoy  water. 
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This  high  silica  content  in  Piney  Point  water  has 
been  documented  by  Mack  and  others  (1971).  These 
authors  note  that  the  Piney  Point-Nanjemoy  is  over- 
lain  by  the  highly  diatomaceous  Calvert  Formation. 
They  hypothesize  that  dissolution  of  the  siliceous 
diatoms  and  vertical  leakage  of  Calvert  water  into 
the  Piney  Point-Nanjemoy  causes  the  high  silica 
content.  One  problem  with  this  hypothesis  is  that 
water  in  equilibrium  with  amorphous  silica  should 
have  silica  concentrations  much  higher  than  those 
observed  in  the  Piney  Point  aquifer  (Garrels,  1975). 
Assuming  that  the  source  of  the  silica  is  dissolution 
of  diatoms  in  the  Calvert  Formation,  an  additional 
mechanism  must  be  limiting  silica  concentrations  to 
the  observed  45-65  mg/L.  It  is  possible  that  simple 
dilution  is  the  limiting  mechanism.  It  is  also  possi¬ 
ble  that  reaction  of  dissolved  silica  with  the  silicate 
minerals  in  the  aquifer  may  be  the  limiting  factor. 

It  is  not  possible  to  test  the  dilution  mechanism 
hypothesis  easily.  It  is  possible,  however,  to  test  the 
silicate  reaction  hypothesis.  Garrels  (1975)  has  cal¬ 
culated  a  phase  diagram  on  which  water  can  be 
plotted  as  a  function  of  dissolved  sodium,  hydrogen- 
ion,  and  silica  in  the  system  Na20— A1203— Si02— H20. 
With  this  diagram,  Garrels  (1975)  has  demonstrated 
that  most  dilute  ground  and  surface  water  falls  in 
the  kaolinite  field  of  stability.  Garrels  (1975)  also 


shows  that  water  tends  to  evolve  along  the  kaolinite- 
montmorillonite  stability  boundary.  This  effectively 
limits  dissolved  silica  concentrations  to  less  than  70 
mg/L.  If  a  similar  mechanism  is  operating  on  silica 
concentrations  of  Piney  Point-Nanjemoy  water,  a 
similar  pattern  should  develop.  Figure  56  shows 
available  analyses  of  Piney  Point-Nanjemoy  water 
plotted  as  a  function  of  dissolved  sodium,  hydrogen- 
ion,  and  silica  in  the  system  Na20— A1203— Si02— H20. 
The  stability  diagram  used  is  after  Garrels  (1975). 
This  figure  shows  clearly  that  the  composition  of 
water  is  indeed  evolving  along  the  kaolinite-mont- 
morillonite  boundary.  Water  from  Region  I  of  the 
Piney  Point-Nanjemoy  aquifer  plots  on  the  lower 
part  of  this  diagram.  Water  from  Region  III  plots 
on  the  upper  part,  and  water  from  Region  II  is  inter¬ 
mediate.  This  trend  is  consistent  with  the  observa¬ 
tion  that  sodium  concentrations  increase  along  the 
flowpath  and  hydrogen-ion  concentrations  decrease. 
The  significant  feature  of  this  diagram  is  that  it 
shows  the  silica  concentrations  adjust  to  the  chang¬ 
ing  log  a  Na+/a  H+  ratio  so  that  the  water  composi¬ 
tions  remain  in  the  kaolinite  stability  field.  This 
strongly  suggests  that  the  stability  relationship  of 
kaolinite  and  montmorillonite  is  the  controlling 
mechanism  which  limits  the  concentration  of  dis¬ 
solved  silica  in  Piney  Point-Nanjemoy  water. 


SUMMARY  AND  CONCLUSIONS 


The  Aquia  Formation  is  a  fine-  to  medium-grained 
glauconitic  quartz  sand  that  has  abundant  shell  ma¬ 
terial.  The  Aquia  Formation  serves  as  a  fresh-water 
aquifer  in  most  of  Southern  Maryland.  The  Aquia 
outcrops  in  a  band  paralleling  the  Fall  Line.  Re¬ 
charge  to  the  Aquia  aquifer  occurs  as  leakage  from 
the  Pleistocene  water-table  aquifer  and  leakage  from 
adjacent  confining  beds.  The  transmissivity  of  the 
Aquia  aquifer  in  Southern  Maryland  is  highest  in 
Anne  Arundel  County.  The  transmissivity  decreases 
to  the  south  and  is  not  usable  as  an  aquifer  in  south¬ 
ern  St.  Marys  County. 

The  upper  sands  of  the  Nanjemoy  Formation  and 
the  sands  of  the  Piney  Point  Formation  in  Southern 
Maryland  are  hydraulically  connected  and  function 
as  one  hydraulic  unit.  The  Piney  Point-Nanjemoy 
aquifer  is  a  fine-  to  coarse-grained  glauconitic  quartz 
sand.  The  sands  of  the  Piney  Point-Nanjemoy  aqui¬ 
fer  are  thickest  in  St.  Marys  County  and  thin  the 
northeastern  direction.  The  sands  are  truncated  in 
the  subsurface  of  northern  Calvert  and  southern 
Prince  Georges  Counties  as  a  result  of  the  Oligocene 
withdrawal  of  the  sea.  The  Piney  Point-Nanjemoy 
aquifer  is  recharged  from  the  Pleistocene  water- 
table  aquifer  by  leakage  through  the  Chesapeake 


Group  confining  bed.  The  highest  transmissivity  of 
the  Piney  Point-Nanjemoy  aquifer  occurs  in  St. 
Marys  County. 

The  subsurface  continuity  of  the  Aquia  and  Piney 
Point-Nanjemoy  aquifers  and  their  confining  beds 
has  been  modified  by  Pleistocene  erosion.  The  Aquia 
aquifer  is  completely  truncated  in  the  Chesapeake 
Bay  near  Kent  Island.  Similarly,  the  Chesapeake 
Group  sediments  are  completely  truncated  near 
Benedict,  Md.  These  Pleistocene  erosional  channels 
are  a  major  control  on  the  head  distribution  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifers.  The 
truncation  of  aquifers  and  confining  beds  by  Pleisto¬ 
cene  channels  increases  the  possibility  of  brackish- 
water  intrusion  in  some  places.  The  Aquia  aquifer 
near  Kent  Island  and  Anne  Arundel  County  in  par¬ 
ticular  may  be  susceptible  to  future  brackish  water- 
intrusion  if  heads  fall  below  sea  level. 

A  digital  finite-difference  quasi  three-dimensional 
ground-water  flow  model  was  constructed  of  the 
Aquia  and  Piney  Point-Nanjemoy  aquifer  system  in 
Southern  Maryland.  The  model  was  calibrated  by 
matching  computed  potentiometric  surfaces  against 
measured  potentiometric  surfaces  for  three  stages 
of  aquifer  development.  Calibration  was  achieved  by 
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varying  vertical  hydraulic  conductivities  of  confin¬ 
ing  beds  within  the  range  of  values  indicated  by  lab¬ 
oratory  tests  of  confining-bed  material.  When  com¬ 
puted  heads  at  each  stage  of  the  aquifer  system’s 
development  matched  measured  heads  at  the  speci¬ 
fied  times,  the  model  was  considered  calibrated. 

After  calibration,  the  model  was  used  to  predict 
future  water  levels  based  upon  possible  scenarios  of 
future  pumpage.  These  simulations  included: 

(1)  A  series  of  simulations  using  estimates  of 
future  pumpage  for  the  periods  1980-1985, 
1980-1990,  and  1980-2000.  This  series  of  simu¬ 
lations  predicts  that  by  the  year  2000  the  po- 
tentiometric  surface  of  the  Aquia  aquifer  near 
Lexington  Park  will  be  more  than  100  ft  below 
sea  level.  It  is  possible  that  these  lowered 
water  levels  will  make  water  production  from 
some  domestic  wells  impractical  near  Lexing¬ 
ton  Park  by  the  year  2000.  These  lowered 
water  levels  may  also  significantly  increase 
the  cost  of  pumping  from  existing  municipal 
wells.  Because  of  these  possibilities,  the  devel¬ 
opment  of  deeper  Cretaceous  aquifers  for 
water  production  might  be  considered  at  Lex¬ 
ington  Park. 

(2)  A  series  of  simulations  in  the  Lexington  Park- 
PNATS  area  to  predict  water  levels  in  the  pe¬ 
riod  1980-2000,  assuming  the  Aquia  aquifer  is 
pumped  1.5  Mgal/d,  2.0  Mgal/d,  2.5  Mgal/d, 
and  3.0  Mgal/d.  This  series  of  simulations  pro¬ 
vides  water  planners  with  drawdown  esti¬ 
mates  for  different  pumping  rates.  The  simu¬ 
lations  were  designed  to  help  planners  esti¬ 
mate  the  maximum  pumpage  rate  that  can  be 
sustained  by  the  Aquia  aquifer.  If,  for  exam¬ 
ple,  water  planners  decide  that  an  additional 
average  drawdown  of  more  than  60  ft  would 
be  undesirable,  the  yield  of  the  Aquia  aquifer 
in  this  area  would  be  limited  to  between  2.0 
and  2.5  Mgal/d. 

(3)  A  series  of  simulations  showing  the  effects  of 
1  Mgal/d  pumpage  from  the  Aquia  aquifer  in 
northern  St.  Marys  County.  One  simulation 
assumed  1  Mgal/d  pumpage  from  a  single 
METCOM  well  field.  Another  simulation  as¬ 
sumed  that  the  1  Mgal/d  pumpage  was  spread 
among  four  METCOM  well  fields.  This  simu¬ 
lation  was  made  to  illustrate  to  water  planners 
the  effects  of  spreading  pumpage  over  a  broad 
area  as  opposed  to  concentrating  it  in  one  area. 
This  simulation  shows  that  spreading  pump¬ 
age  results  in  drawdowns  over  a  larger  area, 
but  that  the  drawdowns  are  less.  Spreading 
pumpage  among  several  well  fields  would 
therefore  minimize  impacts  on  already  exist¬ 
ing  domestic  and  small  industrial  wells. 


(4)  A  series  of  simulations  in  the  Lexington  Park- 
PNATS  area  to  show  the  drawdown  effects  on 
the  Piney  Point-Nanjemoy  aquifer  due  to 
pumpage  from  the  Aquia.  These  simulations 
indicate  that  90  ft  of  drawdown  in  the  Aquia 
aquifer  over  10  years  results  in  4  ft  of  draw¬ 
down  in  the  Piney  Point-Nanjemoy  aquifer. 
This  result  shows  that  Aquia  pumpage  will 
have  drawdown  effects  on  the  Piney  Point- 
Nanjemoy,  but  these  effects  are  not  large. 

(5)  A  series  of  simulations  to  show  the  effects  of 
placing  a  0.5  Mgal/d  well  field  in  southern 
Anne  Arundel  County  near  the  towns  of 
Shady  Side,  Bristol,  and  Fairhaven.  This  simu¬ 
lation  indicates  that  pumpage  of  0.5  Mgal/d 
near  any  one  of  these  towns  will  have  a  moder¬ 
ate  impact  on  water  levels. 

(6)  A  simulation  of  0.5  Mgal/d  from  the  Aquia  for 
the  Chalk  Point  Power  Plant.  This  simulation 
shows  that  about  45  ft  of  drawdown  develops 
after  10  years  of  pumping.  This  relatively 
large  impact  is  partly  due  to  the  low  transmis¬ 
sivity  of  the  Aquia  at  Chalk  Point,  and  partly 
due  to  domestic  and  small  industrial  pumpage 
in  the  same  area. 

(7)  A  hypothetical  simulation  that  assumed  all 
GAP  holders  in  the  Aquia  and  Piney  Point- 
Nanjemoy  aquifers  pumped  their  maximum 
water  allocations  for  a  period  of  10  years. 
Large  drawdowns  occur  in  the  Aquia  aquifer 
near  Lexington  Park,  Leonardtown,  Piney 
Point,  the  Calvert  Industrial  Park,  and  Cove 
Point.  The  drawdown  effects  of  this  simula¬ 
tion  were  not  as  large  on  the  Piney  Point- 
Nanjemoy  aquifer  because  most  GAP  ground- 
water  users  in  Southern  Maryland  pump  from 
the  Aquia  aquifer.  Large  drawdowns  are  also 
predicted  to  occur  in  the  Aquia  aquifer  near 
the  planned  community  of  Prospect  Planta¬ 
tion  in  Queen  Anne’s  County.  The  cone  of 
depression  in  this  area  is  predicted  to  spread 
to  the  subsurface  truncation  of  the  Aquia 
aquifer  under  the  Chesapeake  Bay.  This  could 
possibly  result  in  brackish  water  being  drawn 
into  the  Aquia  at  this  location.  This  simula¬ 
tion  suggests  that  near  Lexington  Park,  Leon¬ 
ardtown,  and  Piney  Point  in  St.  Mary’s  Coun¬ 
ty  appropriation  requests  for  additional  Aquia 
ground  water  should  be  carefully  reviewed. 
The  Piney  Point-Nanjemoy  aquifer  does  not 
appear  to  be  nearing  the  limits  of  production 
capability. 

(8)  A  simulation  of  0.5  Mgal/d  from  the  Aquia 
aquifer  near  Chesapeake  Beach.  This  simula¬ 
tion  results  in  about  20  ft  of  drawdown  occur¬ 
ring  over  5  years. 
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Water  produced  from  the  Aquia  aquifer  generally 
has  good  chemical  quality.  Near  the  outcrop/subcrop 
area  of  the  Aquia  the  water  is  a  calcium  and  mag¬ 
nesium  bicarbonate  type  water  with  relatively  low 
(6. 5-7. 5)  pH.  Downgradient  from  the  outcrop/sub¬ 
crop  area,  Aquia  water  changes  to  a  sodium  bicar¬ 
bonate  type  water  with  relatively  high  (7. 5-8. 5)  pH. 
It  has  been  shown  that  the  Aquia  aquifer  can  be  di¬ 
vided  into  three  regions  based  on  water  chemistry. 
Chemical  models  which  predict  water  composition 
in  each  region  of  the  aquifer  have  been  developed. 
These  models  were  tested  by  comparing  predicted 
changes  of  water  chemistry  to  measured  changes  of 
water  chemistry  for  three  regions  of  the  Aquia  aqui¬ 
fer.  The  matches  obtained  by  these  comparisons 
were  good  in  all  three  cases  and  this  supports  the  ap¬ 
plicability  of  the  models.  Mass  balance  calculations 


with  these  models  suggest  that  the  observed  calcite 
cementation  of  the  Aquia  aquifer  is  post-depositional 
and  has  resulted  from  reaction  of  ground  water  with 
aquifer  material. 

Water  produced  from  the  Piney  Point-Nanjemoy 
aquifer  also  has  generally  good  chemical  quality.  In 
upgradient  areas,  the  water  is  a  calcium  and  mag¬ 
nesium  bicarbonate  type  water  with  relatively  low 
(6. 0-7.0)  pH.  Downgradient,  Piney  Point-Nanjemoy 
water  is  a  sodium  bicarbonate  type  water  with  rela¬ 
tively  high  (7. 0-8. 2)  pH.  The  Piney  Point-Nanjemoy 
aquifer  can  be  divided  into  three  regions  based  on 
water  chemistry  in  a  similar  manner  as  the  Aquia 
aquifer.  The  chemistry  of  these  regions  can  be  de¬ 
scribed  by  the  same  chemical  models  which  were 
developed  for  the  regions  of  the  Aquia  aquifer. 
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Plate  1.— Cross  section  A-A'  showing  major  hy 
drologic  units  in  Southern  Maryland 
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Plate  2. — Cross  section  B-B '  showing  major  hy¬ 
drologic  units  in  Southern  Maryland 
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Plate  3.— Map  showing  the  approximate  pre¬ 
pumping  potentiometric  surface  of  the 
Aquia  aquifer 
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Plate  4.— Map  showing  the  potentiometric  sur¬ 
face  of  the  Aquia  aquifer,  1952 
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Plate  5.— Map  showing  the  potentiometric  sur¬ 
face  of  the  Aquia  aquifer,  May  19-23, 
1980 
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Plate  6.— Map  showing  the  approximate  pre¬ 
pumping  potentiometric  surface  of  the 
Piney  Point-Nanjemoy  aquifer 
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Plate  7.— Map  showing  the  potentiometric  sur 
face  of  the  Piney  Point-Nanjemoy  aqui 
fer,  1952 
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Plate  8.— Map  showing  the  potentiometric  sur¬ 
face  of  the  Piney  Point-Nanjemoy  aqui¬ 
fer,  May  19-23, 1980 
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Plate  8.— Map  showing  the  potentiometric  sur 
face  of  the  Piney  Point-Nanjemoy  aqui 
fer,  May  19-23, 1980 


Plate  9.— Map  showing  the  potentiometric  sur¬ 
face  of  the  water-table  aquifer  of  South¬ 
ern  Maryland 
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Plate  10.— Map  showing  simulated  prepumping, 
1952,  and  1980  potentiometric  surfaces 
of  the  Aquia  aquifer  superimposed  on 
measured  water  levels 
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Plate  10.— Map  showing  simulated  prepumping, 
1952,  and  1980  potentiometric  surfaces 
of  the  Aquia  aquifer  superimposed  on 
measured  water  levels 
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Simulated  pre- pumping  potentiometric  surface  and  water  levels  measured  early  in  aquifer 
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Plate  11.— Map  showing  simulated  prepumping, 
1952,  and  1980  potentiometric  surfaces 
of  the  Piney  Point-Nanjemoy  aquifer 
superimposed  on  measured  water 
levels 


Platell.— Map  showing  simulated  prepumping, 
1952,  and  1980  potentiometric  surfaces 
of  the  Piney  Point-Nanjemoy  aquifer 
superimposed  on  measured  water 
levels 
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Simulated  potentiometric  surface  of  1952  and  water  levels  measured  approximately 
1952  by  Otton  (  1955). 


Simulated  1980  potentiometric  surface  and  water  levels  measured  May  19-23, 1980. 
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Plate  12. — Locations  of  selected  observation 
wells  in  the  Aquia  aquifer 
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•  in  the  Aquia  aquifer.  Number  is  the 
U.S.G.S.  well  number. 
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Plate  13.— Locations  of  selected  observation 
wells  in  the  Piney  Point-Nanjemoy 
aquifer 
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